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ABSTRACT 


A  detailed  program  of  analytical  and  experimental  investigation  into  the 
technical  feasibility  of  the  H-Rho  nuclear  reactor  control  system  and  tetra 
phosphorous  trisulfide  (P^S^)  as  a  reactor  coolant  are  presented  in  this  study. 

The  H-Rho  control  system  investigation  includes  determinations  of  hydrogen 
transport  (permeability)  through  porous  bodies  of  zirconium  hydride  and 
uranium  zirconium  hydride  as  well  as  the  absorption-desorption  of  hydrogen 
within  the  hydride  bodies.  The  effects  of  temperature,  pressure,  particle 
size,  inert  gas  and  irradiation  are  evaluated. 

The  investigation  of  the  coolant  includes  determinations  of  viscosity, 
surface  tension,  thermal  conductivity,  vapor  pressure,  liquid  density  and 
heat  of  fusion.  The  effects  of  irradiation  on  the  P^S^  are  evaluated. 

The  results  of  the  investigation  demonstrate  that  both  the  H-Rho  con¬ 
trol  system  and  the  P^S^  coolant  are  particularly  suited  for  application  to 

the  Terrestrial  Unattended  Reactor  Power  System  (1  URPS). 

(Distribution  Limitation  Statement  No.  2) 
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SECTION  I 

GENERAL  INTRODUCTION 

The  results  of  an  investigation  into  the  technical  feasibility  of  the  control 
mechanism  (H-Rho)  and  into  the  properties  of  the  coolant  (Theon)  selected 
for  the  Terrestrial  Unattended  Reactor  Power  System  (TURPS)  are  presented 
in  this  report.  Portions  of  the  investigations  reported  herein  were  performed  at 
the  Nuclear  Engineering  Test  Facility  at  W  right -Pattersor  Air  Force  Base,  Ohio 
and  at  the  Battelle  Memorial  Institute,  Columbus,  Ohio. 

Under  the  TURPS  program  there  have  been  two  pre-  ious  studies  as  reported 
in  the  following  documents: 

AFWL-TR-64-168,  "Terrestrial  Direct  Conversion  Nuclear  Power  Plant 

Conceptual  Design.  "  March  1965 

AFWL-TR-66-16,  "Feasibility  Study  of  Terrestrial  Thermoelectric 

Nuclear  Power  Plant.  "  June  1966 

Under  the  first  study,  TURPS  was  designed  to  performance  specifications 
developed  by  the  Air  Force  Weapons  Laboratory.  Under  the  second  study  two 
major  constituents  of  TURPS,  H-Rho  control  and  Theon  coolant, were  investi¬ 
gated  through  out-of-core  experiment  and  consequent  analysis.  The  current 
study  initiated  in-core  experimental  investigation  and  continued  the  out  of 
core  work  with  both  H-Rho  and  Theon. 

The  H-Rho  tasks  include  the  performance  and  analysis  of  a  number  of  tests 
designed  to  investigate  the  hydrogen  mass  transport  characteristics  of  porous 
zirconium  and  uranium -zirconium  hydride  bodies.  A  neutron  irradiation  pro¬ 
gram  was  performed  using  similar  uranium -zirconium  bodies  to  determine 
fission  product  retention.  Also  performed  was  a  series  of  thermal  diffusivity 
and  specific  heat  measurements  on  similar  bodies.  From  these  data  the  thermal 
conductivities  of  the  bodies  were  determined. 

In  parallel  with  the  experimental  efforts,  a  major  revision  was  made  in  the 
computer  program  for  the  analysis  of  H-Rho  experimental  data.  The  new  pro¬ 
gram,  HYDRA  II,  retains  the  basic  mathematical  model  in  which  the  body  is  di 
vided  into  a  number  of  discrete  zones  normal  to  the  direction  of  hydrogen  trans¬ 
port.  In  HYDRA  I,  the  analytic  procedure  was  one  of  stepping  sequentially  from 
zone  to  zone  through  the  body.  HYDRA  II  solves  the  equations  for  ali  zones  simul¬ 
taneously. 
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Under  the  Theon  tasks  a  number  of  measurements  were  made  to  deter¬ 
mine  physical  properties  such  as  viscosity,  surface  tension,  vapor  pressure 
and  thermal  conductivity.  A  300-hour,  high  temperature  irradiation  of  the 
coolant  has  been  completed.  An  experimental  program  to  determine  the  com¬ 
patibility  of  Theon  with  various  mefals  and  metal  surface  preparations  was 
completed.  A  dynamic  corrosion  loop  has  oeen  constructed  for  performance 
of  a  4500-hour,  high  temperature  study  of  the  Zircaloy-Theon  system. 

The  following  sections  of  this  report  present  the  technical  approach  and  re¬ 
sults  of  these  investigations.  These  are  divided  into  the  following  topics: 

Section  II  Hydrogen  Transport  in  a  Zirconium  Hydride  Matrix 

Section  HI  Thermal  Conductivity  of  Porous  Zirconium  Hydride 

Section  IV  Physical  Properties  of  Theon  (P4S^) 

Section  V  Irradiation  of  Theon  (P.S_) 

4  d 

Section  VI  Chemical  Compatibility  of  Theon  with  Metals 

For  reasons  of  National  Security,  Volume  II  presents  the  results 
of  the  irradiation  of  zirconium  and  uranium -zirconium  hydride. 
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SECTION  II 

HYDROGEN  TRANSPORT  IN  A  ZIRCONIUM  HYDRIDE  MATRIX 

1 .  Purpose 

This  task  experimentally  determined  the  hydrogen  response  character¬ 
istics  of  the  proposed  TURPS  fuel  form.  As  presently  designed,  the  TURPS 
core  and  hydrogen  reservoir  region  will  consist  of  porous  bodies  of  uranium 
zirconium  hydride  and  zirconium  hydride,  respectively.  Typically,  the 
bodies  are  fabricated  from  mesh-size  particles  using  powder  metallurgy  tech¬ 
niques.  During  the  present  work,  a  number  of  sample  bodies  were  fabricated 
from  spherical  zirconium,  nonspherical  zirconium, and  spherical  zirconium- 10 
percent  uranium  particles.  Selected  bodies  were  subjected  to  the  following  sets 
of  tests: 

(1)  Permeability  measurements. 

(2)  Hydrogen  absorption  transients. 

(3)  Thermal -hydrogen  cycling. 

In  the  performance  of  this  task,  several  major  activities  were  required: 

(1)  Modify  the  existing  HYDRA  code  to  increase  its  analytical  capability 
and  improve  its  operating  characteristics  for  use  as  a  data  analysis  tool. 

(2)  Plan  a  test  program  to  study  the  dynamic  hydrogen  transport  process. 

(3)  Design,  build, and  check  out  an  apparatus  for  investigating  the  hydrogen 
absorption  characteristics  of  porous  zirconium  hydride  and  uranium-zirconium 
hydride  bodies.  Prepare  the  supporting  test  procedures. 

(4)  Perform  the  various  experiments;  reduce  and  analyze  the  data. 

2.  HYDRA  II  Digital  Code 

HYDRA  II  is  a  digital  computer  code  written  to  simulate  the  dynamic  process 
of  hydrogen  transport  from  the  gas  phase  into  a  hydride  bed.  The  basic  mathe¬ 
matical  formulation  is  the  same  as  that  of  the  earlier  computer  code  HYDRA  I 
with  the  inclusion  of  the  continuity  equation  to  describe  the  presence  of  an  inert 
gas.  Further  modifications  included  the  use  of  the  Klinkenberg  equation  in  cal¬ 
culating  the  permeability  and  the  flexibility  of  the  code  to  divide  the  bed  into  as 
many  as  10  zones.  The  nomenclature  used  in  the  development  of  the  hydrogen 
transport  analysis  and  subsequent  experimental  program  i»  presented  in  table  I. 
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Table  I 

NOTATIONS  (USED  IN  HYDROGEN  TRANSPORT  SECTION) 
a  Effective  radius  of  particle 

c  Average  hydrogen  concentration 
c  Solid  heat  capacity 

B 

D  Diffusion  coefficient 
d  Permeability 
AH  Heat  of  reaction 
K  „*.*/.* 

L  Volumetric  heat  loss  rate 
P  Total  gas  pressure 
Pjj  Hydrogen  gas  pressure 
Pj  Inert  gas  pressure 

S  Volumetric  heat  generation  or  absorption  rate 
t 

T  $Ddt 
0 

t  Time 

V  Superficial  gas  velocity 
x  Length 

o.  Constant  used  in  equilibrium  equation 
Constant  used  in  equilibrium  equation 
<  Void  fraction 
0  Temperature 
M  Gas  viscosity 
Pjj  Hydrogen  density 
pf  Inert  gas  density 
Pg  Solid  density 

Surface  concentration 
X  Dummy  variable 


a.  Mathematical  Treatment  and  Solution 

The  basic  mathematical  solution  is  obtained  by  simultaneously  solving  the 
following  set  of  equations. 


i  ■  *¥  I  j*(Tn>  '  K' 
a  n  «  i  ( 


-KT  Jn 


*(X)eKX  dA 
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n 

log  PH  =  ^  (ai  +  ^./e)^*1 


i  =  1 


a  p 


H 


dc  .  *«>HV) 


•if  +  (1  '  €)Hf  +  ax 


=  o 


apT  a(p.v) 

«  +  -i —  = 

at  ax 


v  =  -d  8P 

m  ex 


de  AH  dc  1  ip  ,  I 

dt  p  c  dt  p  c  (1-0  /  M 
^s  s  Ks  s 


(2) 

(3) 

(4) 
(3) 
(6) 


The  solution  of  these  equations  is  performed  in  the  following  manner.  The 
integral  term  in  equation  (1)  is  integrated  by  parts  to  yield 


_  00 
dc.6D  V 

dt  '  2  L 


rp  T 
T  n 


+  (Tn)  -  e 


n*  1 


,  KTn  KX\  n  -  f  d*  KA  A 

■  (4,U)e  Jo  )  dXe  dAj 


6D 

2 


00 

I 


-KT  /»**  .  ... 

n  \  d+  KA 

e  i  dAe  dA 


(7) 


n  =  I  0 

as  vKO)  =  0  for  derivation  and  validity  of  equation  (1). 

The  assumpt’on  is  now  made  that  is  linear  in  each  time  step  so  that 


_  ^i  ^i-1 
dA  '  A  -  A. 

i  i-i 

in  the  ith  time  step,  i  =  1  at  time  zero. 

Substitution  of  equation  (8)  into  (7)  yields 


(8) 
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which  upon  reversing  summation  signs  yields 


d£  .  6D  V  1*1  -  VA  t 
dt  a2  AVV  \-i)  L 


■K(T  -  T.)  -K(T  -  T.  .) 
n  1  n  l-l 


As  T  >  T.  for  all  time  steps,  the  infinite  series 
n —  l 


'«■  1 


-K(T  -  T.) 
n  i 


converges  for  all  time  steps. 

2  2 

As  K  *  f  the  series  to  be  evaluated  has  tho  form 

a 


n  *  i  n 

2 

where  a  >  °. 

The  method  used  in  evaluating  this  series  is  given  in  Appendix  I. 
Setting 

2  2 

-  -  (T  -  T  ) 
m  t*  n  1 


•f  €*  II  1 

c  .  t  e _ 1 _ 

1  A  n2 

n  *  1 

and  substituting  in  equation  (9)  yields 

*' ■  7  i  f^r)  M 

Equation  (’0)  is  then  substituted  into  equation  (3)  and  equation  (5)  into  equa 
tions  (3)  and  (4)  to  give 


6 
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Finite  difference  equations  are  written  for  equations  (11)  and  .12)  in  each  of 
the  zones  into  which  the  bed  has  been  divided.  If  there  a»*e  N  zones,  then  equa¬ 
tions  (11)  and  (12)  yield  2N  equations  and  equation  (2)  yields  N  equations  in  the 
3N  unknowns,  PT,  P,,,  ib in  each  zone  of  the  bed.  The  3N  equations  in  3N  un- 
knowns  are  solved  by  the  Newton -Rhapson  method  which  converges  rapidly.  As 
u,  0,  etc.  ,  are  functionally  related  to  Pp  P^,  4*.  an  outer  iterative  loop  is  |«r- 
formed  on  equation  (fi)  until  the  complete  set  of  equations  (')  to  (6)  converge. 
Usually  about  three  or  four  iterations  ai e  required  to  converge  on  equations  (11), 
(12),  and  (2),  and  four  or  five  iterations  to  converge  on  the  set  of  equations 
(1)  to  (6).  The  process  is  quite  stable  and  allows  use  of  large  time  increments. 

Convergence  is  more  critical  initially,  however,  and  smaller  time  step^ 
must  be  used.  The  time  step  is  adjusted  a  itomatieally  in  the  program  to  in¬ 
crease  accuracy  and  decrease  computer  running  time. 

b.  Computer  Program 
(1)  Input 

The  necessary  input  data  are  given  in  the  order  in  which  they  are  read. 

NZONES  Number  of  zones  that  the  model  is  to  simulate 

NTAB  Number  of  pairs  of  table  entries  of  time  versus  total  pressure 

INPT  Option  (1):  If  INPT=0,  the  standard  output  is  obtained 

Option  (2):  If  ENPTX),  all  intermediate  calculations  are  printed  out. 

This  is  useful  in  determining  if  the  initial  time  increment 
is  too  large . 


XA 

„  .  ,  .  .  -B/R0 

Coefficient  A  in  Ae 

XB 

Coefficient  B  in  Ae 

DP 

Average  particle  diameter 

FD 

Shape  factor  for  diffusion 

DEPTH 

Length  of  bed 

EPS 

Void  fraction 

FF 

Shape  factor  for  flow 

TM  AX 

Maximum  time  the  transient  is  to  be  computed  (seconds) 

DELTAT 

Initial  time  increment 

TA.MB 

Ambient  temperature  (°Fi 

TZERO 

Initial  bed  temperature  (*F) 
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Fi  Fraction  of  inert  gas  present 

RHOS  p  ,  density  of  solid 

s 

CESS  c  ,  heat  capacity  of  solid 

s 

AREA  Cross  sectional  area  of  bed 

XKL  Heat  loss  coefficient 

Ef!3N(J),  J  =  l,  NZONES  S  ,  volumetric  heat  generation  rate 

11  O 

SIGI  Oj,  collision  diameter  of  inert  (A) 

EPSKI  t €  /K)p  Lennaid-Jones  energy  characteristic  of  inert  gas  (°K) 

XMI  Atomic  weight  of  inert  gas 

(TP(J),  Q(J))J  =  1,  NTAB  time  versus  pressure  table 

PIV1  Initial  pressure  of  inert  gas 

PIV2  Pressure  of  inert  gas  at  first  time  step 

QEX  Mixing  factor 

ROMANY  Exhaust  volume 

UMAX  Maximum  time  increment 

DMIN  Minimum  time  increment 

BE  KlinKenberg  constant 

(2)  Output 

The  printed  output  includes  the  temperature,  hydrogen  pressure,  partial 
pressure  caused  from  the  inert  gas,  total  pressure  and  average  hydrogen  par¬ 
ticle  concentration  in  each  zone  as  well  as  the  total  amount  of  absorbed  hydrogen. 

The  FORTRAN  listing  of  HYDRA  II  is  provided  in  Appendix  II. 


jeriment  Planning  Overview 


a-  Purpose 


The  proposed  tests  are  being  performed  to  determine  the  hydrogen  absorp¬ 
tion  rates  of  porous  zirconium  hydride  and  uranium  zirconium  hydride  bodies 
and  to  determine  the  effect  of  thermal  and  hydrogen  cycling  on  these  rates. 


Porous  bodies  were  prepared  from  spherical  zirconium  hydride,  using  a 
number  of  particle  size -void  fraction  combinations,  to  illustrate  the  effect  of 
body  configuration  on  hydrogen  response.  Additionally,  bodies  were  prepared 
from  nonsplu  rical  particles  to  study  the  problems  associated  with  use  of  the 
less  expensive  material.  Finally.,  uranium  zirconium  hydride  spheres  were 
used  so  that  any  variation  in  hydrogen  absorption  characteristics  resulting 
trom  the  presence  of  the  fuel  could  be  observed. 


_  ^  .  .  | ,  -  ..  v-  »r"  o>.  -  -w  ***«*  •  fc  ->  *■  - 
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During  the  previous  experimental  effort  (AFWL-TR-66-16),  absorption 
characteristics  of  relatively  deep  beds  (three  to  four  inches)  of  zirconium 
hydride  powders  were  measured.  The  rates  of  absorption  were  found  to  be 
limited  by  a  pressure  drop  within  the  beds.  In  an  effort  to  observe  both  pres¬ 
sure  drop  and  diffusion  limited  situations,  sintered  bodies  having  a  nominal 
length  of  one  centimeter  were  chosen  for  the  present  studies.  Experimental 
conditions  and  body  characteristics  were  selected  to  allow  study  of  a  large 
range  of  absorption  rates. 

b.  Nature  of  Tests 

The  hydrogen  absorption  test  subjected  a  body  to  a  step -like  change  in  hydro¬ 
gen  pressure.  The  absorption  rate  resulting  from  this  change  was  measured  by 
observing  the  pressure  decay  of  a  known  volume  of  hydrogen  gas.  Experi¬ 
mental  pressures  and  temperatures  were  varied  parametrically.  Selected 
runs  were  repeated  in  the  presence  of  various  amounts  of  inert  gas,  simu 
lating  the  effect  of  possible  gaseous  fission  product  release  in  an  operating 
reactor. 

To  calculate  theoretical  absorption  rates,  the  permeability  of  the  body  to 
hydrogen  flow  at  various  pressure  levels  and  Reynolds  numbers  must  be  known. 
Therefore,  a  series  of  permeability  tests  were  run  on  each  body.  Pressure 
drop,  pressure  level,and  flow  rates  were  measured  over  a  wide  range  of  ex¬ 
perimental  conditions. 

Temperature -hydrogen  cycling  tests  were  conducted  on  a  number  of  selected 
samples.  The  temperature  cycling  induces  a  hydrogen  pressure  cycling  of  the 
bodies.  Cycling  limits  were  set  by  the  range  of  hydrogen  content  variation 
anticipated  in  the  TURPS  core  during  a  typical  load  change  transient. 

c.  Expected  Results 

The  results  of  the  absorption  tests  were  planned  to  provide  the  experimental 
data  needed  to  verify  the  present  analytical  model  used  to  describe  the  hydrogen 
absorption  process.  This  model  is  evaluated  by  using  the  HYDRA  II  code.  The  ex¬ 
perimental  results  provide  more  realistic  body  permeabilities  for  use  in  the  code, 
and  indicate  the  validity  of  applying  literature  values  of  the  diffusion  coefficient. 
The  net  result  is  the  evolution  of  an  accurate  and  reliable  method  of  analyzing 
the  hydrogen  absorption  process,  verified  by  experiments  over  the  range  of 
parameters  of  interest  in  the  TURPS  design. 
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4.  Experimental  Program 

a.  Experimental  Conditions 

The  conditions  varied  in  the  absorption  tests  were  initial  body  temperature, 
pressure  step  and  gas  composition.  Because  a  complete  parametric  study  using 
all  of  the  bodies  was  outside  of  the  scope  of  the  present  program,  only  two  of  the 
bodies  were  so  tested.  The  others  were  tested  at  fewer  conditions. 

For  convenience  in  identifying  the  various  sintered  bodies  refer  to  table  II. 
The  identification  numbers  are  used  in  later  tabulations  of  the  experimental  con¬ 
ditions. 

Table  II 


CHARACTERIZATION  OF  POROUS  BODIES 


Material 

Particles 

Type  of  Test 

Number  of 
Bodies 

Identification 

ZrHx 

Spherical 

Absorption 

5 

SI,  S2,  S3,  S4,  S5 

ZrH 

x 

Spherical 

Cycling 

2 

S6,  S7 

U0.044ZrIix 

Spherical 

Cycling 

3 

FI,  F2,  F3 

ZrHx 

Nonspherical 

Absorption 

2 

Nl,  N2 

ZrHx 

Nonspherical 

Cycling 

1 

N3 

Table  III  shows  the  range  of  experimental  conditions  covered  in  the  tests. 

Table  HI 

NOMINAL  RANGE  OF  EXPERIMENTAL  CONDITIONS 


Type  of  Test 

- r 

Condition  1 

Ranee 

Absorption 

Initial  body  temperature  (8F) 

100  to  1500 

Pressure  step  (atm) 

0. 1  to  1. 1;  0.  5  to  1. 5; 

1.0  to  2.0 

Permeability 

Pressure  level  (psi) 

1 . 0  to  50.  0 

Reynolds  number 

0.0001  to  0.01 

Cycling 

Temperature  cycle  amplitude  (°F) 

200 

Atom  ratio  cycle  amplitude 

0.  10 

b.  Data  Output 


The  raw  data  generated  during  the  absorption  tests  included  the  initial  tem¬ 
perature  and  hydrogen  dissociation  pressure  of  the  body.  From  these  conditions 
the  initial  hydrogen  content  was  taken  from  the  equilibrium  diagram.  The  initial 
pressure  of  hydrogen  and  inerts  in  the  reservoir  tank  was  recorded  as  part  of 


mm.* a* 
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the  pressure-time  trace  on  the  recording  oscillograph.  The  reservoir  tempera¬ 
ture  was  obtained.  From  these  reservoir  conditons,  the  total  gas  leaving  the 
reservoir  as  a  function  of  time  could  be  calculated.  After  correction  for  the 
accumulation  of  gas  in  the  gas  space  of  the  retort,  the  total  hydrogen  and  inerts 
entering  the  porous  body  were  calculated. 


An  estimate  of  the  absorption  rates  to  be  encountered  during  these  tests  was 
obtained  using  the  HYDRA  code.  This  was  done  by  assuming  the  bodies  are  sub¬ 
jected  to  step  changes  in  hydrogen  pressure  at  the  surface.  The  time  constant 
can  be  defined  by 


AC  . 
AC  ~ 


1  -  e 


-t/r 


where 

AC  =  change  in  average  hydrogen  content  of  the  body 

AC^  =  ultimate  change,  after  an  infinite  time,  of  the  average  hydrogen 
content  of  the  body 

t  =  time 

T  =  time  constant 


Figure  1  shows  the  calculated  time  constant  for  a  body  one  centimeter  in 
depth  at  an  initial  temperature  of  1150°  F.  The  pressure  at  the  face  of  the  body 
is  stepped  from  0.5  to  1.5  atmospheres.  No  inert  gas  was  considered.  The 
bodies  were  assumed  to  be  composed  of  spherical  particles  having  average  diam¬ 
eters  of  0.  001  to  0.005  inch.  Void  fractions  varied  from  0. 15  to  0.25. 

Figure  2  shows  the  effect  of  temperature  on  the  predicted  time  constant  for 
the  same  range  of  particle  diameters  and  ihe  same  pressure  step. 

The  presence  of  an  inert  gas  tends  to  decrease  the  total  hydrogen  transport. 

It  may  also  speed  the  transient,  since  less  transport  will  result  from  the  pres¬ 
sure  step.  Figure  3  shows  the  relative  change  in  total  hydrogen  absorption  as 
a  function  of  the  percentage  inerts  for  a  typical  set  of  conditions.  The  curves 
were  generated  for  the  case  where  the  mole  fraction  or  partial  pressure  of  the 
inert  gas  did  not  change  during  the  transient,  which  was  the  situation  in  the 
present  experiments. 

Permeability  was  calculated  from  the  formula 


d 


,  Yfet 

AP 


11 


Temperature  *  1150*  F 
Pressure  step  «  0.  5  to  1.5  atm 

Bed  thickness  «  1.0  cm  Void  fraction  »  0.25 
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Figure  3.  Effect  of  Inert  Ges  on  Hydrogen  Absorption 


0% 


13 


AFWL-TR-67- 1 14 


where 

2 

d  *  permeability  (cm  ) 

V  *  superficial  gas  velocity  (cm/sec) 

M  ■  viscosity  of  hydrogen  (gm/cm-sec) 
t  •  thickness  of  body  (cm) 

2 

AP  »  pressure  difference  (gm/cm-sec  ) 

The  velocity  is  calculated  from  the  measured  mass  flow  rate  and  the  known 

density 


w  *  flow  rate  igm/sec) 

3 

p  ■  hydrogen  density  (gra /cm  ) 

9 

A  *  face  area  of  body  (cm  ) 

The  mass  flow  rate,  pressure  drop,  absolute  pressure,  temperature  and 
body  face  area  were  measured.  Literature  values  of  gas  density  and  viscosity 
were  used  in  the  calculation  of  permeability. 

Each  of  the  bodies  listed  in  table  II  was  subjected  to  a  minimum  of  three 

different  hydrogen  flow  rates  at  three  pressures  in  the  range  1  to  50  psia.  The 

cycled  bodies  were  permeability  tested  both  before  and  after  the  cycling.  The 

flow  rates  depend  upon  the  practical  limitations  of  the  pressure  drop  and  flow 

-4 

rate  measuring  devices,  i.e. ,  pressure  drop  0.02  to  20  psi  and  flow  10  to  10 
gal /min.  Figure  4  shows  the  predicted  pressure  drop  for  a  range  of  Reynolds 
numbers,  average  particle  diameter -shape  factor  products  and  pressure  levels. 
These  data  are  replotted  in  figure  5  to  show  the  corresponding  flow  rates  per 
unit  area  of  body  cross  section. 

Nominally,  the  bodies  were  cycled  between  lTOO0  and  1400°  F.  The  void 
volume  associated  with  the  retort  was  sized  so  that  the  body  will  vary  in  atom 
ratio  (H/Zr)  by  about  0. 10.  The  retort  pressure  will  cycle  between  0.  5  and  1 . 0 
atmosphere. 

c.  Tabulated  Plan  of  Runs 

Table  IV  gives  the  complete  schedule  of  absorption  runs  to  be  made.  Tiie 
letter  prefix  A  indicates  an  absorption  run  while  the  final  digit  indicates  dupli¬ 
cate  runs  at  these  condtions.  In  the  case  of  cycled  bodies,  a  final  8  indicates 
a  test  performed  prior  to  cycling.  A  final  9  indicates  a  post-cycling  test. 


14 


TOO  0 


AFWL-TH-67-114 


U_L 


rj 

O 


(uuj/isd)  mJu*T  jiuq  dojQ  jjnssjjj 


15 


Figure  4.  Anticipated  Pressure  Drops 
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Table  IV 

HYDROGEN  ABSORPTION  TEST  PARAMETERS 


Run 

number 


AO  10 
A020 
A030 
A040 
A050 
A060 
A070 
A080 
A090 
A100 
A110 
A120 
A130 
A140 
A150 
A160 
A170 
A180 
A190 
A200 
A210 
A220 
A230 
A240 
A250 
A260 
A270 
A280 
A290 
A3  00 
A3  10 
A320 


Body  ID 


SI 


Temperature 

(#F) 


1000 

1150 

1300 

1400 

1500 

1000 

1150 

1300 

1400 

1500 

1000 

1150 

1300 

1400 

1500 

1000 

1150 

1300 

1400 

1500 

1000 

1150 

1300 

1400 

1500 

1000 

1150 

1300 

1400 

1500 

1000 

1150 


0.  1  to  1. 1 
0.  1  to  1.  1 
0.  1  to  1.  1 
0.  1  to  1. 1 
0.1  to  1.  1 
0.5  to  1.5 
0.5  to  1.5 
0.5  to  1.5 
0.5  to  1.5 
0.5  to  1.5 
1.0  to  2.0 
1.0  to  2.0 
1.0  to  2.0 
1.0  to  2.0 
1.0  to  2.0 
0. 1  to  1. 1 
0.  1  to  1.  1 
0.1  to  1.  1 
0.1  to  1. 1 
0.  1  to  1 .  1 
0.5  to  1.5 
0.  5  to  1. 5 
0.5  to  1.5 
0.5  to  1.5 
Q.  5  to  1.5 
1.0  to  2.0 
1.0  to  ?.  0 
1.0  to  2.0 
1.0  to  2.0 
1.0  to  2.0 
0.1  to  1.  1 
0.1  to  1 .  1 


None 


None 

Argon 


Argon 

None 

None 


0.0 

0  05 


0. 05 
0.0 
0  0 


(table  cont’d  on  next  page* 
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Table  IV  (cont'  d) 


Run 

number 

Body  ID 

Temperature 

CF) 

Hydrogen 
pressure  step 
(atm) 

Diluent  gas 

Diluent 

pressure 

(atm) 

A330 

S2 

1300 

0. 1  to  1.1 

None 

0. 

0 

A340 

1400 

0. 1  to  1. 1 

A350 

1500 

0. 1  to  1. 1 

A360 

1000 

0.5  to  1.5 

A370 

1150 

0.5  to  1.5 

A380 

1300 

0. 5  to  1.5 

AS90 

1400 

0.5  to  1.5 

A400 

1500 

0.5  to  1.5 

A410 

1000 

1.0  to  2.0 

A420 

1150 

1.0  to  2.0 

A430 

1300 

1.0  to  2.0 

A440 

1400 

1.0  to  2.0 

A4S0 

1500 

1.0  to  2.0 

No 

ne 

0. 

0 

A460 

1000 

0. 1  to  1. 1 

Argon 

0. 

05 

A470 

1150 

0. 1  to  1.1 

A480 

1300 

0. 1  to  1. 1 

A490 

1400 

0. 1  to  1. 1 

A500 

1500 

0. 1  to  1. 1 

A510 

1000 

0. 5  to  1. 5 

A520 

1150 

0,5  to  1.5 

A530 

1300 

0.5  to  1.5 

A540 

1400 

0. 5  to  1 . 5 

A550 

1500 

0.5  to  1.5 

AS60 

1000 

1.0  to  2.0 

A570 

1150 

1.0  to  2.0 

A580 

1300 

1.0  to  2.0 

A590 

1400 

1.0  to  2.0 

A600 

S2 

1500 

1.0  to  2.0 

Ar 

gon 

u. 

05 

A610 

S3 

1000 

0.5  to  1.5 

None 

0. 

0 

A62C 

S3 

1150 

0.5  to  1.5 

None 

0. 

0 

A630 

S3 

1300 

0.  5  to  1.5 

None 

0. 

0 

(table  cont'd  on  nex  t  page) 


18 


AFWL-TR  -67- 1 14 


Table  IV  (cont '  d) 


Hydrogen 

Diluent 

Run 

Temperature 

pressure  step 

■  1 

pressure 

number 

Body  ID 

_ CF) 

latm) 

Diluent  gas 

latm) 

A640 

S3 

1400 

A650 

S3 

1500 

A660 

S4 

1000 

A670 

S4 

1150 

A680 

S4 

1300 

A690 

S4 

1400 

A700 

S4 

1500 

A710 

S5 

1000 

A720 

S5 

1150 

A730 

S5 

1300 

A740 

S5 

1400 

A750 

S5 

1500 

A760 

N 1 

1000 

A770 

N  1 

1150 

A780 

M 

1300 

A790 

N 1 

1400 

A800 

N I 

1500 

A8I0 

N2 

1000 

AS  20 

N2 

1150 

A830 

N2 

1300 

A840 

N2 

1400 

A8  50 

N  2 

1500 

A8S8 

S8 

1  150 

A869 

St> 

1  150 

A873 

s6 

1300 

A879 

S6 

noo 

A888 

S6 

1400 

A889 

S6 

1400 

A898 

S7 

1150 

AS  9  9 

S7 

1  15- 

A903 

57 

\+i 

o 

o 

0.  5  to  1 . 5  N  one  0.0 


i 


i 


i 


0.5  to  I  .  5  None  U.  0 
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Table  IV  (cont'  d) 


Run 

number 

Body  ID 

Temperature 

(°F) 

Hydrogen 

pressure  step 
(atm) 

Diluent  gas 

Diluent 

pressure 

(atm) 

A909 

S7 

1300 

0.5  to  1.5 

None 

0.0 

A918 

A919 

A928 

A929 

A938 

A939 

A948 

A949 

A958 

A959 

A968 

A969 

A978 

A979 

A988 

A9d9 

A998 

A999 

A1008 

A1009 

A1018 

A1019 

A1028 

A1029 


S7 
S7 
N1 
N1 
N1 
N 1 
N 1 
N 1 
FI 
FI 
FI 
FI 
FI 
FI 
F2 
F2 
F2 
F2 
F2 
F2 
F3 
F3 
F3 
F' 


1400 

1400 

1150 

1150 

1300 

1300 

1400 

1400 

1150 

1150 

1300 

1300 

1400 

1400 

1150 

1150 

1300 

1300 

1400 

1400 

1150 

1150 

1300 

1300 


A1038 

i 

1400 

I 

1 

A1039 

F 

3 

1400 

None 

0.0 

A1040 

Selected 

1150 

Krypton 

0.  15 

A 1050 

1150 

Krypton 

0.015 

A1060 

1150 

Krypton 

Selected 

A1070 

1150 

Xenon 

0. 15 

A1080 

1150 

Xenon 

0.015 

A1090 

Selected 

1150 

0.5  to  1.5 

Xenon 

Selected 

(table  cont’d  on  next  page) 
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Table  IV  (cont'  d) 


Hydrogen 

Diluent 

Run 

Temperature 

pressure  step 

pressure 

number 

Body  ID 

(°F) 

(atm) 

Diluent  gas 

(atm) 

A1100  Selected 

A1110 

A1120 

A1130 

A1140 

All  50  Selected 


1150 

1300 

1400 

1150 

1300 

14C0 


0.5  to  1.5  Oxygen  Selected 


0.5  to  1.5  Oxygen  Selected 


Table  V  gives  the  nominal  operating  conditions  of  the  cycling  tests.  The 
prefix  C  indicates  a  cycling  run. 


Table  V 

THERMAL  CYCLING  TEST  PARAMETERS 


Run 

Temperature  cycle 

Pressure  cycle 

number 

Body  ID 

(°F) 

(atm) 

C010  S6  1200  to  1400  0.5  to  1.0 

C020  S7 

C030  FI 

C040  F2 

C050  F3 

C060  N3  1200  to  1400  0.5  to  1.0 

5.  Apparatus  and  Test  Procedures 

a.  Permeability 

(1)  Apparatus 

A  schematic  of  the  main  components  of  the  permeability  system  is  shown  in 
figure  6,  including  valves,  flowmeters,  specimen  holder  assembly  and  pressure 
transducer  connections.  Provisions  were  included  for  flow  and  pressure  measure¬ 
ment.  Calibration  curves  prepared  for  each  flowmeter  were  used  to  determine 
hydrogen  flow  rates  for  a  given  pressure  and  meter  reading.  Upstream 
pressures  were  measured  by  using  the  pressure  gage  and  the  absolute  pres¬ 
sure  transducer.  Pressure  drops  across  the  test  specimen  were  measured 
by  using  the  differential  pressure  transducers.  Output  of  these  three  trans¬ 
ducers  were  recorded  on  an  oscillograph. 
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To  vacuum 


system 


Figure  6.  Hydrogen  Permeation  Test  Apparatus 
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(2)  Procedure 

The  procedure  for  performing  the  hydrogen  permeability  tests  was  as 
follows: 

(1)  Load  the  test  specimen  into  the  holder  assembly,  and  evacuate  the-  sys- 

-3 

tern  to  10  torr  by  opening  valves  through  P^,  excluding  Pg.  In  addition, 
valves  Ag,  Ag  and  A^  were  open,  while  valves  A2,  Ag,  and  Ag  were  closed  (Cf. 
figures  6  and  7). 

(2)  Isolate  vacuum  pump  at  P^  and  check  for  vacuum  leaks.  Close 
valves  opened  in  (1). 

(3)  Depending  on  the  flow  rate  desired,  open  the  valves  at  entrance  and 
exit  sections  of  the  selected  flowmeter. 

(4)  Open  valve  Pj^  and  valve  P^  or  P^  to  initiate  flow.  P^  was  used  to 
obtain  low  pressure  data  while  P2  was  used  for  certain  high  pressure  runs. 

(5)  Make  the  valve  adjustments  required  to  obtain  desired  mass  flow  rate 
and  pressure  level. 

All  permeability  test  data  pertaining  to  each  specimen  were  documented,  in¬ 
cluding:  identification  number,  flow  rates  and  corresponding  pressure  levels  and 
pressure  drops.  In  addition,  it  was  noted  whether  the  test  specimen  had  pre¬ 
viously  undergone  thermal  cycling. 

b.  Hydrogen  Absorption 

(1)  Apparatus 

Basic  components  of  the  hydrogen  absorption  system,  including  arrangement 
of  vacuum-pressure  lines,  are  shown  by  figure  7.  Valves  are  numbered  to  assist 
procedure  instructions.  Two  gas  supply  lines  were  provided,  one  for  hydrogen 
(valve  AjQ)  and  one  for  inert  gas  (valve  Ag). 

This  system  was  so  designed  that  with  an  initial  pressure  differential  of  14.7 
psi  between  the  retort  and  the  surge  tank,  the  maximum  pressure  decay  within 
the  surge  tank,  after  opening  the  solenoid  valve  (A^),  will  be  2  psi.  Pressure 
instrumentation  includes  two  gages  (0  to  100  psi),  one  on  either  side  of  valve  A^ 
and  three  transducers.  The  absolute  pressure  transducer  (0  to  20  psi)  has  a 
dynamic  vacuum  on  the  negative  side  with  the  retort  pressure  on  the  positive 
side.  The  differential  pressure  transducers  can  be  operated  with  either  a 
dynamic  vacuum  or  surge  tank  pressure  on  the  negative  side.  Surge  tank 
pressure  was  monitored  on  the  positive  side  of  the  transducer. 
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Figure  7.  Hydrogen  Absorption  Test  Apparatus 
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Furnace  control  instrumentation  is  shown  in  figure  8.  The  Simplytrol  meter 
was  set  to  shut  off  the  furnace  power  whenever  the  furnace  control  thermocouple 
indicated  1550°  F  or  above.  Retort  thermocouples  are  connected  to  the  millivolt 
recorder.  The  sheath  (specimen)  thermocouple  is  connected  to  the  oscillograph, 
as  are  the  pressure  transducers. 

(2)  Procedure 

The  procedure  for  performing  hydrogen  absorption  tests  is  as  follows: 

(1)  Lo^J  the  specimen  into  a  retort  as  indicated  by  figure  9,  making  cer¬ 
tain  that  the  sheathed  thermocouple  is  in  contact  with  the  sample. 

.Q 

(2)  Connect  retort  to  system.  Evacuate  entire  system  to  10  torr  by 
opening  valves  Aj  through  Ag  (figure  7).  Valve  Cg  of  figure  10,  and  valves  P4 
and  Pg  of  figure  6  remain  closed,  Check-for  vacuum  leaks  by  isolating  vacuum 
pump  (valve  A^). 

(3)  Close  valves  Aj,  Ag,  A4,  Ag,  A^andA^. 

(4)  Heat  specimen  to  test  temperature. 

(5)  After  reaching  the  pressure -temperature  equilibrium  of  the  specimen, 
pressurize  surge  tank  to  a  pressure  14.  7  psi  greater  than  that  of  the  report.  If 
desired,  a  mixture  of  hydrogen  and  inert  can  be  obtained  by  manipulation  of 
valves  Ag  and  A1(). 

(6)  Close  valves  Ag,  Ag  and  Alg. 

(7)  Start  oscillograph,  recording  specimen  temperature  and  the  output  of 
the  pressure  transducers. 

(8)  Open  the  solenoid  valve  A4>  If  the  subsequent  pressure  drop  in  the 
surge  tank  exceeds  0.2  psi,  close  valve  A^. 

(9)  After  pressure  equilibrium  of  surge  tank-retort,  close  valve  A4> 

(10)  Proceed  to  next  test  condition. 

(11)  After  test  completion,  isolate  and  cool  the  retort  to  room  temperature 
and  remove  the  test  specimen. 

All  absorption  test  data  pertaining  to  each  specimen  were  recorded  as 
follows:  identification  number,  weight  before  and  after  test,  dimensions,  com¬ 
position,  and  surface  condition.  Hydrogen  content  was  determined  by  equilibrium 
pressure  temperature  data  and  specimen  weight  change.  Because  the  retort  was 
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Figure  8.  Furnace  Control  Console 
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not  at  a  uniform  temperature,  the  amount  of  hydrogen  in  this  void  volume  was 
difficult  to  calculate  accurately.  However,  an  "equivalent"  void  volume  was 
measured.  With  the  use  of  a  dummy  test  specimen,  the  retort  was  pressurized 
at  room  temperature,  the  specimen  temperature  was  increased  to  T^,  and  the 
corresponding  increase  in  pressure  was  measured.  The  equivalent  volume  was 
then  calculated  by  the  equation: 

P  T. 
e  c  Phc 

where 

*  equivalent  volume  occupied  by  the  gas  at  temperature  Th 

T  *  room  temperature 
c 

Pc  *  pressure  at  Tc 

V  «  gas  volume  at  T 
c  c 

*  specimen  temperature  after  heating 

Ph  *  pressure  at 

In  this  manner,  an  equivalent  volume -temperature  correlation  was  determined 
to  cover  the  temperature  range  of  hydrogen  absorption  testing.  Thus,  the  moles 
of  hydrogen  which  occupy  this  space  during  testing  were  calculated.  Similar 
volume  determination  tests  were  performed  on  the  thermal  cycle  retorts. 

Initial  pressures  and  temperatures  of  the  surge  tank  and  retort  were  recorded, 
as  well  as  tae  oscillograph  traces  showing  subsequent  pressure -temperature 
changes  after  valve  A4  was  opened. 

c.  Thermal  Cycling 

(1)  Apparatus 

The  schematic  of  figure  10  shows  the  basic  components  of  the  thermal  cycle 
system,  including  the  arrangement  of  the  vacuum -pressure  lines.  Each  valve 
was  numbered  to  simplify  procedural  instructions.  Furnace  control  instrumenta¬ 
tion  is  shown  by  figure  8.  The  Simplytrol  meter  is  set  to  shut  off  the  furnace 
power  whenever  the  furnace  control  thermocouple  or  indicator  rises  to  1450°  F. 

The  timer  switch  is  set  to  cycle  the  furnace  between  1200®  and  1400°  F.  The 
sheath  (specimen)  thermocouple  and  the  retort  thermocouples  are  connected  to 
the  millivolt  recorder.  The  pressure  transducers  are  connected  to  the  oscillograph. 
Transducers  operate  with  a  dynamic  vacuum  on  the  negative  side,  and  with  the  surge 
tank  pressure  on  the  positive  side.  A  pressure  line  from  each  retort  is  connected 
to  the  pressure  recorder  (0  to  100  psi). 
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This  system  was  designed  to  cycle  the  ZrHx  (or  ZrUHx)  specimens  between 
1200°  and  1400°  F,  with  a  corresponding  pressure  cycle  of  7.4  to  14.8  psi.  This 
represents  about  0.02  gram  of  hydrogen  movement,  or  six  percent  of  the  initial 
hydrogen  in  the  specimens.  Provision  was  included  for  the  continuous  monitor¬ 
ing  of  temperature  and  pressure  throughout  testing, 

(2)  Procedure 

The  procedure  for  performing  thermal  cycling  tests  is  as  follows: 

(1)  Load  test  specimens  into  retort  as  indicated  by  figure  9  making  certain 
that  the  sheathed  thermocouple  is  in  contact  with  the  sample. 

-3 

(2)  Connect  retorts  to  system.  Evacuate  entire  system  to  10  torr  by 
opening  valves  and  Cg.  Valves  Aj  and  Aj^  of  figure  7  are  closed.  Check 
for  vacuum  leaks  by  isolating  vacuum  pump. 

(3)  Close  valves  C^,  C4  and  Cg.  Heat  specimens  to  1200°  F. 

(4)  After  temperature  equilibrium  is  attained,  adjust  pressure  of  retorts 
to  7. 4  psi. 

(5)  With  valves  C,  and  opened  and  the  remaining  valves  closed,  initiate 
thermal  cycling  by  adjustment  of  the  timer  switch  (Cf.  figure  8). 

(6)  After  test  completion  (200  hours  each  of  continuous  cycling),  the 
specimens  are  cooled  to  room  temperature  and  removed  from  the  retorts. 

All  cycling  test  data  pertaining  to  each  specimen  are  recorded  as  follows: 
identification  number,  weight  before  and  after  test,  dimension,  composition 
and  surface  condition.  Amounts  of  hydrogen  added  to,  or  taken  from,  the  test 
specimen  are  noted.  In  addition,  typical  strip  chart  recordings  of  temperature 
and  pressure  cycles  are  preserved. 

6.  Experimental  Performance  and  Method  of  Data  Reduction 
a.  Permeability 
(1)  Experiment 

A  thorough  discussion  of  the  experimental  apparatus  and  test  procedure 
was  given  in  the  preceding  paragraphs.  In  brief,  a  gas  underpressure  was 
admitted  at  the  end  of  the  specimen  and  allowed  to  flow  through  the  sintered 
body.  The  observed  data  consisted  of  t He  gas  pressure  and  temperature,  gas 
mass  flow  rate,  the  pressure  drop  over  the  length  of  the  sintered  body  and  the 
sintered  body  physical  properties. 
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(2)  Data  reduction 

Reduction  of  the  raw  data  was  performed  by  the  digital  code  PERMA.  The 
permeability  coefficient  was  evaluated  from  Darcy' s  Law  which  states  that  the 
velocity  of  a  gas  (V)  through  a  bed  of  solids  is  directly  proportional  to  the  pres 
eure  drop  (dP/dx),  and  inversely  related  to  the  viscosity  of  the  gas  (p),  where 
the  constant  of  proportionality  is  the  permeability  (d). 

Because  the  mass  flow  rate  (w)  is  related  to  the  velocity  by 

w • pVA, 

and  employing  the  ideal  gas  law  to  evaluate  the  density, 

equation  (13)  becomes 

“  - R  0r)  (^)  (&)  «>«> 

where 

R  *  ideal  gas  constant 

AX  *  specimen  thickness 

A  *  total  cross  section  area  of  specimen  perpendicular  to  the  gas  flow 
p  *  viscosity  of  the  gas  at  temperature  (T)  and  pressure  (P) 

P  ■  pressure  drop  across  the  specimen  for  a  gas  mass  flow  rate  of  w. 

Equation  (13)  defines  the  pressure  drop  equal  to  zero  for  a  gas  velocity 
equal  to  aero.  Thus,  if  one  were  to  plot  the  gas  velocity  versus  the  pressure 
drop  on  a  Cartesiar  graph,  a  straight  line  through  the  origin  would  result  with 
a  slope  proportional  to  the  permeability  constant.  It  was  observed  early  in  the 
data  reduction  that  measurements  at  low  mass  flow  rates  produced  relatively 
erratic  results.  The  cause  was  eventually  traced  to  instrument  error  bands. 
However,  to  evaluate  t?  e  data,  each  data  point  was  assigned  a  weighing  factor 
fW)  based  upon  the  "lever"  or  "moment"  principle. 
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Thus,  the  farther  the  data  point  was  from  the  origin,  the  higher  the  weighing 
factor  (W  =  1).  The  average  permeability  was  then  evaluated  by  the  follow¬ 
ing: 


L 

i 


W  d. 
1  i 


£  W. 


This  technique  of  data  reduction  provided  a.n  excellert  means  of  screening  the  ex¬ 
perimental  data  and  an  immediate  identification  of  questionable  data  points 

The  standard  deviation  (SD)  of  the  experimental  permeability  from  the  average 
value  was  calculated  by 


i=  1 


The  range  of  the  various  parameters  is  given  in  table  VI.  Permeability 
measurements  were  made  before  and  after  the  hydrogen  absorption  tests  and 
before  and  after  thermal  cycling  of  the  specimens. 

Table  VI 


RANGE  OF  PARAMETERS --PERMEABILITY  TEST 


Parameter 

Range 

Pressure  (psia) 

•  35 

Pressure  drop  (psia) 

9,  0  \  to  0. j 

Gas  mass  flow  (gm/run) 

0. 005  to  0.  5 

Specimen  void  fraction  (<r«) 

20  to  35 

Specimen  average  particle  size  (mils) 

1.0  to  6.5 

Type  of  gas 

A.  H2 

b.  Hydrogen  Absorption  Test 
(1)  Experiment 

The  experimental  apparatus  was  divided  into  two  sections  the  retort  section 
which  contained  the  test  specimen  in  an  Inconel  tube  and  furnace  assembly  and 
the  surge  tank  section  which  contained  the  pressurized  hydrogen  gas  at  room 
temperature,  i  hese  two  sections  were  connected  by  sta  r  less  steel  tubing 
through  a  solenoid  valve.  The  temperature  of  the  top  surface  of  the  specimen 
was  measured  by  a  thermocouple  which  protruded  thiou  gh  the  retort  flange  and 
extended  down  to  the  specimen  area. 
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The  hydrogen  pressure  and  specimen  temperature  in  the  retort  were  allowed 
to  come  to  a  steady  state  equilibrium.  The  surge  tank  was  -hen  pressurized  to  a 
level  of  approximately  one  atmosphere  higher  than  that  of  the  retort.  Because 
the  system  volumes  were  predetermined,  the  mass  of  hydrogen  in  each  section 
was  known.  The  solenoid  valve  separating  the  two  sections  was  then  opened  for 
greater  than  0. 1  second  but  less  than  one  second.  The  pressure  drop  in  the 
surge  tank  and  the  pressure  rise  and  subsequent  decay  in  the  retort  were  meas¬ 
ured  by  preesure  transducers  and  recorded  on  an  oscillograph. 

The  obs<  rved  data  consisted  of  the  temperature  and  the  initial  and  final  pres¬ 
sure  in  the  surge  tank,  and  the  initial  temperature  and  pressure  and  pressure 
variation  v'ith  time  in  the  retort. 


Initial  hydrogen  pressure  and  initial  temperature  were  varied  to  observe 
their  effects  on  the  absorption  process  for  a  given  porous  body.  Specimens  of 
varying  particle  size  and  void  fraction  were  tested  in  this  way.  To  study  the 
influence  of  inert  gases,  the  surge  tank  was  initially  pressurized  with  a  mixture 
of  hydrogen  ard  inert  gas.  Gases  investigated  wero  argon,  krypton  and  xenon. 

(2)  Data  reduction 


Tne  reduction  of  the  raw  data  was  accomplished  using  the  digital  computer 
code  RESPND  to  solve  the  hydrogen  continuity  equation.  At  any  time  (t  >0): 


Change  in  H2  mass 
(surge  tank) 


change  in  H2  mass  +  mass  of  H2  absorbed  +  leakage 
v*ctort)  (specimen) 


1.008  ' 


Zr  AX 


91  22 


+  Lt  (15) 


where 

AP  =  time  dependent  pressure  differential  (Pt  -  Pp) 

V  -  volume 

T  =  temperature 

R  =  gas  constant 

w =  weight  cf  zirconium  specimen 

AX  -  change  in  H/Zr  atom  ratio  (Xt  -  X^) 

L  =  leakage  rate 

t  =  accumulated  test  time 
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Equation  (15)  was  used  to  reduce  the  expei  ^mental  data  to  a  time  history  of  the 
percent  of  change  in  hydrogen  content  of  the  specimen  because  of  absorption: 

%  H2  change  =  ^  x  100%  06) 

o 

In  selected  cases,  results  obtained  fiom  the  data  reduction  code  were  com¬ 
pared  with  those  obtained  from  the  digital  computer  code,  HYDRA  II. 

c.  Thermal  Cycling  Test 

(1)  Experiment 

Four  nonuranium -containing  hydride  specimens  were  subjected  to  the  thermal 
cycling  test.  The  procedure  consisted  of  heating  the  specimens  to  approximately 
1000°  F  under  a  dynamic  vacuum.  The  specimens  were  then  rehydrided  to  vary¬ 
ing  initial  hydrogen-to-zirconium  atom  ratios  (H/Zr),  by  proper  selection  of 
initial  hydrogen  pressure  and  specimen  temperature.  The  minimum  and  maxi¬ 
mum  temperature  points  were  then  set  on  a  temperature  controller  and  the 
cycling  commenced. 

(2)  Data  reduction 

The  raw  data  generated  by  the  cycling  tes:  consisted  of  temperature  and 
pressure  time  histories  recorded  on  chart  paper.  These  records  were  then 
scanned  to  obtain  the  minimum  and  maximum  temperatures  and  pressure  varia¬ 
tion  in  order  to  evaluate  the  net  hydrogen  movement  during  a  cycle. 

7.  Experimental  Program 

a.  Results  and  Analysis 

Thirteen  sintered  specimens  of  varying  powder  particle  sizes  and  shapes 
and  void  fractions  were  fabricated  for  this  study.  The  pertinent  body  character¬ 
istics  of  the  spherical  and  nonspherical  zirconium  hydride  bodies  are  listed  in 
table  VII.  The  fuel  bearing  specimens  are  described  in  the  classified  volume  II. 

The  particle  si^es  used  to  prepare  the  sintered  compacts  are  reported  as  a 
function  ot  bracketing  USBS  sieve  numbers.  For  example,  the  particle  diameter 
of  a  powder  sample  that  went  through  an  80 -mesh  sieve  but  did  not  go  through  a 
100- mesh  sieve  (-80  +100)  would  be  in  the  range  of  7.0  to  5.  9  mils.  This  means 
of  determining  the  powder  particle  size  introduced  an  uncertainty  cf  at  most 

8.  5  percent  based  on  the  mean  value.  This  uncertainty  is  further  complicated 
when  two  mesh  fractions  arc  combined.  Although  this  error  in  measured 
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particle  size  did  not  affect  the  data  reduction  of  the  experimentally  determined 
permeability,  it  did  contribute  to  uncertainties  in  the  theoretical  value  for 
permeability  as  calculated  from  the  Ergun  correlation.  Also,  use  of  this 
average  particle  size  in  the  hydrogen  absorption  predictions  (HYDRA)  introduced 
uncertainties  into  the  analysis. 

Table  VII 


SINTERED  BODY  CHARACTERISTICS 


Specimen* 

Powder  mesh  size 

Average 

particle 

diameter 

(mils) 

Void 

fraction 

SI 

100%  of  -270  +400 

1.8 

0.2229 

S2 

100%  of  -270  +400 

1.8 

0.2346 

S3 

100%  of  -150  +170 

3.8 

0.2336 

S4 

100%  of  -150  +170 

3.8 

0.2042 

S5 

100%  of  -400 

1.0 

0.2095 

S6 

100%  of  -80  +100 

6.5 

0.22682 

S7 

100%  of  -80  +10^ 

6.5 

0.2342 

N 1 

80%  of  -80  +100 

20%  of  -325  +400 

5.52 

C. 2492 

N2 

80%  of  -80  +100 

20%  of  -100  +140 

6.2 

0.2733 

N3  80%  of  -140  +200 

20%  of  -325  +400 
*S  3  spherical  powder  specimen 

N  3  nonspherical  powder  specimen 

3.1 

0.2461 

(a)  Permeability 

The  empirical  correlation  for  the  permeability  of  a  packed  bed  as  a  function 
of  the  bed  characteristics  proposed  by  Ergun  states  that  in  the  case  of  laminar 
flow  (Reynolds  number  <  10)  the  permeability  is  predicted  by 


D2  f2 

d  =  Vf 


150  (l-«)‘ 


(17) 
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where 


d 

D 

P 


f 


F 


=  permeability 
=  particle  diameter 
=  void  fraction 
=  shape  factor 


A  comparison  of  the  empirical  or  Ergun  permeability  evaluated  from  equa¬ 
tion  (17)  and  the  experimental  value  evaluated  from  equation  (14)  was  given  in 
table  VIII  for  pre-  and  post-hydrogen  absorption  or  thermal  cycle  test  conditions. 
The  physical  properties  of  the  sintered  bodies  (grain  or  particle  size,  void  frac¬ 
tion,  etc.)  are  a  function  of  the  sintering  temperature  and  the  length  of  time  a 
body  is  subjected  to  this  temperature.  Thus,  as  long  as  the  test  conditions  are 
lower  than  the  minimum  sintering  temperature  (-1700°  F),  one  should  observe 
negligible  change  in  the  body  structure  and  the  effects  of  continual  absorption 
and  desorption  of  hydrogen  and  phase  changes  should  be  self-evident. 


Table  VRI 

HYDROGEN  PERMEABILITY 


Permeability 

x  10+9  (cm2) 

Specimen 

Ergun 

Pretest  ! 

Post- test 

SI 

2.556 

4.981 

S2 

3.069 

4.896 

4.780 

6.725 

S3 

13.481 

16.469 

20.529 

20.259 

S4 

8.355 

11.357 

11.941 

S5 

0.633 

1.  157 
1.313 

1.423 

S6 

35.471 

— 

35.720 

S7 

39.804 

-- 

54.873 

N 1 

35.992 

-- 

7.626 

N2 

63.941 

19.753 

21.097 

N3 

10.846 

4.528 

4.364 

Examination  of  the  pre-  and  post-test  permeabilities  showed  that  essentially 
no  change  occurred.  Variations  are  within  the  experimental  error. 
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In  order  to  compare  the  permeability  of  the  porous  sintered  bodies  with  that 
predicted  by  Ergun  for  a  packed  bed  containing  the  same  void  fraction  and  particle 
size,  an  effective  particle  size  (fpXDp)  evaluated  by  a  rearrangement  of  equation 
10  is  defined. 


Effective  particle  size  (fp 


««v-(lsaV-,d) 


1/2 


A  plot  of  the  effective  particle  size  based  on  the  experimentally  determined 
permeability  versus  the  actual  or  measured  particle  size  as  determined  from 
the  sieve  sizes  is  shown  in  figure  1 1 .  The  slope  of  these  lines  is  particle  shape 
factor.  Theoretically,  the  shape  factor  for  spherical  particles  is  unity,  while 
nonspherical  particles  have  shape  factors  of  less  than  one.  As  noted  on  figure 
11,  the  shape  factor  for  bodies  composed  of  sintered  spherical  particles  is 
slightly  greater  than  one.  This  higher  value  seems  real  and  cannot  be  attributed 
to  any  experimental  error.  It  can  be  explained  by  the  fact  that  the  void  structure 
in  a  sintered  body  is  more  open,  that  is,  the  voids  form  channels,  thus,  lowering 
the  pressure  drop  for  a  given  flow  rate.  The  sintered  bodies  composed  of  non¬ 
spherical  particles  have  shape  factors  less  than  unity.  The  exact  value  of  the 
shape  factor  was  obscured  because  of  the  "ncertainty  in  determining  the  actual 
partical  size  of  the  powder.  However,  the  results  are  in  close  agreement  with 
expected  behavior. 


Prior  to  starting  the  test  program,  it  was  postulated  that  there  might  be  a 
pressure  effect  upon  the  permeability.  As  shown  in  figure  12,  the  hydrogen 
permeability  is  virtually  constant  over  the  range  of  pressures  investigated.  It 
should  be  noted  that  each  data  point  on  the  graph  represents  numerous  runs  of 
varying  hydrogen  gas  mass  flow  rates  and  pressure  drops  at  constant  pressure. 
Variations  can  be  attributed  to  inherent  experimental  errors. 

The  effect  of  the  constituent  gas  on  permeability  was  also  investigated  by 
using  hydrogen  and  argon  for  two  pretest  bodies.  The  results  of  this  experiment 
are  gi/en  in  table  IX.  The  permeability  of  the  porous  sintered  bodies  was  ap¬ 
parently  independent  of  the  type  gas  used. 

b.  Hydrogen  Absorption 

The  primary  objective  of  this  study  was  the  evaluation  of  the  hydrogen  absorp¬ 
tion  response  of  porous  zirconium  hydride  bodies  for  varying  conditions.  In  per¬ 
forming  this  study,  it  was  necessary  to  examine  the  effect  of  numerous  combina¬ 
tions  and  permutations  f  the  pertinent  parameters.  Rather  than  inundate  the 
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Figure  11.  Hydrogen  Permeability  for  Porous  ZrHx  Sintered  Bodies 
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reader  with  a  myriad  of  data  plots  which  may  obscure  the  facts,  the  findings  of 
this  study  will  be  discussed  and  specific  plots  will  be  introduced  to  substantiate 
the  experimental  observations. 


Table  IX 

EFFECT  OF  CONSTITUENT  GAS  ON  PERMEABILITY 


Specimen 

Permeability  x 

109  (cm2) 

Hydrogen  | 

Argon 

S2 

4.896 

4.346 

S3 

16.469 

16.940 

S3 

20.529 

16.370 

(1)  Temperature  effect 

Increasing  the  initial  steady  state  temperature  of  a  specimen  results  in  the 
following: 

(1)  An  increase  in  the  hydrogen  mass  diffusion  coefficient 

(2)  An  increase  in  the  H/Zr  atom  ratio  differential  for  a  constant  pres¬ 
sure  change. 

Thus,  the  temperature  effect  would  be  one  of  increasing  the  rate  of  absorption 
with  increasing  specimen  temperature.  Table  X  illustrates  these  effects  by  com¬ 
paring  the  percent  change  in  the  H/Zr  atom  ratio  after  a  test  time  of  30  seconds 
as  a  function  of  initial  specimen  temperature.  The  pressure  function  in  each 
case  was  a  ramp  change  in  hydrogen  pressure  from  0.  5  to  1. 5  atmospheres. 

Thus,  an  increase  in  the  rate  of  hydrogen  absorption  with  increasing  initial  tem¬ 
peratures  was  observed.  This  relationship  is  not  a  simple  one  because,  as 
shown  in  figure  13,  two  different  mechanisms  are  prevalent.  In  the  first  region 
(lower  temperature  region)  the  variation  of  the  percent  of  hydrogen  absorbed  is 
almost  linear  with  temperature.  The  second  region  exhibits  an  extreme  tem¬ 
perature  sensitivity.  The  particle  diameter  determines  the  slope  of  the  linear 
portion  of  the  curve.  For  small  diameters,  mass  diffusion  of  hydrogen  proceeds 
rapidly  and  is  less  sensitive  to  temperature.  However,  as  the  particle  diameter 
increases,  the  temperature  and,  consequently,  the  diffusion  coefficient  have  an 
increasing  effect  (steeper  slope)  upon  the  hydrogen  absorption  response.  Figure 
14  illustrates  the  same  phenomena  for  the  nonspherical  particle  specimens.  This 
figure  is  a  clear  example  of  the  effect  of  particle  size  as  predicted  earlier.  In 
the  low  temperature  range,  the  small  particle  size  showed  more  rapid  response 
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Hydrogen  Absorbed  After  60  Seconds  {%) 
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Figure  13.  Hydrogen  Absorption  Response  as  Function  of  Initial  Specimen  Temperatur 
and  Particle  Size 


Hydrogen  Absorbed  After  60  Seconds  (%) 
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Figure  14.  Hydrogen  Absorption  Response  as  Function  of  Initial  Specimen  Temperature 
and  Size 
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because  the  rate  of  absorption  was  limited  by  the  diffusion  process.  At  higher 
temperatures,  the  order  was  reversed  because  diffusion  was  no  longer  an  im¬ 
portant  rate-determining  mechanism  and  body  permeability  becomes  dominant. 
Figures  15  and  16  show  the  percent  change  in  hydrogen  concentration  for  the 
spherical  and  nonspherical  particle  specimens  as  a  function  of  test  time. 

Table  X 


HYDROGEN  ABSORPTION  RESPONSE  FOR  HYDROGEN  PRESSURE 
RAMP  OF  0.  5  TO  1.5  ATMOSPHERES  AS  FUNCTION  OF 
INITIAL  SPECIMEN  TEMPERATURE 


scimen 

Particle 

diameter 

(mils) 

Change  in  H/Zr  Atom  Ratio  After  60  Seconds  (%) 

1000*  F 

1  1150°  F 

1  1300°  F 

I  1400°  F 

|  1400°  F 

SI 

1.8 

0.919 

1. 520 

1.875 

2.  135 

2.657 

S2* 

1.8 

ND 

ND 

2.083 

2.715 

3. 540 

S3* 

3.8 

0.018 

0.  042 

ND 

1.652 

3.086 

S4 

3.8 

0.633 

1.273 

1.824 

2.  115 

2.  573 

2.516 

S5 

1.0 

1.206 

1.877 

1.947 

2.281 

2.739 

S6 

6.5 

ND 

0.541 

1.368 

ND 

UD 

S7 

6.  5 

ND 

0.337 

1.397 

1.897 

ND 

N 1 

5.52 

0.419 

ND 

1.810 

1.9S9 

2.633 

N2 

S.  2 

ND 

1.425 

1.727 

2.003 

2.850 

N3 

3.  1 

ND 

1. 601 

1.888 

2.  942 

ND 

^Extensive  surface  contamination  of  specimen  in  the  form  of  an  oxide  film. 

ND  *  no  data 

(2)  Pressure  effect 

It  was  previously  mentioned  that  increas.ng  the  initial  specimen  temperature 
resulted  in  an  increase  in  the  initial  H/Zr  atom  ratio  differential  for  a  constant 
pressure  change.  Decreasing  the  initial  pressure  had  the  same  effect.  Figure 
17  illustrates  the  initial  surface  concentration  (H /Zr  atom  ratio)  differential  for 
a  one -atmosphere  hydrogen  pressure  change  as  a  function  of  the  initial  system 
pressure  and  specimen  temperature.  Thus,  higher  rates  of  absorption  for 
lower  initial  system  pressures  could  be  expected.  Table  XI  shows  the  experi¬ 
mental  hydrogen  absoro'ion  response  for  a  hydrogen  pressure  change  of  one 
atmosphere  as  a  function  of  initial  hydrogen  pressure  and  initial  specimen 
temperature. 
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Figure  16.  Hydrogen  Absorption  Response  for  Nonspherical  ZrH  as  Function  of  Initial 
Specimen  Temperature 
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HYDROGEN  ABSORPTION  RESPONSE  FOR  HYDROGEN  PRESSURE 
RAMP  OF  ONE  ATMOSPHERE  AS  FUNCTION  OF  INITIAL  HYDROGEN 
PRESSURE  AND  INITIAL  SPECIMEN  TEMPERATURE 


Specimen 

Pressure  ramp 
(atm) 

Change 

in  H/Zr  Atom  Ratio  After  60  Seconds  (%) 

1000°  F 

IBE3QI 

I  1300°  F 

1400°  F 

1500°  F 

SI 

0.1  to  1.1 

1.931 

3.004 

2.  980 

3.706 

4.783* 

0.5  to  1.5 

0.919 

1.520 

1. 875 

2.  135 

2.657 

1.0  to  2.  0 

ND 

0.938 

1.  173 

1.380 

1.642 

S2 

0.1  to  1.1 

ND 

ND 

3.  126 

4.579 

5.037* 

0.5  to  1.5 

ND 

ND 

2.083 

2.715 

3.  540 

1.0  to  2.  0 

ND 

ND 

1.431 

1.730 

1. 765 

*ZrH  specimen  was  initially  in  |3  phase  and  transformed  to  the  6  phase  by  the 
end  of  the  transient. 

ND  =  no  data 


As  expected,  the  transient  with  the  highest  initial  specimen  temperature  and  the 
lowest  initial  hydrogen  pressure  produced  the  largest  change  in  nydrogen  absorp¬ 
tion.  Figures  18  and  19  depict  the  percent  change  in  hydrogen  concentration 
versus  time  as  a  function  of  the  initial  specimen  temperature  and  hydrogen  pres¬ 
sure  for  a  constant  hydrogen  pressure  ramp  change  (one  atmosphere). 

(3)  Particle  size  effect 

Referring  to  table  X,  it  is  noted  that  decreasing  the  particle  diameter  of  the 
Zill  specimen  increases  the  absorption  of  hydrogen.  However,  as  discussed  in 
the  section  on  temperature  effects,  the  influence  of  the  particle  size  effect  dimin¬ 
ished  with  increasing  temperature  (see  figures  13  and  14).  This  results  from  the 
fact  that  the  diffusion  coefficient  becomes  sufficiently  high  that  the  rate  at  which 
hydrogen  is  absorbed  depends  on  other  processes. 

The  sphere  has  the  lowest  surface-to-volume  ratio  of  all  geometric  figures. 
Thus,  for  a  given  solids  volume,  the  nonspherical  particle  specimens  should 
have  a  larger  surface  area.  Because  the  hydrogen  absorption  phenomena  is  a 
strong  function  of  the  surface  area,  one  would  expect,  body  permeabilities 
being  equal,  that  the  nonspherical  particle  specimens  would  exhibit  better  hydro¬ 
gen  absorption  characteristics  when  compared  to  the  spherical  particle  specimens. 
Examination  of  the  experimental  data  presented  in  table  X  shows  that  for  approxi¬ 
mately  the  same  specimen  properties  and  test  conditions,  the  percent  change  in 


Hydrogen  Absorbed  (%) 


Pressure  ramp  change  =  1  atmosphere 


Initial  Hydrogen 
Pressure 
_ (atm) _ 


F’igure  18.  Hydrogen  Absorption  Response  for  Spherical  ZrH  as  Function  of  Initial 
Specimen  Temperature  and  Hydrogen  Pressure 


hydrogen  concent  rat  inn  is  higher  )o;  the  nonspherical  particles.  This  is 
especially  true  at  low  temperatures  where  diffusion  is  more  important.  One 
can  also  conclude  that  the  cold  pressing  and  sintering  operations  have  not 
caused  extensive  alteration  of  the  individual  characteristics  of  the  powder 
particles. 

(4)  Inert  gas  effect 

When  an  increase  in  system  pressure  occurs  at  the  face  of  the  specimen, 
the  rate  at  which  hydrogen  will  flow  into  the  specimen  is  determined  by  the  gas 
pressure  gradient.  In  the  case  of  the  pure  hydrogen  gas  transients,  the  absorp¬ 
tion  of  hydrogen  by  the  specimen  results  in  a  continual  pressure  gradient.  When 
the  pressure  within  the  specimen  is  equal  to  the  system  pressure,  no  gas  flow 
will  occur  and  a  new  equilibrium  is  attained.  This  assumes  that  the  specimen 
temperature  is  also  in  equilibrium  with  its  surroundings. 

When  the  increase  in  system  pressure  was  from  the  addition  of  hydrogen  and 
an  inert  gas,  the  initial  response  followed  that  of  the  pure  hydrogen  case.  How¬ 
ever,  because  the  inert  gas  by  definition  did  not  react  with  the  ZrH,  no  inert  gas 
was  absorbed-,  This  resulted  in  a  buildup  of  gas  pressure  within  the  specimen. 
When  the  pressure  within  the  specimen,  which  was  primarily  inert  gas,  equaled 
the  system  pressure,  the  net  gas  flow  into  the  specimen  was  essentially  zero. 
However,  because  the  void  channels  within  the  specimen  contained  a  much 
higher  percentage  of  inert  gas  constituent  than  was  present  in  the  system  or  at 
the  face  of  the  specimen,  this  concentration  gradient  resulted  in  a  molecular 
diffusion  of  the  constituent  gases,  such  that  the  net  flow  of  hydrogen  was  into 
the  specimen  while  the  net  flow  of  inerts  was  out.  Thus,  a  continual  system 
pressure  decay  or  hydrogen  absorption  was  observed  but  at  a  much  slower  rate. 

The  experimental  data  for  the  inert  gas  test  are  given  in  table  XII  for  varying 
test  conditions  and  inert  gas  constituents. 

Examination  of  table  XII  shows  that  the  temperature,  pressure  and  particle 
size  effects  are  still  evident.  Superimposed  upon  these  effects  is  that  of  the 
inert  gas,  which  tends  to  suppress  hydrogen  absorption  for  the  reasons  stated. 
Referring  to  the  krypton  gas  data,  decreasing  the  initial  inert  gas  fraction  de¬ 
creases  the  inert  gas  effect  (see  figure  29).  The  results  illustrated  in  figure  20 
are  typical  for  almost  all  of  the  inert  gas  runs. 
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(1)  HYDRA  II  vs  experimental  data 

Of  major  concern  in  the  preparation  of  HYDRA  II  was  the  validity  of  the  mass 
transport  model.  As  such,  the  diffusion  and  flow  equations  are  elaborately  de¬ 
veloped  while  the  heat  transfer  equations  are  included  in  a  rather  simple  lumped 
model.  One  wouid,  therefore,  expect  HYDRA  II  to  most  accurately  describe  the 
transport  process  in  those  cases  where  heat  transfer  effects  are  not  controlling. 
In  other  cases,  it  was  clear  that  a  more  accurate  description  of  the  specimen 
temperature  distribution  was  necessary  to  improve  the  analytical  predictions. 
Such  was  found  to  be  the  case  in  the  present  study. 

The  input  parameters  which  most  heavily  influence  the  prediction  of  hydrogen 
absorption  were: 

d  permeability  of  the  body 

fp  shape  factor  for  diffusion 

L  heat  loss  term 

pCp  volumetric  heat  capacity 

The  value  of  d  was  taken  directly  from  the  results  of  the  permeability  tests.  The 
other  values  were  varied  within  their  probable  range  to  obtain  best  agreement 
with  the  data.  The  effect  of  each  is  discussed  in  the  following  paragraphs. 

(1)  The  particle  diffusion  shape  factor  (fp):  The  inverse  of  the  shape 
factor  was  employed  in  the  mass  diffusion  rate  equation.  Thus,  decreasing  its 
value  increased  the  rate  of  hydrogen  absorption  or  the  initial  slope  of  the  percent 
hydrogen  absorbed  versus  te3t  time  curve.  However,  a  point  was  reached  where 
decreasing  the  shape  factor  no  longer  affects  the  absorption  rate.  At  this  point, 
diffusion  was  essentially  instantaneous  and  the  absorption  rate  was  controlled  by 
thermal  effects.  In  a  perfectly  insulated  container,  the  absorption  curve  would 
be  flat  (slope  =  0)  after  a  few  seconds  because  the  absorption  process  continues 
only  if  a  specimen  heat  loss  exists. 

(2)  The  heat  loss  factor  (L):  The  value  of  the  system  heat  loss  term  con¬ 
trols  the  rate  at  which  the  specimen  exchanges  heat  with  its  environment.  Thus, 
its  effect  was  prevalent  throughout  the  transient.  In  the  present  version  of 
HYDRA  II,  the  heat  loss  was  taken  as  proportional  to  the  temperature  differential 
between  the  average  specimen  temperature  and  the  ambient  temperature  which 
was  assumed  to  be  the  initial  specimen  value.  Estimates  of  the  proportionality 
have  been  made.  The  primary  mode  of  system  heat  loss  was  thought  to  be  across 
the  hydrogen  gap  between  the  specimen  and  the  retort  wall.  As  a  result,  the  heat 
loss  factor  was  expected  to  be  on  the  order  of  the  hydrogen  conductivity  value 
(approximately  10  cal/cm-sec-°C). 
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(3)  Specimen  specific  heat  \pC)g:  The  specimen  specific  heat  controlled 
the  average  temperature  of  the  specimen  for  a  given  rate  of  heat  adoition  or  loss. 
Thus,  its  value,  along  with  the  system  heat  loss,  determined  the  point  in  the 
transient  where  the  absorption  curve  leveled  out  (slope  »  0).  In  the  current 
version  of  HYDRA  II,  the  value  for  the  specific  heat  is  a  constant  input  param¬ 
eter.  In  reality,  the  specific  heat  of  zirconium  hydride  is  a  function  of  tempera¬ 
ture  and  H/Zr  atom  ratio.  Thus,  the  specific  heat  is  more  accurately  taken  as 
a  time  variant  function. 

By  proper  selection  of  these  parameters,  it  is  possible  to  obtain  a  fair  agree¬ 
ment  between  HYDRA  II  and  the  experimental  data  (see  figures  21  and  22).  in 
general,  HYDRA  II  followed  the  transients  controlled  by  the  diffusion  process  such 
as  those  run  at  1150°  F  or  with  large  particle  sizes  quite  well.  Where  rates  were 
dominated  by  heat  transfer  considerations,  the  success  of  the  model  was  less 
dramatic,  as  shown  by  the  results  obtained  with  small  particle  sizes  at  1300°  or 
1400°  F. 

The  digital  code  has  been  converted  to  run  on  the  IBM  360  Model  65  system 
which  has  a  faster  running  time  and  larger  storage  capacity  than  the  IBM  7094. 

The  average  computing  time  for  a  60-second  transient  is  approximately  20 
minutes.  The  more  severe  a  transient,  the  longer  the  running  time.  This  re¬ 
sults  from  the  fact  that  for  a  transient  in  which  the  system  parameters  (pressure 
and  temperature)  are  continually  varying,  all  of  the  time  history  terms  in  the 
diffusion  equation  significantly  affect  the  absorption  rate  and  must  be  evaluated 
at  every  time  step. 

Efforts  to  reproduce  the  experimental  results  obtained  with  the  ir.ert  gas 
were  not  completely  successful.  The  code  predicts  much  slower  absorption 
rates  and  faster  inert  gas  pressure  buildup  in  the  specimen  voids  than  were 
experienced  in  the  experiment.  In  general,  predicted  hydrogen  absorption  was 
about  one -half  of  that  observed.  The  mathematical  model  used  does  not  account 
for  these  possibly  significant  processes: 

(1)  Adsorption  of  the  inert  gas  by  the  porous  body 

(2)  Molecular  diffusion  of  the  inert  gas  from  the  body 

(3)  Separation  of  the  gas  mixture  during  the  flow  process. 

These  and  other  causes  can  be  checked  analytically  and  experimentally  and 
should  the  inert  content  of  the  TURPS  fuel  element  gas  space  prove  to  be  large. 


53 


Hydrogen  Absorbed  (%) 


Hydrogen  Absorbed  (°}e) 


AFWL-TR-67-i  14 


Figure  22. 


Hydrogen  Absorption  Response  for  Nonspberical  ZrH 


x 


5  5 


AFWL-TR  -67  -  I  !4 


such  action  should  be  pursued.  Section  III  is  a  discussion  of  gaseous  fission 
product  release. 

c.  T^rmal  Cycling  Test 

The  thermal  cycling  test  data  are  presented  in  table  XIII.  The  minimum 
and  maximum  points  for  the  system  parameters  are  given  for  a  typical  eyrie 
in  the  test  history.  A  trace  of  the  variation  of  the  system  parameters  during 
these  typical  cycles  is  shown  in  figure  23. 

The  permeation  rate  of  hydrogen  through  a  metal  is  quite  significant  at  the 

3  2 

thermal  cycling  temperature  (approximately  5  cm  (STP)/cm  /hr /mil  at  780 

mm  H  ).  Thus,  a  hydrogen  loss  and  the  resulting  continual  hydrogen  pressure 

& 

decrease  throughout  the  course  of  the  test  was  observed.  It  was  also  noted  that 
setting  the  retort  heater  power  too  high  resulted  in  the  override  of  the  tempera¬ 
ture  controller  set  points.  Both  of  these  caused  a  variance  in  the  retort  hydro¬ 
gen  pressure  and  the  shape  of  the  system  parameter  trace  (see  figure  24).  Re¬ 
ferring  to  figure  24,  it  was  observed  that  the  leakage  of  hydrogen  from  the  retort 
lowered  the  hydrogen  pressure  to  the  point  where  the  specimen  (S6)  was  cycled 
through  the  /3-6  phase  equilibrium. 

Visual  examination  of  the  specimens  at  the  conclusion  of  the  thermal  cycling 
test  revealed  no  change  in  appearance.  As  previously  discussed,  "here  was  no 
discernible  difference  in  the  permeability  values.  The  cnly  difference  noted 
was  that  for  Specimen  N3,  the  hydrogen  absorption  rates  for  tlm  1 30u°  and  1150°  F 
initial  specimen  temperature  runs  were  higher  for  the  pretest  condition.  Examin¬ 
ation  of  the  rates  at  1400°  F  revealed  no  change  between  the  ^«re-  ar.d  post -test 
results.  The  cause  may  well  have  been  an  increase  in  soecimei  contamination 
during  the  long  cycling  run.  Such  an  effect  v/o*Td  be  most  evident  at  low  tem¬ 
perature. 

d.  Effect  of  Testing  on  Bodies 

As  previously  mentioned,  the  porous  bodies  were  prepared  by  cold  pressing 
and  sintering.  They  were  then  placed  in  a  stainless  steel  capsule  and  hydrided. 
The  density  change  on  hydriding  provided  a  very  tight  expansion  fit.  Figure  2  5 
shows  a  typical  body  prior  to  testing.  The  small  pitted  areas  were  caused  during 
handling  of  the  green  (unsintered)  body  ar.d  have  no  effect  on  the  results.  Fig¬ 
ures  26  and  27  show  the  tightness  of  the  fit  for  a  typical  cross  section.  The 
maximum  and  minimum  gaps  are  shown. 
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Figure  24.  Thermal  Cycle  Test 


1 


AFWL-TR-67-1 14 


Following  absorption  testing,  the  bodies  were  inspected  for  cracking  and 
other  effects.  Most  of  the  nonuranium-contaming  bodies  showed  no  effect. 

Three  such  bodies  are  shown  in  figure  28.  Some  showed  surfa'  ?-  .’  lamination 
apparently  from  oxygen  inadvertently  introduced  into  the  system.  Even  thcae 
were  intact  and  showed  only  the  surface  cracks  and  imperfections  that  were 
noted  prior  to  testing  (figure  29). 

Those  bodies  which  were  thermal-hydrogen  cycled  for  200  to  400  cycles 
appeared  free  from  physical  damage.  A  typical  body  is  shown  in  figure  30. 

There  was  no  visible  spallation  and  no  measurable  loss  of  material  from  these 
bodies. 

8.  Conclusions  and  Recommendations 

The  following  conclusions  and  recommendations  are  made: 

(1)  Within  the  limits  of  experimental  error,  the  hydrogen  absorption  test¬ 
ing  and  thermal -hydrogen  cycling  had  no  effect  on  the  permeabilities  of  the  bodies. 
It  is  concluded,  therefore,  that  the  basic  pore  size  and  structure  were  unaffected 
by  these  tests. 

(2)  The  permeability  measurements  show  that  previous  application  of  the 
Ergun  correlation  to  sintered  bodies  gave  conservative  values  in  that  predicted 
values  were  found  to  be  lower  than  experimental  values.  Shape  factors  for 
bodies  formed  from  sintered  spheres  indicate  the  behavior  is  very  close  to  that 
of  loosely  packed  spheres.  It  is,  however,  recommended  that  the  permeability 
of  all  future  porous  bodies  be  experimentally  determined  to  minimize  analytic 
uncertainties.  This  appears  necessary  because  of  the  wide  variations  that  can 
be  introduced  by  varying  the  pressing  and  sintering  condition. 

(3)  The  permeability  of  the  porous  bodies  was  found  to  be  independent  of 
the  gas  used  or  its  pressure  in  the  range  investigated.  Flow  through  the  bodies 
was  laminar  and  obeyed  Darcy’ s  law. 

(4)  Expected  temperature,  pressure  and  body  characteristic  effects  were 
observed  in  the  hydrogen  absorption  tests.  Because  of  sharply  increasing  diffu¬ 
sion  coefficients,  absorption  rates  were  found  to  increase  rapidly  with  tempera¬ 
ture.  Varying  the  initial  pressure  of  the  transient  resulted  in  changing  the  initial 
driving  function  (concentration  change).  The  available  concentration  change  in¬ 
creases  as  the  initial  pressure  is  lowered,  thus  resulting  in  a  higher  rate  of 
absorption. 
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The  inert  gas  results  exhibited  higher  absorption  rates  than  anticipated, 
although  the  net  effec*  of  the  inerts  was  to  reduce?  the  hydrogen  absorption,  at 
least  in  the  short-term  transient.  There  appeared  to  be  no  basic  difference 
oetween  the  various  gases  used,  i.e.  ,  argon,  krypton  and  xenon. 

(5)  Good  analytical  agreement  with  the  absorption  experiments  was 
obtained  where  the  absorption  rate  was  controlled  by  transport  processes. 
Poorer  agreement  was  obtained  in  those  cases  where  heat  transfer  processes 
were  rate  controlling.  Such  results  are  understandable  in  light  oi  the  rather 
simple  lumped  heat  transfer  model  being  used.  The  results  confirm  the  validity 
of  the  mass  transport  model.  A  mere  complete  thermal  model  should  be  used 
in  future  TURPS  analysis  where  appropriate. 

(6)  With  the  exception  of  some  surface  contamination  on  <x  few  of  the 
samples,  the  testing  had  no  effect  on  the  physical  appearance  of  the  porous 
bodies.  This  included  bodies  that  were  thermal -hydrogen  cycled  tor  up  to  105 
cycles.  No  physical  damage  of  any  kind  wa«  attributed  to  the  testing. 
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SECTION  III 

THERMAL  CONDUCTIVITY  OF  POROUS  ZIRCONIUM  HYDRIDE 


1.  Purpose 

The  thermal  conductivity  of  the  porous  fuel  body  and  reservoir  matrix  of 
TURPS  is  an  important  system  characteristic.  Because  maximum  fuel  temper¬ 
atures  are  major  design  constraints,  it  is  important  to  have  accurate  knowledge 
of  thermal  conductivity  to  enable  optimum  design.  Prior  to  this  investigation, 
the  Bruggeman  variable  dispersion  model  for  porous  matrices  was  used. 

From  this  model,  knowledge  of  the  solid  material  and  gas  conductivities  yields 
a  prediction  of  the  matrix  conductivity  in  the  following  equation: 


where  the  subscripts  1,  2  and  m  refer  to  the  solid,  gas  and  composite  matrix, 
respectively,  and  «  is  the  void  or  gas  volume  fraction.  The  objective  of  this 
study  was  to  test  the  validity  of  this  model  for  TURPS  sintered  matrix  material 
and  to  provide  empirical  data  for  reactor  design  use. 

The  technique  chosen  for  evaluation  of  thermal  conductivity  was  the  flash 
method  (reference  1)  of  obtaining  thermal  diffusivity  together  with  separate 
evaluation  of  specific  heat.  This  method  provided  a  means  of  avoiding  trouble¬ 
some  heat  loss  considerations  encountered  in  conventional  steady-state  methods 
and  elimination  of  gross  hydrogen  concentration  gradients  from  impressed  tem¬ 
perature  gradients. 

It  was  recognized,  however,  that  the  thermal  diffusivity  measurements 
might  be  impaired  by  hydrogen  dissociation  upon  impingement  of  the  heat  source 
(light  flash)  and  therefore  the  experiments  were  conducted  in  two  stages.  The 
first  stage  used  unhydrided,  sintered  zirconium  bodies  composed  of  differing 
spherical  particle  size  distribution  to  yield  a  range  of  matrix  densities.  Also, 
the  void  space  was  varied  to  contain  argon,  helium  and  vacuum  to  enable  deter¬ 
mination  of  the  gas  contribution  to  thermal  conductivity.  The  second  stage  used 
hydrided,  sintered  zirconium  bodies  in  which  the  effects  of  vacuum  and  argon  gas 
presence  in  the  void  could  be  examined  up  to  650°  F,  a  temperature  range  in 
which  the  hydrogen  dissociation  pressure  is  negligible.  Following  this,  a  hydro¬ 
gen  atmosphere  was  introduced  to  investigate  the  behavior  cf  the  hydride  material 
at  higher  temperatures. 
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The  thermal  diffusivity  is  related  to  thermal  conductivity,  specific  heat  and 
density  by  the  following  equation: 

a  (thermal  diffusivity)  = 

P 

where: 

k  =  thermal  conductivity 
p  =  density 

=  specific  heat 

in  this  investigation,  the  thermal  diffusivity  was  obtained  from  the  time  required 
for  a  temperature  pulse  to  reach  one -ha1  f  of  its  maximum  value  (asymptotic 
maximum  in  a  perfectly  insulated  body)  at  the  rear  surface  of  a  slab  whose  front 
surface  receives  an  impulse  of  radiant  energy.  According  to  Parker  (reference 
1),  if  the  energy  pulse  duration  is  short  compared  to  the  temperature  rise  time, 
and  if  heat  losses  are  negligible,  the  lollowing  relationship  is  applicable: 

1-37  L2 
*2- 
T 1  /  2 

where  L  is  the  slab  thickness  and  T  is  the  time  f°r  the  temperature  to  reach 
one -half  of  its  maximum  rear  surface  value. 

2.  Experiment  Description 

All  samples  tested  were  made  of  sintered,  spherical  zirconium  particles. 

The  particle  size  distributions  used  id  the  dimensions  of  each  specimen  are 
given  in  table  XIV.  The  densities  of  both  hydrided  and  unhydrjded  bodies  were 
calculated  from  the  post-sintering  dimensions  and  weight  (before  hydriding). 

It  is  of  interest  to  aote  that  although  the  sintering  was  done  at  two  conditions, 

900°  C  for  two  hours  and  1100°  C  for  four  hours,  only  those  bodies  sintered  at 
the  latter  condition  possessed  suitable  strength  for  handling  during  the  dif¬ 
fusivity  measurements. 

After  unsuccessful  attempts  at  using  pressure  contact  for  the  chromel- 
alumel  thermocouple  against  the  rear  surface  of  the  test  specimens,  a  spot 
welding  technique  was  developed.  Because  of  the  porous  nature  of  the  surfaces, 
spot  welding  was  difficult,  particularly  for  the  hydride  material.  Furthermore, 
it  was  found  essential  to  weld  each  wire  separately  to  the  surface  so  that  a  bead- 
induced  temperature  lag  would  b<-  avoided.  It  was  during  the  spot  welding  develop¬ 
ment  that  the  900°  C  sintered  bodies  were  deemed  not  suitable  for  testing.  Inter - 
particle  adherence  was  insufficient  to  allow  a  suitably  strong  bend  between  the 
thermocouple  wires  and  the  surface. 
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Corre  ctions  for  heat  loss  from  the  test  specimens  were  made  according 
to  a  one -dimensional  heat  transfer  model  (references  1  and  2).  Because  the  heat 
loss  was  nearly  equal  for  specimens  tested  in  vacuum  and  in  argon  or  helium 
atmospheres,  it  was  concluded  that  the  main  path  was  conduction  to  the  insu¬ 
lator  and  the  thermocouple  wire  holders  pressed  against  the  rear  surface 
of  the  pellet.  The  effect,  was  small,  however,  causing  only  about  a  four  per¬ 
cent  decrease  in  the  observed  response  time. 

The  accuracy  and  reliability  of  this  thermal  conductivity  determination 
method  was  estimated  by  comparison  with  a  known  material  (Zircaloy-2  speci¬ 
men)  and  by  checking  the  reproducibility  of  measurements.  Reproducibility 
differences  of  less  than  five  percent  and  Zircaloy-2  results,  in  agreement 
within  10  percent  of  accepted  literature  values,  indicate  an  overall  accuracy 
of  about  10  percent. 

3.  Resuits 

a.  Unhydrided  Zirconium 

A  summary  of  the  porous  zirconium  body  test  results  is  presented  in 
table  XV.  The  measured  thermal  diffusivities  at  various  temperatures  for 
vacuum,  argon  and  helium  gas  presence  are  listed  for  each  specimen  studied. 
The  values  of  thermal  conductivity  were  derived  with  the  aid  of  zirconium  heat 
capacity  data  of  figure  31.  Figure  32  contains  gas  conductivities  for  use  in  the 
analytical  prediction  of  porous  body  thermal  conductivity  by  the  Bruggeman 
model. 

The  experimentally  determined  thermal  conductivities  for  76,  78, and  80 
percent  density  zirconium  matrices  are  presented  graphically  in  figure  33  to¬ 
gether  with  analytic  predictions  with  helium  and  argon  gas  environments  (one 
atmosphere  pressure).  Also  included  is  a  best  literature  result  for  pure  solid 
(100  percent  dense)  zirconium.  It  is  apparent  that  the  analytic  model  grossly 
underestimates  the  absolute  value  of  matrix  conductivity,  while  overestimating 
the  gas  effect.  This  failure  is  readily  understandable  in  that  the  model  assumes 
point  contact  of  spherical  particles  whereas  photomicrographs  of  the  porous 
bodies  show  a  clear  connection  between  the  particles.  The  effectiveness  of 
this  particle  interconnection  is  illustrated  by  the  high  values  of  conductivity 
in  vacuum  environment. 

The  observed  sensitivity  of  thermal  conductivity  with  body  void  fraction 
(or  density)  is  shown  in  figure  34.  It  should  be  noted  that  the  gross  parameter 
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(void  fraction)  is  dependent  upon  the  particle  size  distribution  and  the  sintering 
process  used  i.i  specimen  fabrication.  Therefore,  this  graph  shoula  be  inter¬ 
preted  cautiously  until  further  study  of  microscopic  structure  effects  is  per¬ 
formed. 


Table  XIV 

SPECIMEN  FABRICATION  DATA 


Specimen 

number 

Composition 
(particle  grading) 

Length 
(in. ) 

Diameter 
(in. ) 

Weight 

(gm) 

Density 

% 

X 

Unhydrided  Zirconium 

6-9  70%,  -170  +270  Mesh 

20%,  -80  +270  Mesh 
10%,  -400  Mesh 

0.  071 

0.  481 

1.0482 

76.  2 

6-7  70%,  -400  Mesh 

30 %,  -230  +325  Mesh 

0.  072 

0.  470 

1.  0388 

78.  1 

— 

3-7  80%,  -170  +270  Mesh 

20%,  -400  Mesh 

0.  074 

0.  494 

1. 0863 

72.  0 

— 

3-5 

100%,  -400  Mesh  0. 065  0.  479  0,  9975 

Hydrided  Zirconium  ZrHx 

80.  0 

6-3 

75%,  -100  +120  Mesh 
25%,  -80  +270  Mesh 

0.  056 

0.  493 

0.  8316 

73.  1 

1.  59 

6-8 

70%,  -400  Mesh 

30%,  -230  +325  Mesh 

0.  047 

0.  471 

0.  6850 

78.  6 

1.  59 

6-10 

70%,  -170  +270  Mesh 

20%,  -230  +325  Mesh 
10%,  -400  Mesh 

0.  071 

0.  480 

1. 1100 

79.  2 

1.  78 

3-6 

100%,  -400  Mesh 

0.  070 

0.  479 

1.  0732 

79.  7 

1.  59 

3-8  i 

80%,  -170  +270  Mesh 

20%,  -400  Mesh 

0.074 

0.  494 

1.  0863 

72.  4 

1. 78 
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Table  XV 

EXPERIMENTAL  RESULTS 


Temperature  (°  F)  ! 

Gas 

I  o  (cm2/ sec) 

k  {Btu/hr-ft-#  F) 

Specimen  6-S;  p 

=  76,  4  percent 

70 

Vacuum 

0.  044 

in 

CO 

330 

Vacuum 

0.  045 

3.9 

730 

Vacuum 

0.  044 

4.  2 

1080 

Vacuum 

0.  045 

4.  4 

70 

Argon 

0.  046 

3.  6 

300 

A.rgon 

0.  045 

3.9 

710 

Argon 

0.  045 

4.  2 

70 

Helium 

0.  049 

3.  8 

250 

Helium 

0.  049 

4.2 

600 

Helium 

0.  048 

4.4 

1040 

Helium 

0.  049 

4.8 

Specimen  6-7;  p  =  78.  1  percent 

70 

icuum 

0.  058 

4.7 

350 

Vacuum 

0.  058 

5.  2 

690 

Vacuum 

0.  054 

5.2 

1112 

Vacuum 

0.  054 

5.  5 

70 

Argon 

0.  058 

4.7 

325 

Argon 

0.  056 

5.  1 

660 

Argon 

0.  054 

5.  2 

Argon 

0.  057 

5.  7 

70 

Helium 

0.  061 

5.  0 

250 

Helium 

0.  059 

5.  2 

450 

Helium 

0.  058 

5.  4 

900 

Helium 

0.  058 

5.7 
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Table  XV  (cont'd) 


Temperature  (°  F) 

Gas 

IBSfBfR 

k  (But/hr-ft-0  F) 

Specimen  3-7;  p 

8  72  percent 

70 

Vacuum 

0.  051 

3.  8 

350 

Vacuum 

0.  047 

4.0 

710 

Vacuum 

0.  045 

4.  0 

1050 

Vacuum 

0.  046 

4.  2 

70 

Argon 

0.  051 

3.  8 

300 

Argon 

0.  048 

3.9 

680 

Argon 

0.  047 

4.  2 

1010 

Argon 

0.  046 

4.  2 

70 

Helium 

0.  052 

3.  9 

300 

Helium 

0.  050 

4.  1 

550 

Helium 

0.  051 

4.  5 

900 

Helium 

0.  050 

4.  5 

Specimen  3-5; 

;  p  =  80  percent 

70 

Vacuum 

0.  061 

6.3 

275 

Vacuum 

0.  064 

7.  2 

700 

Vacuum 

0.  059 

7.3 

70 

Argon 

0.  064 

6.  6 

300 

Argon 

0.  062 

7.  1 

660 

Argon 

0.  058 

7.  1 

70 

Helium 

0.  065 

6.  7 

190 

Helium 

0.  066 

7.3 

460 

Helium 

0.  059 

7.0 

b.  Hydrided  Zirconium 


The  diffusivity  measurement  results  and  the  derived  values  of  thermal  con¬ 
ductivity  for  the  hydrided  porous  zirconium  matrix  specimens  are  contained  in 
table  XVI.  The  heat  capacity  data  of  Beck  as  presented  in  AFWL-TR -66- 16  and 
reproduced  in  figure  35  was  verified  over  a  limited  range  of  temperatures  by 
calorimetry  measurements  on  specimens  of  zirconium  hydride  prepared  simul¬ 
taneously  with  the  diffusivity  test  bodies. 
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Figure  33.  Conductivity  of  Porous  Sintered  Zirconium  Bodies 
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Table  XVI 

EXPERIMENTAL  RESULTS  FOR  HYDRIDED  SPECIMENS 


Temperature  (°  F) 


Gas 


a  (cm  / sec) 


k  tBtu/hr  -ft-0  F) 


Specimens  6-3;  p  -  73.  1  percent;  x  r  1.  59 


70 

Vacuum 

0.  038 

3.  5 

305 

Vacuum 

0.  033 

3.  7 

655 

Vacuum 

0.  030 

4.  2 

70 

Argon  ( 1  atm) 

0.  040 

3.  6 

300 

Argon  (I  atm) 

0.  035 

3.9 

650 

Argon  ( 1  atm) 

0.  03  2 

4.  4 

Specimen  3-6;  p  =  79.  7  percent;  x  r  1.  59 

70 

Vacuum 

0.  047 

4.  6 

345 

Vacuum 

0.  043 

5.  2 

690 

Vacuum 

0.  038 

5.  9 

70 

Argon  (1  atm) 

0.  050 

4.  9 

330 

Argon  (1  atm) 

0.  044 

5.  3 

700 

Argon  (1  atm) 

0.  041 

6.  3 

Specimen  3-8;  p  -  72.  4  percent,  x  ;  1.  78 

70 

Vacuum 

0.  040 

4.  0 

290 

Vacuum 

0.  038 

4.  6 

640 

Vacuum 

0.  032 

4.  7 

70 

Argon  (1  atm) 

0.  044 

4.  5 

280 

Argon  (1  atm) 

0.  040 

4.  7 

640 

Argon  (1  atm) 

0.  035 

5.  2 

Specimen  6-10;  p-  79. 

2  percent,  x  :  1.  78 

70 

Vacuum 

0.  044 

5.  0 

300 

Vacuum 

0.  043 

5.  7 

600 

Vacuum 

0.  038 

5.  9 

800 

Vacuum 

0.  035 

6.  2 

900 

Vacuum 

0.  034 

6.  3 

800 

H2  (1  atm) 

0.  043 

7.  5 

Specimen  6-8.  P  -  78.  6 

percem,  x  '  1 .  59 

70 

Vacuum 

0.  038 

3.  8 

305 

Vacuum 

0.  037 

4.  4 

625 

Vacuum 

U.  03  4 

5.  0 

70 

OS’s  Ar  -  5°o  H  (1  atm) 

0.  04.) 

4.  2 

3  00 

95tT=  Ar  -  5°c  H  (1  atm) 

0.  040 

4.  8 

600 

95^  Ar  -  o’c  H  (1  atm) 

0.  038 

5.  5 

930 

95<rc  Ar  -  5rc  H  (1  atm) 

0.  036 

6.  5 
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(1)  Low  teniperature  evaluation  (no  significant  hydrogen  dissociation 
pressure) 

The  thermal  conductivity  of  porous  bodies  of  about  73  and  80  percent  of 
maximum  density  is  compared  in  figure  36  to  solid  hydride  material  conductivity 
obtained  by  other  investigators.  It  is  apparent  that  without  the  presence  of  gas 
in  the  void  space,  the  porous  bodies  possess  about  50  to  60  percent  of  the  solid 
material  conductivity,  a  result  similar  to  that  obtained  for  the  unhydric'ad  zir¬ 
conium.  Likewise,  the  tabulated  data  show  that  argon  gas  exerts  a  small  effect, 
less  than  10  percent  on  the  body  conductivity. 

The  relationship  between  conductivity  and  void  fraction  for  the  hydrided 
bodies  is  weaker  than  that  observed  for  the  unhydrided  zirconium.  This  small 
dependence,  illustrated  in  figure  37,  is  probably  a  result  of  increased  particle 
contact  during  the  hydriding  process. 

(2)  High  temperature  evaluation 

Measurements  of  thermal  aiffusivity  at  temperatures  where  the  hydrogen 
dissociation  pressure  is  significant  were  undertaken  with  a  mixed  argon-hydrogen 
atmosphere  and  with  a  pure  hydrogen  atmosphere.  During  the  energy  pulse  tran¬ 
sient,  the  imposed  partial  hydrogen  pressure  was  maintained  at  the  equilibrium 
value  for  the  corresponding  initial  temperature  and  hydrogen  content  of  the  hydride. 
Because  of  the  small  sample  temperature  rise  incurred  during  the  transient 
(<  5°  F),  no  adjustment  was  deemed  necessary. 

The  significant  finding  of  these  experiments  was  a  marked  reduction  in  the 
maximum  temperature  rise  recorded.  This  observation  can  be  explained  by 

«  > 

the  dissociation  of  hydride  in  the  specimen  upon  receipt  of  the  radiant  energy 
impulse.  Unfortunately,  because  of  the  resultant  small  noisy  signal,  detailed 
analysis  was  not  possible.  Future  efforts  are  required  to  yield  data  amenable 
to  interpretation  by  application  of  hydrogen  transport  analysis. 

4.  Conclusions  and  Recommendations 

As  indicated  by  the  data,  the  Bruggeman  model  failed  to  predict  accurate 
thermal  conductivity  of  the  porous  zirconium  or  zirconium  hydride  matrix.  One 
may  conclude  that  the  sintering  operation  yields  sufficient  interparticle  connec¬ 
tion  and  that  the  presence  of  gas  in  the  void  space  is  not  necessary  to  yield  good 
thermal  conductivity.  This  result  adds  further  emphasis  to  the  importance  of 
thorough  sintering  of  the  particle  matrix.  Not  only  is  the  body  strength  increased, 
but  the  thermal  conductance  is  greatly  enhanced.  From  the  results  of  hydride 
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conductivities  made  at  low  temper  ature,  together  with  the  effect  of  hydrogen 
gas  on  the  conductivity  of  specimen  6-10  at  800°  F  (a  20  percent  increase),  it 
is  reasonable  to  expect  thermal  conductivities  of  7  to  9  Btu/hr-ft-°F  at  1500 
F  for  porous  bodies  of  70  to  80  percent  density. 

In  view  of  the  desirability  for  measurements  at  higher  temperatures  where 
significant  hydrogen  pressures  exist,  together  with  an  investigation  of  the  im¬ 
portant  effects  of  the  sintering  process,  additional  thermal  conductivity  studies 
are  recommended.  These  efforts  should  encompass  a  range  of  specimen  fabri¬ 
cation  techniques  (particle  size  and  sintering)  and  should  emphasize  testing  at 
TURPS  operating  conditions.  Also,  analysis  of  the  thermal-hydrogen  transport 
phenomena  by  adaption  of  the  techniques  presented  in  Section  II  of  this  report 
will  enable  better  application  of  the  thermal  diffusivity  method  to  this  task. 
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SECTION  IV 

CHEMICAL  AND  PHYSICAL  PROPERTIES  OF  THEON  (P^) 

1 .  Purpose 

It  was  the  purpose  of  this  task  to  establish,  by  means  of  analysis  and/or 
measurement,  the  specific  characteristics  of  Theon  required  to  permit  future 
TURPS  reactor  core  designs  using  this  material  as  the  coolant.  The  necessity 
for  undertaking  such  an  investigation  stems  solely  from  the  lack  of  such  informa¬ 
tion  in  the  existing  literature. 

The  particular  characteristics  or  properties  of  interest  are  as  follows: 

(1)  Viscosity 

(2)  Surface  tension 

(3)  Thermal  conductivity 

(4)  Vapor  pressure 

(5)  Liquid  density 

(6)  Heat  of  fusion 

In  performing  the  investigation,  two  major  areas  of  activity  were  required. 
The  first  concerned  the  preparation  of  purified  Theon  to  supply  the  require¬ 
ments  for  property  evaluation  and  other  tasks  such  as  the  Theon  irradiation 
and  chemical  compatibility  studies.  The  second  major  area  concerned  the  de¬ 
sign  and  performance  of  the  various  physical  property  measurement  experiments. 

2.  Preparation  of  Theon 

a.  Preparation  of  Purified  Stoichiometric  Theon  (P^S^) 

Commercial  grade  P^S^  contains  trace  amounts  of  elements  introduced  during 
the  svnthesis  procedure  and,  supposedly,  a  slight  excess  of  free  sulfur  added  to 
increase  the  stability  of  the  material.  Removal  of  these  impurities,  along  with 
any  of  the  higher  homologues  of  phosphorus  and  sulfur  (P4S-,  P4S7>  P4Sjq.  .  . ), 
is  desirable  to  lessen  the  likelihood  of  significant  changes  in  the  physical  and 
chemical  properties  of  the  reactor  coolant.  The  preparation  pr  uecure  con¬ 
stituted  the  present  standard  production  procedure  for  pure  Commercial 

grade  P^^  was  purified  and  converted  to  stoichiometric  P^S^  through  a  series 
of  solvent  extraction  and  distillation  steps  as  follows: 


1 
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(1)  The  starting  material  was  heated  for  two  hours  at  250°  C  under  an 
argon  atmosphere.  Sulfur  and  the  higher  sulfides  of  phosphorus  are  not 
stable  toward  phosphorus  when  melted,  thus  only  ?4S3  and  P  remained  in  the 
melt.  The  addition  of  excess  phosphorus  was  not  necessary  because  chemical 
analyses  indicated  an  excess  of  phosphorus  already  present  in  the  as -received 
material. 

(2)  The  molten  material,  consisting  of  P4S3,  red  phosphorus,  white 
phosphorus,  and  inorganic  impurities,  was  allowed  to  cool  to  room  temperature 
with  the  flow  of  argon  continued. 

(3)  Dissolution  of  the  solid  cake  was  accomplished  by  the  addition  of  CSg 
to  the  reaction  vessel.  The  CSg  dissolved  the  P4S3  and  any  white  phosphorus 
present,  but  left  the  red  phosphorus  and  other  solid  impurities  behind. 

(4)  After  filtration  to  remove  suspended  solids,  the  solution  was  cooled 
to  dry-ice  temperature  where  the  P4S3  and  white  phosphorus  separated  as 
fine  crystals.  The  crystals  were  filtered  off,  washed  with  ether  and  vacuum 
dried. 

(5)  The  removal  of  white  phosphorus  was  accomplished  through  vacuum 
distillation  where  the  purified  P4S3  was  heated  to  360°  C  under  vacuum.  The 
white  phosphorus,  being  more  volatile,  came  off  in  the  first  fraction  and  was 
discarded.  The  remaining  fraction,  consisting  of  P4S3  free  of  white  phosphorus, 
was  cooled,  dissolved,  recrystallized,  filtered,  washed  and  dried  as  in  step  4. 

(6)  Confirmation  of  the  complete  removal  of  phosphorus  was  afforded  by 
heating  samples  of  purified  P4S3  to  approximately  500°  C  for  one  hour  and  noting 
the  absence  of  red  or  white  phosphorus  in  the  condensed  vapor. 

b.  Preparation  of  Nonstoichiometric  Theon  (p4S3> 

(1)  Preparation 

Theon  forms  containing  five  percent  excess  sulfur,  one  percent  excess 
sulfur  and  five  percent  excess  phosphorus  were  prepared  for  use  in  the  static- 
corrosion  studies.  The  compositions  were  formulated  by  the  addition  of  pure 
sulfur,  and  phosphorus  to  stoichiometric  P4S3.  Wet  chemical  analyses  were 
performed  to  confirm  the  excess  sulfur  and  phosphorus  concentrations. 

(2)  Verification  of  Theon  composition 

A  weighed  sample  of  Theon  was  placed  in  an  Erlenmeyer  flask  and  reacted 
with  red  fuming  HNC>3  under  an  argon  atmosphere  in  an  ice  bath.  The  phosphor¬ 
us  content  was  thereby  oxidized  to  phosphoric  acid  and  the  sulfur  was  oxidized 
to  sulfuric  acid. 
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The  phosphorus  was  precipitated  with  a  magnesia  mixture,  collected, 
calcined  at  850°  C  and  weighed  as  Mg2P20^. 

The  sulfur  was  then  converted  to  SOg,  which  passed  through  an  iodine  solu¬ 
tion  and  was  titrated  with  standard  thiosulfate  solution.  Sulfur  content  evalua¬ 
tion  was  less  consistent  than  the  phosphorus  determination  and  repeated 
analyses  were  necessary. 

This  method  yielded  9S.  8  percent  accountability  of  stoichiometric  con¬ 

stituents  (P  and  S)  and  thus  enabled  confidence  in  the  analysis  of  excess  sulfur 
and  phosphorus  systems.  For  comparison,  it  should  be  noted  that  as-received 
P^Sg  showed  a  slight  excess  of  phosphorus  and  approximately  2.3  percent  other 
impurities. 

c.  Produce  Specification  Testing  for  Stoichiometric  Theon  (P^S^) 

Procedures  for  verifying  the  purity  of  stoichiometric  include  detection 

of  free  phosphorus,  emission  spectroscopy  and  wet-chemical  analysis.  Infrared 
spectroscopy  has  been  used  to  identify  P^S^,  but  adequate  sensitivity  for  impurity 
inspection  has  not  been  obtained  yet. 

(1)  Detection  of  free  phosphorus 

Free  phosphorus  is  detected  by  heating  P^S.^  to  500°  C  at  the  bottom  of  a 
Vycor  glass  tube  and  observing  the  vapor  condensate  in  the  cool  upper  tube 
region.  Concentrations  of  a  few  parts  per  million  will  yield  the  typical  red 
phosphorus  color.  White  phosphorus  is  observed  by  breaking  the  tube  after 
it  has  cooled  to  room  temperature.  Spontaneous  ignition  or  fuming  indicates 
free  white  phosphorus. 

(2)  Emission  spectroscopy 

Small  quantities  of  impurities  can  be  detected  by  emission  spectroscopy 
techniques.  Theon  material  purified  by  the  method  described  previously  showed 
about  50  pnm  silicon  content.  Less  than  1  ppm  was  observed  for  any  other  im¬ 
purity  elements  detected  (Al,  Ca,  Cu,  Zr,  Fe). 

( 3)  Wet  chemistry 

As  described  previously,  wet  chemical  determination  of  phosphorus  and 
sulfur  content  will  assure  Theon  content  of  these  components  to  within  0.2 
percent.  However,  present  techniques  have  not  been  directed  toward  impurity 
detection  by  this  means. 
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3.  Physical  Properties  Measurements 

Measurements  of  viscosity,  surface  tension, and  thermal  conductivity  were 
made.  In  addition,  vapor  pressure  and  liquid  density  values  were  obtained  for 
P^Sj  conditions  approximating  TURPS  operation.  Of  further  interest  was  the 
heat  of  fusion;  the  heat  content  measurements  about  the  normal  melting  point 
yielded  this  quantity. 

a.  Viscosity 

To  eliminate  the  necessity  for  flow  corrections,  a  comparative  method  was 
used  with  a  liquid  of  known,  comparable  viscosity  over  the  temperature  range 
of  interest.  Preliminary  experiments  showed  that  bismuth  was  a  suitable  com¬ 
parator,  and  its  viscosity  was  available  from  references  8,  9  and  10.  Measure¬ 
ments  of  the  effluent  times  for  bismuth  and  P4S3  in  the  same  reservoir-capillary 
configuration  (figure  38)  at  each  temperature  point  of  interest  enabled  P4S3 
viscosity  determination  from  the  following  equation  which  is  a  simplification  of 
an  approach  suggested  by  Kaaschow  (reference  11). 


0* 


P  t 
Ko  o 


where: 

n 


1  +  k 


(i?) 


*  required  viscosity  of  P4S3 

nQ  *  known  viscosity  of  bismuth 

p,  pQ  *  respective  densities  of  P4S3  and  bismuth 

t,  t  *  respective  effluent  times  of  P.S,  and  bismuth 
O  a  4  0 

k  .  mv 

„2  'g  r4 

m  ■  a  calibration  constant 

V  *  volume  of  P4S3 

I  *  capillary  length  =  mean  pressure  head 

r  *  capillary  radius 

g  *  gravitational  acceleration 


Two  capillary-reservoir  systems,  A  and  B,  and  two  independently  prepared 
P4S3  samples,  1  and  2,  were  employed  in  this  experiment.  The  result'  of  sys¬ 
tem  A  with  samples  1  and  2,  and  system  B  with  sample  2  are  listed  in  table 
XVR.  These  values  are  plotted  in  figure  39  and  the  best  fit  experimental  line 
is  compared  to  the  analytically  predicted  curve.  It  is  seen  that  the  measured 
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variation  with  temperature  is  very  close  to  the  predicted  behavior  although  the 
absolute  value  prediction  is,  as  expected,  only  a  rough  approximation  of  the 
actual  observed  values. 


Table  XVR 

P.S,  EXPERIMENTAL  VISCOSITY 
4  3 


1/T  (°R_1) 

Temperature 

<#F) 

Test 

Viscosity  (centipoise) 

A  1 

La  2 

1  B  2 

0.943 

600 

1.160 

1.261 

0.909 

640 

1.124 

1.135 

1.057 

0.877 

680 

0.824 

0.905 

0.865 

696 

-- 

0.82? 

0.847 

720 

0.531 

0.820 

760 

0.491 

-- 

0.810 

775 

-- 

0.601 

0.794 

800 

0.485 

0.490 

0.772 

835 

0.534 

0.769 

840 

-- 

0.489 

0.763 

850 

0.  ?92 

0.735 

900 

0.285 

0.425 

0.395 

0.712 

945 

0.437 

0.398 

b.  Surface  Tension 

Surface  tension  measurements  of  P^S^  to  510®  C  require  a  pressurized  system 
adapted  to  the  particular  physical  and  chemical  characteristics.  The  method  used 
in  this  investigation  was  the  capillary-rise  technique.  The  simple  apparatus  shown 
in  fig..  ~e  40  provided  the  basic  experimental  tool.  Careful  calibration  with  known 
materials  assured  sufficient  accuracy. 

2 

The  governing  liquid  rise  relationship  is:  2wrycos0  =  tr  hpg 
where 

r  a  radius  of  the  capillary' 

h  *  height  of  liquid  in  capillary  above  the  surface  of  liquid  in  which 
capillary  is  immersed 
®  =  angle  of  contact  of  liquid  and  capillary 

p  *  density  of  liquid 
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g  3  gravitational  acceleration 
y  3  surface  tension  (dynes/ cm) 

The  values  of  h  and  0  are  obtained  by  measurement  of  the  P^S.^  system,  and 
the  radius  of  the  capillary  is  carefully  determined  prior  to  use. 

The  average  diameter  of  the  capillary  is  determined  by  classical  mercury 
displacement  methods.  Tubing  lengths  of  best  diameter  uniformity  are  found  by 
examining  with  a  Kodak  optical  comparator  and  are  used  in  the  regions  of 
capillary  rise  measurement. 

In  preparing  the  surface  tension  unit,  the  P^S^  was  admitted  through  the  top, 
the  assembly  evacuated  and  the  Vycor  cylinder  tapered  and  sealed.  The  quantity 
of  was  governed  by  the  known  diameter  characteristics  of  the  im’-viduai 

capillary  tube. 

The  sealed  surface  tension  unit  was  then  placed  in  a  controlled  temperature 
furnace  and  subjected  to  the  temperature  range  of  interest  while  measurements 
of  the  capillary  rise  and  liquid  contact  angle  were  made  with  an  optical  level. 
Figure  41  shows  the  arrangement  of  this  experiment.  Figure  42  presents  the 
experimental  data  together  with  a  straight  line  fit  whose  equation  is 

y  (dynes/cm)  3  81.5  -  0.398  T  (T=#R) 

It  is  of  interest  to  note  that  extrapolation  of  the  data  to  zero  viscosity  suggests 
a  critical  temperature  of  about  1590*  F,  whereas  the  theoretical  prediction  was 
about  I860*  F. 

c»  Thermal  Conductivity 

'Pie  thermal  conductivity  of  P^S^  was  measured  by  the  flash  method  of  obtain¬ 
ing  thermal  diffusivity.  Separate  determination  of  heat  capacity  and  density 
enabled  evaluation  of  the  thermal  conductivity  from  the  following  relationship 

k  3  oCpp 

where 

k  2  thermal  conductivity 

a  •*  measured  diffusivity 

C  *  measured  heat  capacity 

p 

y  *  measured  density 

Table  XVGI  contains  the  measured  quantities,  while  figure  4H  presents  the 
thermal  conductivity  variation  with  temperature.  It  is  cf  interest  to  note  ’hat 

sr, 


Seal  point 


Figure  40.  Sj  Surface  Tension  Unit  Figure  41.  Surface  Tension  Measurement  Apparatus 
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the  conductivity  at  950°  F  is  0.45  Btu/hr-ft-°F,  significantly  greater  than  the 
prior  predicted  value  of  0.09  Btu/hr -ft-9F.  Such  a  value  invites  recognition  of 
aa  a  promising  forced  circulation  coolant  with  a  Prandtl  number  at  950°  F 

close  to  that  of  water  (  1.0). 


Table  XVIII 

THERMAL  CONDUCTIVITY  DETERMINATION  OF 

4  3 


Temperature 

Density 

3 

Heat  capacity 

Diffusivity 

2 

Thermal 

conductivity 

■49b  ■ 

(am/cm  ) 

(cal/gm-*C) 

(cm  / sec) 

(cal/sec-cm-°C) 

200 

1.665 

0.  184 

0.88  x 

10-3 

0.270  x  10~3 

250 

1.616 

0.  217 

1.02  x 

10'3 

0.358  x  10‘3 

300 

1. 570 

0.254 

1.20  x 

10~3 

0.478  x  10~3 

350 

1.525 

0.298 

1.46  x 

10'3 

0.663  x  10*3 

400 

1.483 

0.350 

i.80  x 

Iff3 

0.934  x  10"3 

450 

1.445 

0.410 

2.  19  x 

10'3 

1.30  x  10"3 

500 

1.412 

0.481 

2.61  x 

10-3 

1.77  x  10~3 

d.  Vapor  Pressure 

The  vapor  pressure  of  P^S^  was  measured  with  the  configuration  depicted  in 
figure  44.  As  shown,  argon  gas  was  used  to  balance  the  P4S3  vapor  pressure 
and  the  U-tube  extension  was  connected  to  a  mercury  manometer.  Measurement 
of  the  argon  pressure  to  maintain  equal  liqu’d  levels  in  both  sides  of  the  U-tube 
provided  P4S3  vapor  pressure  values  up  to  about  940°  F. 


Figure  44.  U-T.ibe  for  Vapor  Pressure  Measurement 


AFWL-TR  -67-114 


The  experimental  data  are  shown  in  figure  45.  Regression  analysis  yielded 
the  following  equation: 


P  =  6.86  - 


2.70  x  10‘ 


P  -  mm  Hg 
T  =  °K 

e.  Liquid  Density 

The  liquid  density  of  ?4Sg  was  measured  by  determining  the  volume  of  a 
known  weight  of  P^S^.  The  volume  change  with  temperature  was  noted  by 
observing  the  liquid  height  in  a  highly  uniform  and  accurately  measured  Vycor 
cylinder.  A  small  correction  for  the  Vycor  thermal  expansion  was  included. 

The  experimental  results  are  shown  in  figure  46  together  with  a  curve  which 
gives  the  best  representation  of  the  recently  obtained  and  previous  data. 

f.  Heat  of  Fusion 

The  heat  content  of  solid  and  liquid  P^S^  has  been  reinvestigated  about  this 
melting  point  (172°  C)  to  obtain  accurate  evaluation  of  the  heat  of  fusion.  The 
new  data,  together  with  that  obtamed  previously  (reference  12)  is  presented  in 
figure  47.  From  where  the  solid  and  liquid  heat  content  intercepts  the  172°  C 
temperature  line  it  can  be  seen  that  the  heat  of  fusion  is  nearly  10  cal/gm. 

The  experiments  were  performed  with  an  adiabatic  ice  calorimeter  (figure  48) 
and  a  vertical  drop -tube  furnace.  The  calorimeter  had  been  calibrated  to  +  3 
percent  accuracy. 

The  solid  P^S^  heat  content  data  between  50°  and  170°  C  can  be  expressed  as 

.  „  ,  ,  ,  v  6.  44  x  10‘3  T  ■‘•2.72 

AH  (cal/gm)  =  e 

where 


Differentiating  this  equation  with  respect  to  temperature  gives  the  solid  material 
heat  capacity  versus  temperature 


C  (cal/gm-°C) 
P  B 


*  [e.44  x  10'3J  [c 


6.44  x  10  ‘  T  +  2. 


1 


Vapor  Pressure  (mm  Hg) 


nsity  (gm/cm 
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Figure  46.  Liquid  Density  of  P.S 
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Similarly,  the  liquid  heat  content  follows  the  relationship 

A„.  .  3.  18  x  10*3T  +  3.425 

AH  (cai/gm)  *  e 

and  the  specific  heat  is  given  by 

_  ,  in-3  r  3.  18  x  10'3T  +  3.  4251 

Cp  (cal/gm-'C)  *  3. 18  x  10  Le  J 

between  180*  and  500°  C 
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SECTION  V 

IRRADIATION  OF  THEON  (P^) 

1 .  Purpose 

The  purpose  of  this  task  was  to  investigate  the  behavior  of  P^S^  in  an  operat¬ 
ing  reactor  environment  of  neutron  and  gamma  fluxes.  A  concurrent  analytic  - 
experimental  program  was  performed.  The  analytic  effort  enabled  study  of  the 
energy  exchange  mechanisms  and  thereby  yielded  information  for  predicting 
radiation  effects.  Actual  performance  data  were  obtained  by  irradiating  P^S^ 
in  the  Air  Force  NETF  reactor  at  cold  (<  200°  F)  and  elevated  temperatures 
approximating  TURPS  operation. 

The  major  efforts  under  this  task  included: 

(1)  Analytic  study  of  radiation  effects 

(2)  Design  and  fabrication  of  the  test  apparatus 

(3)  Post -irradiation  examination 

2.  Analytic  Study  of  Radiation  Effects 

The  radiation  damage  of  P^S^  results  from  the  neutron  and  gamma  fields 
in  the  TURPS  core.  To  estimate  the  magnitude  of  this  effect,  it  is  necessary 
to  calculate  the  energy  deposited  in  P^S^  by  neutrons  and  gammas  and  the 
fraction  of  this  energy  which  results  in  decomposition.  At  this  point,  it  should 
be  recalled  that  from  the  data  on  heat  of  formation,  it  is  evident  that  to  decompose 
a  P^Sg  molecule  ( i.  e .  ,  to  form  P  „  +  S^)  at  least  1.5  ev  of  energy  must  be  de¬ 
posited  for  each  P^S^  molecule. 

Gamma  radiation  interacts  mostly  with  the  bound  electrons,  producing  free 
electrons,  ions,  and  X-rays.  The  X-rays  lose  energy  by  collisions  with  electrons, 
converting  all  their  energy  to  heat.  The  ions  and  electrons  can  recombine, 
liberating  more  X-rays  (ultimately,  heat),  or,  if  the  lifetimes  of  the  ions  are 
large,  they  may  induce  or  accelerate  chemical  reactions  leading  to  molecular 
decomposition  products.  In  water,  ions  persist  for  a  long  enough  time  that 
ionic  reactions  constitute  an  important  source  of  stable  molecular  decomposition 
products.  However,  ions  are  not  nearly  so  stable  (have  shorter  lifetimes)  in 
PS-  as  m  water,  so  that  ion- ion  reactions  will  contribute  less  significantly  to 

*1  O 

radiation  decomposition  processes  which  terminate  : n  molecular  residue. 
Therefore,  it  will  be*  assumed  that  gamma  radiation  will  no'  contribute  signifi¬ 
cantly  to  the  decomposition  of  P^S., . 
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Neutron  damage  is  caused  primarily  by  elastic  collision  of  incident  neutrons 
with  atomic  nuclei  in  P4Sg.  In  a  fast  neutron  collision,  the  chemical  binding 
energy  of  the  scattering  atom  can  be  neglected.  For  example,  in  an  elastic 
collision  with  sulfur,  a  neutron  will,  on  the  average,  lose  six  percent  of  its 
energy.  This  means  that  a  1  Mev  neutron,  upon  colliding  with  sulfur,  will 
yield  a  940-kev  neutron  and  a  60-kev  sulfur  atom.  Only  when  the  incident 
neutron  energy  is  as  low  as  10  ev  does  the  chemical  binding  energy  of  the  atom 
become  important. 

To  determine  the  nature  of  neutron  damage,  it  is  necessary  to  follow  the 
recoil  and  determine  the  manner  in  which  it  loses  its  energy.  If  ionization  is 
the  mode  of  energy  loss  most  of  the  energy  will  be  converted  to  heat  as  is  the 
case  with  gamma  rays  and  X-rays.  On  the  other  hand,  elastic  collisions  will 
lead  to  continued  decomposition.  The  velocity  of  particles  below  which  ioniza¬ 
tion  becomes  negligible  is  given  by  (reference  13; 

V  *  jj-  *  *  0.22  x  109  cm/sec 

For  sulfur  this  velocity  corresponds  to  an  energy  of  800  Kev,  which  is  greater 
than  the  recoil  energy  resulting  from  an  elastic  collision  with  a  fast  neutron. 
Ionisation  is  not,  therefore,  the  mode  of  energy  loss  of  the  scattered  particle. 
Instead,  it  will  lose  energy  by  elastic  collisions,  thus  causing  further  decompo 
sition.  It  can  be  assumed  that  in  a  primary  collision  the  scattered  particle  will 
lose  all  its  energy  by  breaking  chemical  bonds. 

The  equilibrium  concentration  of  any  chemical  decomposition  product  will 
depend  on  its  production  rate  and  its  chemical  recombination  rate.  It  is  assumed 
that  all  chemical  decomposit  on  follows  the  reaction 


P4S3  radiation  +  3S 

It  is  also  assumed  that  the  back  reaction  rate  is  g^ven  by  the  following  second - 
order  rate  equation: 

-B[p43tg 


where 


B-B0e 


-E/kT 


and  [S^  3  4re  the  concentration*  of  and  5^  ,  re»pectivel 


and 


B  and  E  are  constants, 
o 
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The  temperature  dependence  of  B  is  that  of  a  characteristic  Arrhenius  reaction, 
where  the  chemical  reaction  is  exponentially  dependent  on  the  energy  of  the 
potential  barrier.  This  equation  holds  for  most  reactions  where  no  inter¬ 
mediate  products  are  present,  and  is  assumed  to  hold  for  the  present  case. 

The  net  P^  formation  rate  is  then  given  by 

d  [P,l 

-3t  *  A  -  B  f  J  [SJ  "fi> 

where  A  is  the  production  rate  of  1'^.  The  solution  to  equation  (18)  is 

M  ■  £ 

where  _ 

a=^AB 

if[P4]  =  [sO- 

The  corresponding  transient  time  (time  for  to  achieve  lie  of  its  equili- 
brium  value)  is  given  by  the  relation 


2  at 

e  -  1 


2  at  ,  , 
e  +  1 


(i'») 


T  = 


e  +  1 
e  -  1 


(20) 


and  the  equilibrium  concentration  is 


[pJ  ■£ 

o 


Equation  (19)  is  the  solution  tc  equation  (18)  when  all  the  free  sulfur  and  phcs- 
pnorus  formed  results  from  the  decomposition  . P^S^.  If,  however,  an  excess 
of  sulfur  is  present  and  if  the  concentration  o,  suit  r  does  not  change  signifi¬ 
cantly  during  the  decomposition  of  P^S^,  another  solution  of  equation  {  18)  results 


-  e'  B  tS3 

I  \ 


(22) 


The  transient  time  is 


B  f  S ' 


and  the  equilibrium  concentration  is 


A 


o : 


(21) 


AFWL-TR-67-1 14 


To  evaluate  the  decomposition  of  P^S^  in  any  reactor  environment,  it  is  nec¬ 
essary  to  evaluate  A,  and  E. 

The  production  rate.  A,  can  be  estimated  by  calculating  the  energy  deposited 
in  the  coolant  by  neutron  slowing  down  scattering. 


00 

$°SN 


2  M  3 

- dE  ev/ sec  cm 

(M  +  1) 


2  M 

(M  +  l)2 


(25) 


(20) 


in  the  rotation  of  the  multigroup  specification  of  energy  space. 

N  -  atom  density  in  P.S0  coolant 

4  3 

a  -  scattering  cross  section 
& 

$  =  neu-ron  flex 

M  =  atomic  mass  of  scattering  nucleus  (30.  9  for  P,  32.1  for  S) 

E.  =  is  the  midDOint  energy  of  group  i 

Ei 

♦j  •  \  <l>  (E)  dE 

Ei-! 

For  the  TURPS  reactor 

an  a 

Q  =  2.58  x  10fiU  evisec-cm 

Since  three  sulfur  atoms  must  be  removed  from  P4Sg  to  form  P^,  the  energy  of 
decomposition  is  1 . 5  ev,  and  this  yit  ’  .  a  producti'  n  rate  of 


A  =  1.72  x  1016/cm2-sec 

3  is  called  the  frequency  factor  in  the  Arrhenius  rate  equation.  It  is  related 

O  ' 

to  the  degrees  of  freedom  of  reacting  molecules  and  to  the  temperature.  A  con¬ 


servative  estimate  of  Bq  is  obtained  from  two  comple).  molecule",  reacting  which 

yields 


;,n-15  2 

-  Id  cm  /moiccule-sec 


(reference  14) 


o 


The  activation  energy  must  be  determined  from  experiment  and  is  at  the 
present  unknown.  I*  is  estimated  that  10  keal/moie  is  a  highly  conservative 
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value  in  view  of  the  activation  energy  of  other  reactions.  Consider,  for  instance, 
the  activation  energy  for  these  typical  secord-order  Arrhenius  reactions  (data 
from  reference  14): 

H  +  H2 - ^  H9  +  H  E  =  7.  5  kcal 

Cl  +  H2 - ^HC1  +  H  E*  5.5  kcal 

CH„  +  C„Hr - -*C,HQ  E  =  8.0  kcal 

3  o  o  4  9 

Finally,  computation  of  the  back-reaction  constant  obtains 

„  ir.-15  -  10  kcal/kt 

B  -  10  e 

“18  3 

B  =  10  cm  /molecule -sec 

These  values  yield,  when  substituted  into  equations  (20)  and  (21/ 
t  =  5.0  seconds 

17  3 

[PJ  =2.3x10  molecules /cm 

o 

The  concentration  of  P.S0  in  the  coolant  is 

4  3 

21  3 

[P4S3J  -  1.30  x  10  molecules /cm 

so  that  represents  0.017  percent  of  the  coolant  molecules. 

If  an  excess  of  one  weight  percent  sulfur  is  assumed,  equations  (33)  and  (24) 
yield 

T  =  1.1  x  10  ^  seconds 

[P^]  =  1.92  x  1014  molecules/cm3 

'  o 

so  that  an  excess  of  sulfur  reduces  the  phosphorus  concentration  considerably . 

The  phosphorus  concentrations  derived  are  those  which  would  exist  in 
solution  with  P4Sg.  They  should  not  be  taken  as  the  concentrations  which  exist 
in  the  vapor  i.°ar  the  condensing  surfaces.  To  determine  the  amount  of  free 
phospnorus  which  accumulates  in  the  vapor  is  a  problem  yet  to  be  investigated, 
it  will  depend  upon  the  thermodynamic  and  chemical  properties  of  both  phases, 
and  it  may  exceed  the  concentration  in  liquid  phase. 
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3.  Design  and  Fabrication  of  Test  Apparatus 

a.  Low  Temperature  Apparatus 

Investigation  oi  P4Sg  behavior  in  a  reactor  environment  of  neutron  and 
gamma  fluxes  was  pursued  by  irradiation  tests  performed  in  the  Air  Force 
NETF  reactor.  Capsules  containing  purified  stoichicmetric  P  .Sq  were  subjected 

19  q  •* 

to  an  integrated  neutron  dose  of  about  10  nvt.  The  samples  were  irradiated  at 
ambient  temperature  (gamma  heating  to  about  200°  F)  and  at  1000°  F.  Glass  en¬ 
capsulation  served  to  contain  P^Sg  with  and  without  the  presence  of  Zircaloy 
specimens  to  test  the  effect  on  metal  reaction  on  the  radiation  stability  of  P^S,, . 

During  the  first  irradiation  in  which  the  test  apparatus  was  under  water  the 
heated  capsules  developed  thermocouple  performance  problems  and  were  not 
able  to  complete  the  required  300  hours  irradiation.  The  unheated  capsules 
performed  satisfactorily.  Subsequent  examination  of  the  first  set  of  heated 
capsules  revealed  water  leakage  and  destruction  of  the  specimens. 

The  results  of  the  unheated  capsule  irradiation  are  presented  later  in  this 
report.  The  unheated  capsule  test  apparatus  consisted  of  the  glass  vial  as 
shown  in  figure  49  further  protected  by  an  aluminum  outer  container. 

b.  High  Temperature  Apparatus 

The  desire  to  obtain  higher  fast  neutron  fluxes  than  are  available  in  the 
NETF  Bulk  Shield  Facility,  together  with  the  obvious  advantage  of  avoiding  the 
water  environment  which  led  to  the  failure  of  the  first  tests,  prompted  use  of 
the  NETF  West  Test  Cell  for  the  second  high  temperature  P^Sg  irradiation. 

The  air  environment  of  the  test  cell  enabled  use  of  a  conventional  furnace  ar¬ 
rangement  with  doubly  encapsulated  P^Sg  specimens  located  inside  *he  heated 
zone.  Further,  a  large  part  of  the  heating  was  supplied  by  core-evolved  gammas 
and  only  supplemental  electrical  heating  by  the  furnace  was  required. 

Figure  50  illustrates  the  furnace  and  capsule  configuration.  Six  P^Sg  cap¬ 
sules  were  contained  inside  a  Hoskins  furnace  core.  Because  of  the  large 
gamma  heating  contribution  and  relatively  low  heat  loss  to  the  surrounding  air, 
only  minimal  insulation  was  required  at  the  top  and  bottom.  Figures  51  and  52 
are  photographs  of  the  completed  furnace  assembly  and  its  installation  op  the 
NETF  West  Test  Cell  ooUy. 

Each  of  the  six  P^Sg-containing  quartz  capsules  was  doubly  encapsulated  in 
titanium,  a  material  which  incurs  low  activation  and  has  a  short  half-life.  The 
P^Og  specimens  were  composed  of  stoichiometric  P^g,  P^Sg  Wlth  °-5  percent 
excess  sulfur,  and  P4Sg  with  1  percent  excess  sulfur.  Three  vials  contained 
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(1)  Aluminum  plate  10  inch  OD  x  1/16 -inch  thick  (2,  top  and  bottom) 

(2)  Fir*  brick  6  inch  OD  x  1  / 4 -inch  thick  (2,  top  and  bottom) 

(3)  Furnace  core  (SIC  wound  with  nichrome  wire) 

(4)  Chromel-alumel  thermocouple*  with  ceramic  insulator* 

(5)  Outer  titanium  capsules  (-5- 1/2  inches  high  x  1  inch  OD  x  0.012  inch  wall 
Six  tubes  placed  as  shown  on  3 -inch  ♦  circles  60  degrees  apart 

(6)  Inner  titanium  capsules  (-5  Inches  high  x  3 /4-inch  OD  x  0. 040  inch  wall 

(?)  tubes  (6)  containing  10  grams  P^Sg  each.  15-millimeter  diameter 

x  1-millimeter  wall,  -4  inches  in  height 

(8)  Aluminum  holddown  rods  (4)  spaced  90  degrees  apart 

(9)  Glass  wool 


l  , 


Figure  50.  Furnace  Configuration  o'i  High  Temperature  P^S3  Irradiation 
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Figure  52.  Installation  of  High  Temperature  P^Sg  Irradiation 
Teat  Configuration  on  NETF  Weat  Teat  Cell  Dolly 
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these  compositions  only;  and  the  three  remaining  vials  contained  a  presulfided 
Zircaloy-2  coupon  in  addition.  The  excess  sulfur  compositions  were  studied 
since  analytic  predictions  show  greatly  reduced  dissociation  for  such  material. 
The  reasons  for  presulfiding  of  the  Zircaloy-2  are  given  in  Section  VI  of  this 
report. 

Figure  49  shows  the  quartz  vial  and  double  titanium  containment  sleeves 
prior  to  assembly.  Each  quartz  vial  was  approximately  four  inches  long  and 
contained  approximately  10  grams  P4Sg. 

Because  of  uncertainties  in  the  West  Test  Cell  gamma  flux,  the  actual 
operating  reactor  power  level  during  the  test,  and  the  heat  transfer  calculations, 
a  duplicate  furnace  with  mockup  test  specimens  was  installed  and  operated  until 
equilibrium  temperatures  were  achieved.  After  verification  that  the  configura¬ 
tion  could  maintain  internal  capsule  temperatures  at  about  950°  F,  the  mockup 
was  discarded  and  the  actual  P4S3  test  assembly  installed.  Subsequent  opera¬ 
tion  of  the  assembly  proceeded  smoothly  for  the  desired  total  of  300  hours  with 
approximately  250  watts  supplemental  electrical  heating  required  to  maintain  the 
capsule  walls  at  900*  F.  The  corresponding  internal  temperature  was 
about  950*  F. 

4  Post-Irradiation  Examination 

a.  Low  Temperature  Specimens 

The  post -irradiation  examination  of  the  low  temperature  P4Sg  included 
melting  point,  viscosity  and  visual  examination  as  well  as  inspection  of  the 
Zircaloy-2  specimens  contained  in  two  of  the  capsules.  Gold  foil  flux  monitors 
(some  cadmium  covered)  enabled  evaluation  of  the  thermal  and  epicadmium 
neutron  fluxes  for  each  capsule  tested.  Table  XIX  presents  the  flux  data. 


Table  XIX 

NEUTRON  FLUX  EXPOSURE  OF  THE  LOW 

TEMPERATURE  P.S,  CAPSULES 
4  3 


Capsule 

number 

Neutron  flux 

2 

(n/cm  -sec) 

Time  integrated  flux 
(300-hr  nvt)(n/cm2) 

Thermal 

i  Epicadmium 

Thermal 

I  Epicadmium 

2 

0.35  x  1012 

-0. 15  x  1011 

0.4  x  1018 

-0.2  x  1017 

3 

0,94  x  1012 

-0.5  x  1011 

1.0  x  1018 

-0.5  x  1017 

6 

1.2  x 1012 

-0.6  x  1011 

1.3  x  1018 

17 

-0.6  x  10 

9 

2.8  x  1012 

-1.2  x  1011 

2.8  x  1018 

17 

-1.3  x  10  ' 

(1)  Melting  point  evaluation 

Measurement  of  the  normal  melting  point  of  the  irradiated  P4S3  both  as 
irradiated  and  after  thermal  cycling  showed  essentially  insignificant  changes 
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from  unirradiated  P^Sg.  Table  XX  lists  the  values  for  each  specimen,  the  ex¬ 
perimental  accuracy  being  +  2.0°  C.  It  is  seen  that  all  measurements  are  quite 
close  to  the  accepted  melting  point  of  172°  C. 


Table  XX 

MELTING  POINTS  OF  IRRADIATED  P4Sg 


Capsule  number 

Melting  point 
as-irradiated 

TO 

After  thermal 
cycling 

2 

174 

172.5 

3 

173 

170.5 

6 

172.5 

171. 5 

9 

172.5 

172.5 

(2)  Viscosity  determination 

The  evaluation  of  liquid  viscosity  dependence  on  temperature  for  the  irradiated 
P^Sg  was  made  in  a  similar  manner  to  that  described  for  the  unirradiated  deter¬ 
minations  presented  previously.  Table  XXI  contains  the  measurements  for  the 
temperature  range  of  600°  to  950°  F. 

Table  XXI 


VISCOSITY  MEASUREMENTS  OF  IRRADIATED  P4Sg 


Specimen  number 
viscosity  (centipoise) 


(°F> 

1  2 

!  3 

6 

9 

600 

1.177 

620 

1.072 

0.879 

640 

-- 

-- 

1.131 

680 

-- 

0.796 

-- 

720 

-- 

-- 

0.612 

840 

-- 

0  525 

-- 

890 

0.378 

... 

920 

-- 

--- 

0.513 

950 

0.324 

0.412 

Similarity  between  the  viscosities  of  irradiated  and  umrradiated  P.S.,  is 

*  V/ 

shown  in  figure  53  which  contains  the  values  of  table  XXI  plotted  over  the  best 
fit  line  from  the  unirradiated  data  (from  figure  39).  It  is  apparent  that  no 
significant  viscosity  property  change  has  occurred. 
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Irradiated  P.S„  Data 
4  3 

0  Capsule  2 
V  Capsule  3 
□  Capsule  6 
A  Capsule  9 


.Unirradiated 


i/T  x  iO3  rB'1) 

- - - 1 - 1_ _ _ _ 

700  80S 

Temperature  <*F! 

Figure  53.  Viscosity  e?  Irradiated 
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(3)  Visual  inspection  of  P^S^  and  Zircaloy-2  coupons 

Visual  examination  of  the  P^S^  subsequent  to  irradiation  and  upon  melting 
during  the  melting  point  and  viscosity  tests  revealed  normal  P4S3.  The  initial 
coloring  and  darkening  upon  temperature  raising  was  identical  to  that  of  unir¬ 
radiated  material. 

Capsules  2  and  3  contained  coupons  of  Zircaloy-2  about  0.5  by  0.25  by  0.030 

inch.  As  expected,  these  specimens  exhibited  no  change  in  appearance  since 

the  surrounding  P^S3  was  at  low  temperature  and  not  molten.  Also,  the  P^S3 

in  these  capsules  exhibited  the  same  properties  as  the  material  in  capsules  6 

and  9, which  contained  P.S„  alone. 

4  3 

f4)  Analysis  and  summary 

The  fast  neutron  flux  exposure  of  the  irradiation  specimens  was  about  0  003 
of  the  TURPS  fast  flux.  This  yields  an  analytical  prediction  of  the  equilibrium 
concentration  of  about  0.001  percent  of  the  Theon  molecules  at  950°  F. 

Such  a  small  composition  change  is  in  agreement  with  the  observed  absence  of 
discernible  radiation  effects. 

As  previously  shown,  the  analytic  model  indicates  only  0,02  percent  [pJ 
concentration  at  TURPS  operating  conditions.  Such  a  small  dissociation  should 
yield  essentially  no  coolant  property  changes. 

b.  High  Temperature  Specimens 

Post-irradiation  examination  of  each  of  the  six  Theon  specimens  irradiated 

at  approximately  9503  F  included  melting  point,  viscosity  and  visual  studies  as 

well  as  inspection  of' the  Zircaloy-2  specimens  contained  in  three  of  the  capsules. 

Gold  foil  flux  monitors  enabled  determination  of  the  average  thermal  neutron  flux 

exposure,  while  the  titanium  capsule  material  activation  from  the  Ti46  (rj,  p) 

46 

Sc  reaction  provided  a  measure  of  the  fast  neutron  flux  abov-e  2  .vlev.  Since  the 

test  w  as  conducted  in  the  nonmoderating  medium  of  the  test  cell  (air  and  heavy 

element  structure),  the  fasi  *rd  thermal  neutron  flux  distributions  were  quite 

12  2 

uniform  in  the  experimental  configuration.  Average  values  of  2.2  x  10  17/cm  - 

13  2 

sec  for  the  fast  flux  02  Mev)  and  1.0  x  10*  n/cm  -sec  for  the  thermal  flux 
were  obtained  with  only  about  10  percent  deviation  across  the  test  assembly. 

Table  XXII  contains  the  capsule  identification  and  contents  data. 
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Table  XXII 

HIGH  TEMPERATURE  THEON  IRRADIATION  SPECIMENS 


Capsule  No. 

Theon  Composition  \ 

Zr-2  Coupon 

1 

Stoichiometric  P^S^ 

No 

2 

Stoichiometric  P.S0 

4  3 

Yes 

3 

P4S3  +0.5  wt%S 

No 

4 

P  .S0  +  0.  5  wt  %  S 

Yes 

5 

P4S3  +  1.0  wt  %  S 

No 

6 

P,S0  +  1.0  wt  %S 

4  3 

Yes 

(1)  Visual  inspection  of  Theon  and  Zircaloy-2  coupons 

Visual  examination  of  the  Theon  subsequent  to  irradiation  showed  the  pres¬ 
ence  of  dark  deposits  on  the  walls  of  the  quartz  containers.  Most  of  this  mate¬ 
rial  remained  attached  to  the  walls  upon  transfer  of  the  Theon  to  the  viscosity 
test  apparatus.  However,  specimens  3  and  4  contained  sufficient  insoluble 
material  within  the  Theon  matrix  that  viscosity  measurements  were  prevented 
because  of  clogging  of  the  small  venturi  tube  in  the  test  apparatus.  This  mate¬ 
rial  is  believed  to  be  a  reaction  product  with  the  glass  co  tainer  surface. 

Despite  this,  however,  all  of  the  Theon  specimens  exhibited  the  normal  color 
change  from  bright  yellow  to  dark  brownish  red  upon  melting  and  increasing  of 
the  temperature  to  950°  F. 

The  Zr-2  coupons  exhibited  a  uniform  dark  coating  that  was  well  adhered  to 
the  surface.  The  presence  of  zirconium  sulfide  was  verified  by  the  X-ray  dif¬ 
fraction  study. 

(2)  Viscosity  determination 

The  evaluation  of  liquid  Theon  viscosity  dependence  on  temperature  for  the 
high  temperature  irradiated  specimens  was  made  in  a  similar  manner  to  that 
described  for  toe  unirradiated  measurements  presented  previously.  Table  XXIII 
contains  the  results  of  these  measurements.  An  unirradiated  specimen  (7)  was 
included  to  provide  a  check  of  the  method  and  to  obtain  values  for  Theon  contain¬ 
ing  with  one  weight  percent  exces'  sulfur. 
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Table  XXIII 

VISCOSITY  MEASUREMENTS  OF  THEON 
IRRADIATED  \T  950°  F 


Temperature 

t*F) 

Specimen  Number 

Viscosity  (centipoise) 

1 

2 

5 

6 

i 

610 

1 .  149 

625 

1 . 255 

630 

1.034 

660 

1 . 283 

680 

0.987 

750 

0.  656 

765 

0.605 

7  85 

0.514 

790 

0.  572 

905 

0.448 

0.396 

925 

0.450 

930 

0.451 

935 

0.  420 

Similarity  between  the  viscosities  of  irradiated  and  unirradiated  is 

apparent  in  figure  54  which  contains  the  values  of  table  XXIII  plotted  over 
the  best  fit  line  from  the  umrradiated  data.  No  significant  viscosity  change 
is  observed.  In  addition,  the  umrradiated  with  one  weight  percent 

excess  sulfur  exhibits  the  «arne  viscosity  as  stoichiometric  P.S,. 

(3)  Melting  point  evaluation 

Measurements  of  the  normal  melting  point  of  Theon  irradiated  w  ithout 
the  presence  of  Zircaiov-2  she  .ed  essentially  insignificant  changes  from 
umrradiated  Theon.  However,  the  melting  points  obtained  from  samples 
irradiated  w  ith  Zircaloy-2  coupons  were  2*  to  3°  F  lower.  This  Entering  of 
the  melting  point  may  be  attributed  to  the  presence  of  dissolted  metal  sulfide 
in  the  Th  con.  This  observation  agrees  with  indications  of  limited  solubility 
of  the  metal  surface  sulfide  layer  obtained  daring  s  o  t'jce  treatment  ones'  - 
gat  ions  presented  in  Section  VI  of  this  report. 
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Theon  irradiated 
at  950°  F 

*  .  O  Specimen  1 

A  Specimen  2 
□  Specimen  5 
V  Specimen  6 
O  Specimen  7 


0.  101 _ I _ I _ i _ I _ I _ I _ I 

0.95  0.90  0.85  0.80  0.75  0.70 

1/T  x  103  ( °R " 1 ) 

I _ I _ I _ I _ I 

600  700  800  900  1000 

Temperature  (°F) 

Figure  54.  Viscosity  of  Theon  Irradiated  at  950°  F 
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(4)  Analysis  and  summary 

The  fast  neutron  flux(>2  Mev)  exposure  of  the  high  temperature  irradiation 
test  was  about  0.  1  of  the  TURPS  fast  flux  (about  30  times  the  fast  flux  received 
by  the  low  temperature  test).  This  yields  an  analytical  prediction  of  the 
equilibrium  [VJ  concentration  for  stoichiometric  of  about  0.002  percent 

of  the  number  of  molecules  at  950°  F.  Such  a  small  composition  change  is  in 
agreement  with  the  observed  absence  of  coolant  property  changes.  Further¬ 
more,  as  shown  in  the  analytic  study,  the  addition  of  excess  sulfur  greatly 
reduces  the  equilibrium  JpJ  concentration,  but  since  all  coolant  property 
changes  were  negligible,  distinction  of  the  sample  behavior  with  excess  sulfur 
present  was  not  possible. 

The  excess  sulfur  specimens  did,  however,  contribute  an  interesting  obser¬ 
vation  for  the  capsules  containing  a  Zircaloy-2  coupon.  Small  amounts  of  free 
phosphorus  were  detected  in  the  specimens  containing  stoichiometric  P^Sg  and 
0.5  weight  percent  excess  sulfur  in  F^Sg,  but  no  phosphorus  was  observed  in 
the  one  weight  percent  excess  sulfur  specimen.  Thus,  it  appears  that  reaction 
of  the  Theon  with  the  Zircaloy-2  surfaces  yielded  some  free  phosphorus  when 
less  than  0.  5  weight  percent  excess  sulfur  was  present,  but  the  one  weight 
percent  excess  sulfur  condition  prevented  the  freeing  of  phosphorus.  Quanti¬ 
tative  evaluation  of  the  rate  and  extent  of  Theon  reactions  with  Zircaloy  as  a 
function  of  Theon  composition  (percent  excess  sulfur)  will  be  the  subject  of 
future  investigations. 

Both  the  low  and  high  temperature  irradiation  tests  have  supported  the 
analytic  prediction  of  a  very  high  recombination  of  P^S^  dissociation  products 
in  a  radiation  field.  Insignificant  property  changes  were  noted  with  and  without 
the  (.'^esence  of  Zircaloy-2,  a  promising  Theon  containment  material  for  TURPS. 
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SECTION  VI 

CHEMICAL  COMPATIBILITY  OF  THEON  WITH  METALS 

1.  Purpose 

Previous  investigation  of  purified  compatibility  with  Zircaloy  indicated 

no  evidence  of  corrosive  attack.  However,  earlier  studies  of  metals  with  com¬ 
mercial  grade  P^Sg  showed  small  rates  of  corrosion  on  Zircaloy  and  a  much 
larger  attack  of  austenitic  stainless  steels  such  as  SS  304.  The  purpose  of  the 
present  task  was  to  evaluate  the  performance  of  structural  materials  through 
static  corrosion  tests  of  SS  304,  SS  310,  Hastelloy-N,  Zircaloy-4  and  Zircaloy- 
2  with  purified  P^Sg.  In  addition,  the  effects  of  excess  sulfur,  excess  phos¬ 
phorus  and  a  small  amount  of  chlorine  added  to  the  P^Sg  were  to  be  investi¬ 
gated. 

A  second  major  objective  of  the  compatibility  studies  was  the  performance 
of  a  dynamic  corrosion  loop  test  of  4500 -hour  duration  at  simulated  TURPS 
conditions  incorporating  the  material  selected  from  the  static  screening  tests. 

2.  Static  Corrosion  Tests 

a.  Experiment  Description 

Five  materials  were  evaluated  during  this  investigation  based  on  anticipated 
corrosion  resistance  in  the  P4Sg  environments.  Specimens  of  Zircaloy-4, 
Zircaloy- 2,  Hastelloy-N  and  Types  310  and  304  stainless  steels  were  selected 
for  exposure  to  Theon  environments  consisting  of  P4Sg  (stoichiometric), 

P4S3  +  5%S,  P4Sg  +  1%S  and  P^Sg  +  5%P.  In  addition,  specimens  were  tested 
in  stoichiometric  P4Sg  and  P4Sg  +  5%S  with  a  small  amount  of  chlorine  gas 
(5  ppm)  to  duplicate  conditions  expected  after  five  years  of  reactor  operation. 
Each  metal  was  tested  in  a  capsule  constructed  of  the  same  material  tc  eliminate 
foreign  corrosion  products.  Samples  were  suspended  from  brackets  and  im¬ 
mersed  in  liquid  P4Sg,  in  the  vapor  phase  and  approximately  at  the  liquid-vapor 
interface. 

All  capsule  loading  and  final  welding  was  completed  in  an  inert  (argon) 
atmosphere.  Each  capsule  was  evacuated  to  0.  025  mm  Hg  vacuum  and  sealed. 
The  chlorine  gas-containing  capsules  were  evacuated  and  backfilled  from  a  mani¬ 
fold  with  measured  quantities  of  gas  (figure  55).  All  tests  were  run  at  1000°  F 
for  500  hours.  The  test  apparatus  is  shown  in  figure  56.  The  specimens  after 
exposure  are  shown  in  figures  57  through  63. 
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b.  Visual  Inspection  of  Test  Samples 

(1)  Zirealoy-4 

This  material  showed  excellent  performance  in  all  environments  tested. 
X-ray  diffraction  of  the  surface  after  test  indicated  the  presence  of  an  adherent 
zirconium  sulfide  layer.  The  samples  immersed  in  the  liquid  and  at  the  inter¬ 
face  showed  the  greatest  tendency  to  form  this  barrier. 

In  no  case  was  spalling,  pitting  or  cracking  of  the  material  noted.  Further, 
the  presence  of  5  ppm  chlorine  had  no  adverse  effects  on  the  material. 

(2)  Zircaloy-2 

When  tested  in  unchlorinated  P^S^,  this  material  displayed  the  same  cor¬ 
rosion  resistant  characteristics  as  Zircaloy-4.  Visual  examination  indicated 
that  the  material  was  susceptible  to  some  uniform  corrosion  prior  to  the  for¬ 
mation  of  the  zirconium  sulfide  barrier  which  was  present  on  all  samples. 
However,  the  samples  exposed  to  and  chlorine  developed  a  nonadherent 

zirconium  sulfide  layer,  the  cause  of  which  is  unknown  at  present. 

(3)  Hastelloy-N 

Although  the  gross  corrosion  rates  in  the  nonchlorine  bearing  environments 
are  not  large,  a  spalling  type  of  attack  was  noted  on  several  samples  (figure 
57 --Hastelloy-N  E-6;  figure  60- -Hastelloy-N  E-ll).  The  material  showed 
good  resistance  in  the  P4S^  +  1%S  media  and  in  both  P4S3 -chlorine  environ¬ 
ments.  Selective  pitting  and  cracking  or  intergranular  attack  were  not  noted 
either  visually  or  by  metallographic  examination. 

(4)  Type  310  stainless  steel 

The  specimens  tested  exhibited  significant  corrosion  in  all  environments 
with  the  greatest  attack  occurring  in  the  liquid  phase  samples.  The  mode  of 
attack  was  one  of  erosion  with  possibly  some  exfoliation  occurring  during  test. 
The  rate  of  attack  was  lowest  in  stoichiometric  P^S.^.  The  introduction  of 
five  percent  excess  sulfur  or  phosphorus  yielded  accelerated  attack  in  the 
liquid  phase. 

(5)  Type  304  stainless  steel 

This  material  displayed  similar  corrosion  properties  to  those  mentioned  for 
Type  310  material,  although  somewhat  greater  susceptibility  was  noted.  Again, 
the  samples  tested  in  stoichiometric  exhibited  the  least  corrosion.  When 


127 


AFWL-TR-67- 114 


tested  in  stoichiometric  +  Cl  and  P  ,S3  +  5%S  +  Ci.  the  samples  in  the 
liquid  phase  formed  a  thin,  nonadherent  scale.  Samples  of  this  scale  sub' 
mitted  for  analysis  by  X-ray  diffraction  indicate  the  presence  of  nickel  and 
iron  oxide.  Chlorides  were  not  detected. 

c.  Quantitative  Evaluation  of  Test  Results 

The  weight  changes  incurred  by  all  corrosion  specimens  are  presented  in 
tables  XXIV  to  XXIX.  Normalization  to  the  sample  area,  however,  is  not 
directly  applicable  for  the  Zircaloy  pieces  since  considerable  surface  roughness 
existed  prior  to  testing.  (It  is  estimated  that  about  a  30  to  40  percent  effective 
area  increase  was  present. )  For  comparison  purposes,  the  liquid  and  vapor 
phase  data  are  plotted  in  figure  64.  The  interface  sample  data  are  omitted  be¬ 
cause  of  uncertainties  in  liquid  level  position  for  the  various  capsules. 

(1)  Zircaloy 

The  values  for  liquid  phase  sulfide  film  buildup  on  Zircaloy  agree  well  with 

predictions  made  from  previous  data  (reference  12).  The  empirical  formula 

in  the  referenced  document,  w  =  aeD  ,  yields  about  40  mg/dm  for  a  500- 

2 

hour  test  at  1000°  F.  The  average  observed  value  of  nearly  60  mg/dm  is 

2 

actually  about  45  mg/dm  when  corrections  for  surface  roughness  are  made, 

In  accordance  with  previous  tests,,  the  vapor  phase  specimens  experienced 
extremely  small  weight  change. 

(2)  Stainless  steel 

The  steel  samples  experienced  relatively  large  weight  losses  when  exposed 

to  the  liquid  Theon  environment.  The  absence  of  an  adherent  protective  sulfide 

2 

layer  caused  changes  of  200  to  300  mg/dm  for  stoichiometric  or  excess  sulfur 
conditions.  The  presence  of  chlorine  appeared  to  accelerate  the  corrosion. 
However,  uncertainty  in  the  amount  of  chlorine  added  to  the  steel  capsules  pre¬ 
vents  quantitative  estimation  of  the  effect. 

It  is  interesting  to  note  that  the  vapor  phase  specimens  showed  very  little 
change.  Further  tests  are  necessary  to  assure  reproducibility  of  this  behavior, 
but  indications  of  applicability  to  the  vapor  region  of  the  TURPS  system  arise 
from  this  result. 
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(3)  Hastelljy 

The  Hastelloy  specimens  showed  lower  average  corrosion  rates  than  the 
steel  specimens,  but,  as  indicated  previously,  irregular  spalling  was  ob¬ 
served.  Hence,  despite  lower  observed  total  weight  change,  local  attack  may 
be  most  significant. 

3.  Dynamic  Corrosion  Test 

The  static  corrosion  test  established  Zircaloy  as  the  most  promising 
structural  material  for  containment.  The  primary  aim  of  the  present 

effort  w?  s  to  evaluate  the  performance  of  Zircaloy  in  a  long-term  test  with 
P4S3  conditions  approximating  TURPS  operation.  To  this  end.  the  design  and 
manufacture  of  a  dynamic  corrosion  test  apparatus  capable  of  4500  hours  oper¬ 
ation  was  completed.  To  aid  the  study  of  Zircaloy  performance,  an  additional 
short-term  investigation  was  conducted  to  explore  differences  of  behavior  be¬ 
tween  Zircaloys-2  and  -4  and  effects  of  pretreatment  of  the  Zircaloy  surfaces. 

a.  Surface  Treatment  Investigation  of  Zircaloy- 2  and  Zircaloy  -4 

The  effects  of  surface  preparation  on  Zircaloys-2  and  -4  were  studied  by 
exposing  variously  treated  test  coupons  to  P^  at  1000°  F  for  periods  of  U0 
and  540  hours.  All  specimens  were  degreased  in  Oakite  and  chemically 
cleaned  in  a  10  percent  HNO^-IO  percent  HF  solution.  In  addition,  air  and 
vacuum  firing  operations  were  performed  on  selected  samples.  Subsequent 
to  the  initial  preparation,  coupons  of  each  Zircaloy  and  treatment  type  were 
presulfided  for  24  hours  at  1000°  7  in  pure  sulfur  and  in  pure  P^Sg.  Following 
the  presulfiding,  immersion  of  the  specimens  in  P^Sg  at  1000°  F  was  conducted. 
In  order  to  minimize  container  reaction  effects,  the  exposure  was  conducted  in 
thoroughly  cleaned  Vycor  capsules. 

Study  of  the  test  coupons  and  P^S^,  after  the  desired  exposure,  was  made 
by  weight  change  determination,  metallographic  observation.  X-ray  diffraction 
and  emission  spectroscopy.  The  following  significant  conclusions  were  drawn 
from  this  investigation: 

(1)  No  difference  was  noted  between  the  behavior  of  Zircaloy  2  and 
Zircaloy -4. 

(2)  Pretreatment  of  the  Zircaloy  by  air  firing,  vacuum  firing  or 
chemically  cleaning  alone  did  no*  affect  thp  subsequent  compatibility  with 
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(3)  Presulfiding  of  Zircaloy  with  either  P4S3  or  sulfur  decreased  the 
subsequent  contamination  of  P^Sg.  The  degree  of  discoloration  of  P^S^  for 
the  170-  and  540-hour  tests  was  much  less  than  that  observed  in  the  initial 
presulfiding  in  P^S^. 

(4)  There  was  no  difference  in  the  behavior  between  samples  pre¬ 
sulfided  with  sulfur  or  P.S.,. 

4  3 

(5)  Indications  of  limited  solubility  of  the  presulfided  layer  during  sub¬ 
sequent  P4S3  exposure  were  observed.  This  conclusion  results  from  the  obser 
vation  that  the  170-hour  specimens  showed  small  weight  reductions  whereas 
the  540-hour  specimens  incurred  weight  increases. 

b.  Dynamic  Corrosion  Test  Apparatus  Design 

A  detailed  design  of  a  dynamic  corrosion  test  loop  was  prepared.  Photo¬ 
graphs  of  the  completed  loop  are  shown  in  figures  65  through  67.  Figure  65 
presents  the  bare  Zircaloy  structure  while  figure  66  illustrates  the  installation 
of  the  main  heaters,  guard  heaters  and  thermocouples.  Figure  67  contains  the 
complete  assembly  with  insulation  attached  prior  to  insertion  into  the  inert 
atmosphere  outer  container. 

To  assure  proper  fabrication  of  the  Zircaloy  loop,  extensive  radiographs 
of  welds  and  leak  check  tests  were  performed.  Figure  68  shows  the  radio¬ 
graph  of  the  high  velocity  test  section  and  two  thermocouple  wells.  It  should 
be  noted  that  the  one -half  of  the  high  velocity  section  length  was  constructed 
of  Zircaloy-2  and  the  other  half  of  Zircaloy-4.  This  enabled  further  evalua¬ 
tion  of  the  relative  performance  of  the  Zircaloys. 

c.  Dynamic  Corrosion  Test  Operation 

Because  the  surface  treatment  test  results  indicated  that  presulfiding  of 
Zircaloy  minimizes  the  reaction  with  P4S3,  the  test  apparatus  was  operated 
for  approximately  two  days  with  an  initial  loading  of  P4S3,  plus  one  weight 
percent  of  excess  sulfur.  The  initial  heatup  proceeded  smoothly  with  all 
heaters  and  thermocouples  functioning  satisfactorily.  Small  valve  leaks  in 
the  pressure  gage  and  fill  lines  were  remedied  by  replacement  with  steel 
bellows-type  valves  capable  of  withstanding  higher  temperatures  for  longer 
periods  of  time.  The  test  apparatus,  during  this  presulfiding  period,  vented 
with  ease  and  drained  in  a  clean  and  straightforward  manner.  No  attempts 
were  made  to  obtain  accurate  heat  balance  or  flow  determinations,  the  objec¬ 
tive  being  short-term  operation  at  950°  to  1000°  F  with  all  Zircaloy  surfaces 
in  contact  with  P^Sg  liquid  or  vapor.  Following  the  draining  of  the  apparatus. 
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it  was  recharged  with  fresh  stoichiometric  P4Sg,  When  venting  was  attempted, 
it  was  immediately  evident  that  the  vent  lines  were  clogged  with  condensed 
material.  Internal  pressure  measurements  indicated  the  presence  of  non¬ 
condensibles,  as  pressure  readings  were  higher  than  anticipated  from  the  known 
P4S3  vapor  pressure- temperature  relationship.  Inspection  of  the  vent  line 
disclosed  the  presence  of  red  phosphorus. 

Because  T*ed  phosphorus  (melting  point  500*  C)  could  not  be  melted  in  the 
vent  line,  physical  removal  by  coring  with  a  drill  was  attempted.  Although 
most  of  the  line  was  cleared  in  this  manner,  a  weld  penetration  caused  breakage 
of  the  drill  bit.  The  drill  bit  was  firmly  lodged  in  the  line,  and  attempts  at 
removal  were  unsuccessful.  It  was  concluded  that  removal  of  the  bit  necessi¬ 
tate:.  removing  the  outer  container  top  flange  and  cutting  and  re -welding  the 
vent  line.  In  addition,  disassembly  and  rewiring  of  thermocouple  and  heater 
circuits  was  required.  At  this  point,  operations  were  suspended  because  of 
funding  limitations. 

Failure  of  the  apparatus  operation  is  attributed  to  the  presence  of  free 
phosphorus  in  the  P^S^.  Initially,  the  free  phosphorus  appears  in  the  gaseous 
or  white  phosphorus  form,  but  at  the  conditions  of  this  test,  transformation  to 
red  phosphorus  occurs  rapidly.  The  red  phosphorus  is  liquid  at  the  operating 
test  temperature  (-550s  C)  and  would  aot  clog  the  passages  in  the  apparatus. 
However,  the  vent  line  heating  was  designed  only  for  melting  (plus  some 
margin,  i.  e. ,  5<XP  F)  and  the  vent  line  could  not  be  freed.  Hence,  removal 
of  noncondensibles  and  test  operation  could  not  proceed. 

d.  Conclusions  and  Recommendations 

The  decision  to  operate  fee  test  apparatus  with  stoichiometric  P4S3  appears 
to  have  been  most  unfortunate,  since  fee  indications  are  that  fee  presence  of 
one  percent  excess  sulfur  would  have  suppressed  fee  free  phosphorus.  There 
has  been  no  indication  from  prior  experimental  work  that  P4S3  plus  one  percent 
excess  sulfur  differs  in  corrosive  properties  from  the  stoichiometric 
material.  Moreover,  computations  indicate  that  the  presence  of  a  small 
quantity  of  excess  sulfur  substantially  increases  the  radiation  resistance 
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1  High  velocity 
teat  cectioi! 

2  Thermocouple  well 

3  Thermocouple  well 


Figure  68.  Radiographs  of  High  Velocity  Venturi  Teat  Section  and  Two  Tfcermocotmle  Walla 


Resumption  of  the  test  can  be  accomplished  without  major  changes  to  the 
apparatus.  Some  changes  in  the  fill  system  and  vent  line  configuration  should 
be  made  to  minimize  susceptibility  to  clogging  during  normal  venting;  venting, 
of  course,  being  required  to  blow  off  the  argon  blanket  used  during  the  fill 
process.  The  significant  finding  is  that  the  test  should  have  used  P4S3  with 
some  small  fraction  of  excess  sulfur.  Use  has  been  made  here  of  the  terms 
"corrosion"  and  "radiation  effects"  but,  in  fact,  the  work  to  date  indicates 
that  P4S3  has  negligible  corrosive  effect  on  presulfided  Zircaloy  and  that  the 
radiation  effect  on  P^Sg  plus  excess  sulfur  is  also  negligible. 
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APPENDIX  1 

APPROXIMATION  OF  AN  INFINITE  SERIES 

•  2  2 
V  e‘n  a  2 

Evaluation  of  ^  - 2 —  ^or  °  *  an  infinite  series  of  the 


n  *  1 


m 

form  ^  fn  (n)  can  be  approximated  by  a  finite  sum  plus  an  error  term, 
n  *  1 

Letting  flx)  be  the  continuous  function  which  agrees  with  f  (n)  at  the  integer, 
then  a  good  approximation  is  given  by 

f-1 


I  '  X  fn(n)  +  K(,,  +  $  t(x)dx 


n  *  1  n  *  1  t 

Derivation  of  this  formula  is  straightforward  from  figure  69.  Clearly,  the  sum 

•0  m 

of  the  series  is  bounded  between  ^  fix)  dx  and  ^  f  (x-1)  dx  and  as  A,  *  A, , 

1  1 


■  ■  a  • 

fn  ^  M  dx  *  J  f  (x-1)  dx  -  ^  fn 


T  *1 


A2  *  A2,  etc. ,  then 


(n) 


r.  =  1 


m  mm  a 

I  fn(n)+  Z  fn(n)*  J  f(x)dx+J 


(n)  + 
n  *  1  2 


f  (x)  dx+  ^  f  (x-1)  dx 
1  2 


(n)  -  fn(l)  *2  \  f  (x)  dx 


m 

$ 


n  *  1 


S‘  n 


f  (1) 


(n) 


m 

+  J  f  (x)  dx 
l 


n  *  1 


from  equating  Aj  to  Aj  and,  as  the  approximation  of  Aj _j*A|  ^  becomes  much 
better  as  I  —  •  then  the  approximation 
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-  <1 

l  fn  (n)  *  ^  fn  (n)  +  ^  f  (I)  +  ^  f(x)dx 
n  *  1  n s  1  I 

is  quite  good  for  t  large. 

In  the  particular  case  of  interest,  the  function  f  (n)  goes  to  zero  so  rapidly 


that  *  *  5  is  quite  adequate.  It  is  still  necessary,  however,  to  evaluate  \  f(x)dx 


which  in  this  case  is 


2  2  I 

e  ”a  x  2  2  2 

- x — dx.  The  transformation  y  =  a  (x  -  I  )  gives 


2  2  2  2  2  - 
P  ■»  x  .  -ate 
\  e  dx  k  \ 


e  L  dy 

,  A  2 .2,' 
(y  +  o  I  ) 


The  integral  term  is  calculated  by  Gaussian  quadrature  with  weight  function 

e  ^  on  the  interval  (0,-).  The  seven-point  quadrature  formula  is  used  as  the 
2  2  -3/2 

function  (y  +  a  t  )  goes  to  zero  quite  rapidly. 

This  approximation  was  used  in  calculating  dc/dt  for  the  given  analytic  sur¬ 
face  concentrations  used  in  checking  the  earlier  HYDRA  code  and  excellent 
agreement  was  obtained. 
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APPENDIX  II 


FORTRAN  LISTING  OF  HYDRA  II 


*  FORTRAN1  STAB t  DECK «L ST  IN 

*  INCODE  |BMF 
CHYDRA 

C  THE  NEW  HYDRA//THE  TWO  HEADED  MONSTER// 

COMMON  5 1 G  .  EPSK «  XM . DEN 1 « DENE « C  2 « CON « SRA T  .EPS 
DIMENSION  S1GC2) .XM<2) .EPSK<2>  *C0N<2> 

DIMENSION  XMATA(30«30) «BX(30)«XMAT(30*3Q) *  YMUE ( 1 1 > 

DIMENSION  XLAM!  101  •  1  1  )  .PSI  <2  •UliVI  11)«T(2  *11)*C(2  ill) 

1  .0(501)  *P<  2  .12). DC BAR  <  l  1  )  » FH (  2  «12).PI(2  .12). 

2  TP(lOl) • ESSN ( 1 1 ) .DOFT ( 1 1 ) «SEE( 101.11) .CAY1 < 11). 

3  CAY2< 1 1 ) .DAFT ( 101.11  )  »CEE<4« 1  I ) . S A ( 4 ) . SB ( 4 ) 

ECiUl  VALENCE  (SEE!  1  )  .XMATA(  l  )  ) 

C  max  ZONES  is  10 

DATA  PlE.SA .SB.XMH.R/3. 1 4 1 5927  * -3.34 1 5 ♦ 38.6433 . -34 » 2639 . 

I  9.2821 ,-31298.2.23574. t . -6028. t 0 • .2. 0 16 . 8 . 3 1 4E7/ 

HCON* 1,1 5*1 .987*4. 186E10 
S I 0*8,4972 

EP$K*38. 

XM.x2.016 

ROMH.R/XMH 

DCBAR.o. 

C  GET  INPUT 

22  READ  <  5. 101 1NZ0NES.NTAB. INPT 

READ  (5. 100)XA,XB.DP. FD.DEPTH.EPS. FF.TMAX.DELTAT. 

1  TAMB.TZERO.FI . RHOS. CESS. AREA , XKL . (ESSN( J) .J*l .NZONES) 

READ  (5. 100)SIGI .EPSK I • XM I 

READ  (5.100)  (TP(J).O(J  )  *J*1.NTAB) 

C  DEMAX  LE  ?ERO  KEEPS  DELTAT  FIXED 

r.  WHEN  ROMANV  LESS  THAN  ZERO  KEEP  PERCENTAGE  OF  INLET  INERTS  CONSTANT 
READ(5. 100)  PIV1 .PIV2.0EX. ROMANV. DMAX.DMIN 
WRi'.E  (6.105) 

105  FORMAT! 12Hi HYDRA  INPUT) 

WRITE  (6.106) 

106  FORMAT (24H0PRESSURE  TABLE  (P.VS.T)) 

WRITE  <6.1 07 ) ( Q ( U  ) . TP ( J ) . J* 1 .NTAB ) 

107  FORMAT (1H  2E15.6) 

WRITE  (6.108) 

108  FORMAT! 15HC INERT  GAS  DATA) 

51*100. *P 1 V 1 /Q ( 1 ) 

IF(ROMANV.LT.O. )  81*100. *Fl 

WRITE  (6.1 09  JBI.SIGI. EPSK 1 . XM 1 

109  FORMAT ( 23H  FRACTION  OF  INERT  GAS*F 10.6. 7HPERCENT » 

1  2X . 7H , S I GMA  *F5 . 3 . 2X . 

1  7H.EPS/K*F4.1 «2X. I4H, MOLECULAR  WT *F6«  2 ) 

WRI TE ( 6. 127 )  PI VI .P1V2. QEX . ROMANV 
127  FORMA? ( 6H0P IV1*E15.6»2X.6H.PIV2*  E 1 5 .6 » 2X ♦ 3H . 0*  F7.5.2X. 

1  BH.VOLUMN*  E15.6) 

PI V2*P I V2*6. 894 7E4 
WRITE  (6.110) 

110  FORMAT ( 16H0T J ME  PARAMETERS) 

WRITE  (6* 11 1 ) TMAX. DELTAT. DMAX 

111  FORMAT ( 1  OH  MAX  TIME*E6.2.  4H«DT*F6»3.6H.DMAX*  E15.6) 

WRITE  (6.114.1 

112  FORMAT ( 32HOO 1 FFUS I  ON  COEFFICIENT  CONSTANTS! 

WRITE  (6.113) XA « X3 

113  FORMAT  (  3H  A*El5.7.2X.3H.B*Ei5. 

WRITE  (6.114) 

114  FORMAT (27H0REMA I NING  INPUT  QUANTITIES) 

WRITE  (6. 1 15)NZ0NES«DEPTH, EPS. TZERO. AREA 

U2  FORMAT! 1 7H0NUMBER  OF  ZONES* I  2 • 2X« 14HLENGTH  OF  BED-F5.2.2X. 

1  S4HV01D  FRACTI0N*F4.2.2X. 13MINITIAL  TEMP«F9.4 . 2X .  I  7HX-SFCT  I ONAi. 

2  AREA.-.F6.3) 


148 


WRITE  (6*  130lXKL«RH0S»CESS«0P.FD.FF 

1  JO  FORMAT ( 17H0HEAT  LOSS  CGEFFb  E 1 5. 7 .2X «SHRHOS«  Eli. 7/ 

1  6M0CESS*  E15.7.2X.3HOP-  E 1 5.7 .2X . JHFD»  £ 15.7.2X*3HFF»  Elb.71 
WRITE  (6.1  16M  J.ESSNI  J)  .  J«1  .NZONES) 

116  FORMAT ( 3H  S< I 2 • 2H I s£ 1 3.7 1 

C  INITIALIZATION 

S1G( 2 ) bS IGl **2 

EPSK(2)bEPS<I 

XM<2)mxMI 

RAT«XM<  2 1/XMC 1 ) 

SRATbSqRT  (RATI 
OENIbSQRT  < 8.*< 1 .+RAT1 1 
DEN2«SQRT  <8.*( 1.+1./RATJ1 
0040  I  a  1 .2 

40  CONI  I )b2.665E-5*SGRT  {XM(I1)/SIGCI1 

C2« 15. *1.987/1 4. *XM(2» )*4.184E7 
CESS»CESS*4 . 1 86E7 
XKL*XKL*1 .8E7 
TZEROb ( TZERO+460 .1/1.8 
TAMB  BjTAMd  *460 .1/1*8 

xz*nzones 
NZONESbNZONES+1 
YMUE (NZONES IbO. 

OELXbOEPTH/XZ 
A«.5*DP*FD 
PIA»PIE**2/A**2 
NZMIbnzONES-1 
NZM2“2*NZM1 
NHZM1b3*NZM1 
NLbNZM2+1 
NZT1*NZ0NES+1 
DO  200  J* 1 . NNZM 1 
DO  200  KbI.NHZMI 
200  XMAT (J.K)bO. 

401  EPSM!aEPS-l. 

EPDTbEpS/DELTAT 

D*(DP*fF*EPS/EPSM1 )**2*EPS/( 150.*DELX1 

DD-D/OELX 

P( 1 . 1 }bG< 1 1 

PI < 1 .1 »«0( 1 1#(P1V1/Q( 1 1 1**0EX 
IF<R0MANV.LT.0.1  PI < 1 . 1 JBF1*P( 1 » 1 l 
P(  1  « 1 JbP( 1. 1 J*6.8947E4 
PI  ( 1 . 1 >bPI (l.i ) *6. 8947E4 
PH(1.1)bP(1.11-PI<1,1) 

VIN«0. 

ROC ANb- AREA /ROMANY 
JS«2 

CALL  FAZEDA (PH( 1 , 1 1 .TZERO.CZERO) 

WRITE <6* 1321  TZERO.CZERO 

132  FORMAT ( 14M0 INITIAL  TEMPb  E  15  »  7 » 2X.22HIN1 T I AL  CONCENTRAT 1 ONb 

1  E15.7} 

EPDXbepS»OELX 

AREXbapEA#DELX 

AREYbAREA/DELX 

ADTMbrhOS/ ( 1 C  008*CZERO >9 1.221*1. 008 
EPSRMbePSM 1 *ROMH  * ADTM*6 • /P I E* *2 
DCCONb-EPSRM*P I E  **2/6 • 

ASCONb-DELX *EPSM 1 * ADTM* ARE A 
CLASSbADTM/(RH0S*CESS1 
CLESSb-1 ./(RH0S*CESS*EPSM1 1 
DHPCONbI./  ALOG( 10.1 


DO  1  Ms  1 tNZONES 
DAFT ( 2  «M ) »0 • 

ESSN I M ) sESSN ( M )  *  1 t E7 
XL AM ( 1 ,M)*0. 

PSII ltM>-CZERO 
P( 1 «M)sP( 1*1) 

PHI  1 »M ) sPHI  1«1) 

Pill tM J-PI I  1 « 1 ) 

T I  1  tMjsTZERO 
Cl  1  *M J  mCZERO 
I  CONTINUE 

PHI  1 .N2T1 
PIll.NZTl 
PI  1  «NZT1 
TIMEsOELTAT 
C  START  TIME  COUNTER 
I  *2 

CSTARsq. 

INOW-2 

23  JCAP*InOW-1 

DO  26  J* JS« NTAB 
IS- J 

IFI  <TP| J)-TIME)*ITPI J-1)-T!ME>  )  30t30t28 
28  CONTINUE 

C  IF  EXIT  HERE  WILL  EXTRAPOLATE 
30  PI  I t 1)«QI IS-i  >  +  tOt IS)-Q(IS-1 )  )  *1  T IME—TPI IS-1 ) )/ 


S-PHI 1 «NZ0NES) 
) -PI  II tNZONES) 
)«  PI  1 tNZONES) 


I  TP  I IS)-TP< IS-) 


PI  I t 1 j*PI 1 1 l ) *6# 8947E4 

PI V2-P I V2+ROCAN*V I N*P I  I  1-1 1 1 )#DELTAT 
IF I  PI V2«LT • 0* )  PIV2-0* 

PI(Itl)-P(I t) )*IPIV2/PIItl)I**QEX 
IF  IROMANVtLT  «0t )  PH  Itl)«FI*PI  It  I) 

PHI  1 1 1 )*Pt I 1 1 ) -P I  I  I t 1 ) 

LITER-0 

DO  3  J»2 tNZONES 
PHI  I . JJ-PHI  1-1 t J) 

PH  ItJ)»PH  I-ltJ) 

PSI  I  1  t,J)  -PS  1 1  l-l tJ> 

PI  1 1 J)«PH| I tj)+PI( 1 1 J) 

3  TII tJ>-TII-ItJ) 

TIItl)sTII-ltl) 

PHIItNZTl  ) -PHI  I  — I tNZT 1  > 

PIlItNZTI  ) -P I  I  1 - 1 tNZT )  > 

PIItN2Ti  .'*  PI  I  - 1  tNZT  1  ) 

20  DO  6  M-2 tNZONES 

OOFT I M ) -XA/EXP  IXB/I lt987#TI ItM)  )  ) 

XLAMI iNOWtM  )-T IME*DOFT IM ) 

DO  SO  1 1 -It  INGW 

ALPS-P I  A* I XL AM  I  I  NOW t M ) -XL AM I  I  I tM>  ) 

CALL  SERVAL  I  ALPS t SEE  I  I  I tM ) ) 

50  CONTINUE 

6  CONTINUE 

DO  31  m-2 tNZONES 
DO  37  J— 1 t A 

57  CEEI JtM)-SA I J)+SB( J)/TI I tM) 

CSMs { SEE  I INOW  t M ) -SEE  I JCAP  t M )  )/|XLAM{ I NOW t M ) -XLAM {  ICAPtM) ) 

CAY2 I M ) —I  DOFT I M ) *CSM#EPSRM 
CAYHM)--PSI  I  t-l  tM)#CSM 

52  DO  33  j«2t*JCAP 

53  CAY 1 1M)-CAY  H  M  >  + I  SEE I J  tM ) -SEE I J-l tM> )#DAFT( JtM) 

51  CAYHM)—  D0FTIM)*CAYHM)*EPSRM 

NCuT-0  150 


ISUR*0 

I F  <  F  1  )  400.400.402 
C  ONLY  DO  WHEN  CHANGE  T 
400  DO  202  Jal.NZMl 

201  XMAT ( J. J ) > 1 . 

jj«j+nzmi 

202  XMAT ( JJ. JJ ) *CAY2 ( J+ 1  ) 

GO  TO  8 

402  DO  300  J*NL.NHZM1 

JJ* J-NZM2+1 

300  XMAT ( J • J ) *C Ay2 ( J J ) 

8  ABSORD-0. 

00177  M-1.NZM1 
!F(M.SE.t»  GO  TO  214 
TEMP«T( I . 1 ) 

PRESH*PH< 1.1) 

PRES 1 *p 1 (1.1) 

GO  TO  21S 

214  TEMP* ,5# (T ( I »M)+T( 1 «M+1 ) ) 

PRESH* ,5* (PH( I .M ) +PH ( l «M+l ) ) 

PRES  I  * ,5* ( P I ( I «M)+Plll «M  + 1 ) ) 

215  CALL  PROPiTEMP. PRESH. PRESI .XMUM.XKBED* 1SUR. IVIS) 

V  <  M )  .DD/XMUM 

IF(M-i)  2 10.210.211 

210  YMUE< 1 )-XMUM*DELX*T( I . 1 ) 

V*2.*V 

GO  TO  177 

211  YMUE(M)*XMUM*DELX*(T< I »M)+T< 1 «M+1 ) ) 

177  YMUE<M)*2.*D/YMUE(M> 

V(NZ0NES)*0. 

IF ( F I >59.59.60 
59  CONTINUE 

DO  213  J*1.NZM1 
JJ* J+NZM1 
8X< J)*6» 

DO  212  K*1.4 

212  BX ( J ) »0X ( J ) +CEE  <  K . J+ 1 >*PSI < l .J+l )**(<- 1 > 

BX( J) *-1.01325*10. **B*< J) 

XMAT ( J. JJ) *0, 

DO  204  K*2.4 

XK*K— I 

204  XMAT ( J. JJ ) *XMAT ( J. J J ) +XX+CEE ( K . J+l )*PSI C I • J+l )**CK-2) 

XM AT ( J , J J ) - BX  <  J ) *XMAT ( J . J J ) /DHPCON 

8X1 J)*-BX< J)-PH< I. J+l ) 

213  CONTINUE 

XMAT<NZONES* 1 )«(PH( I . 2 ) *YMUE<2 ) +PMC  I. 1 )*YMUE( 1 )+EPDt/T( I .2) ) 

DO  205  K-2.NZM1 
JD«NZM1+K 

C  JO  IS  ROW  NUMBER 

XMAT  (  JO « K  >  *  < PH (  I  .K+l  )  *  <  YMUE  ( K+  1  )  ♦  YMUE  ( K )  l+EPOT/TI  I  »K+l  )  ) 

XMAT( JO-l .K )«-YMUE(K)*PH( I «K+l > 

205  XMAT( JO.K-1 ) ■— YMUE  <  K ) *PH ( I .K) 

DO  206  J*1  .NZM2 

DO  206  K* 1 . NZM2 

206  XMAT  A ( J.K ) *XMAT <  J.K ) 

DO  209  Jal.NZMl 

J J*NZM 1 ♦  J 

BX ( JJ ) *PS I < I . J+l )«CAY2( J+l ) 

l  -( .5*:yMUS< J+l >*IPMt I .J+2)«42-Ph< 1  .  J+l )*«2)> 

1  YMUE ( J)*(PH ( I . J+l )«*2-PH( I .JI++2) l-CAYl l J+l  )-EPOV*IPM< I . J*l )/ 

2  T< I . J+l )-PH( I-l. J+l )/TI I-l.J+1 )) ) 

BX( JJ)*-3X( JJ) 
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209 


CONTINUE 

BX(NZONtS)«PSI I I «2>*CAY2(2) 

1  -I .5*YMUE(2)*IPH( l.3)**2-PH( l .21**2)- 

1  YMUE I  1 ) *PHI  1  til* (PHI  1 «  2 ) —PH ( l  »  1  > )-CAY 1  1 2 ) -EPOT* ( Ph I  I «2)/T{ l «c> 

2  —PH 4  I— 1 »2)/T( |-1 *2) ) ) 

BXINZM1+1 ) ■ — BXI NZM l  +  l ) 

CALL  MATINVIXMATA.NZM2.BX) 

C0NE*0, 

00  207  J-1.NZM1 
PH( I • J+l J*PH( I • J+l )+BX< J) 

CONtaCONE+ABS  (faXI JI/PHI I . J+ 1) ) 

JJaJ+NZMl 

PSI < I » J+l >aPSt I I .J+1 )+BX( JJ) 

CONE*CQNE+ABS  (BX( JJI/PSI ( I « J+1 ) ) 

207  CONTINUE 

PHdtNZTl  )  a  PH  I  I  . NZONCS  ) 

PI  I *NZTI)«PH( I «NZT 1 ) 

NCUT.NCUT+1 
IFINCUT-IO)  66.66.62 
66  IFlCONE-l.E-9)  62.62.S9 

60  CONTINUE 

C  INERTS 

00  30S  J*I.NZMt 
JJaJ+NZMl 
JJ J* JJ+NZM1 

XMAT I J«JJ)aQHPC ON/PHI  I «J+l ) 

XMATI J.JJJ) «0. 

DO  302  K»2.4 
XKaK-) 

302  XMATI J,JJJ)a-XK*CEE(K.J+l >*PS1 I  I .J+1 > **(K-2 ) +XMAT I J, JJJ) 

BXI J)»— OMPCON*ALOO(PH( I .j+1 I/I .0132SE6) 

Kaj+l 

BXIKK ) •— EPOT* (PHI  I .K)/T( 1 .KI-PHI t-1 .KI/Tl 1-1 .X ) ) —CAY l  IK  )  — 

1  CAY2I K  )  *PSI  I  1  .K I +.9*  I  YMUE  HO*  I  PH  I  I  .K+l  I+PHI  I  «K>  )*(P(  1  .X+l  )- 

2  P(1.«))-VMUEIJ)*IPh( l.K I+PHI J.JI )*IP( I.K)-P(l.J)I) 

00  306  Ka 1*4 

300  BXI  JlaBXI  JI+CEEIK.J+1  )*PSI  1 1.  J+1  )**IK-l ) 

BXI JJ) ■— EPOT* I PI  1 1 • J+I )/T( I . J+1 ) -PI  1 1 -I. J+1 )/TI I-i * J+l j )  +  .  j*l 

1  YMUE I J+1 )*(PI I  1 • J+2 ) +PI I  I « J+1 ) )*IP( 1.J+2I-PI I .J+1 1 1 -YMUE I J>*( 

2  PI  I  I. J+1 1+Pt 1 1 «J))*(P( 1. J+l I -PI  I .J) ) 1 

309  CONTINUE 

00  303  JaZ.NZMl 
JJaJ-l+NZM| 

J+I 

KKaJJ+NZMl 

XMATI JJ. Jla— ,9+YMUEI J)*(P( I.Kl+PI I  I tKI-PMI I «J| 1 

XMATI J J+l. JlaEPOT/TI I .X 1 -.9*1 YMUE IK )* I PHI  I .X+| 1 -PI (l.K 1— PI  l.K)) 
1  —YMUE! J)*IPI I  I »X)+P( I .Kl-PHI 1 • J))) 

XMATI JJ+l.J-1 I  a-. 9* YMUE I J)*(PI  I.JI+PK I.Jl-PHI l.K) ) 

XMATI JJ. JJ+1 ) a-. 9* YMUE I J)*IP1 (l.K) +p I  I  I • J ) ) 

XMATI J j+1 .JJ+1 ) a .9* I YMUE (K)*«PI( 1 »K  +  ! i+PI I  I .X) )+YMUE( J)*l 

1  PI  I  1 .K)+PI( I • J) > ) 

XMATI jj+1. JJ)aXM ATI JJ. JJ+1 ) 

XMAT IKK* J)a  — .9*YMUE I J )• I PH (l.K) +PH I  I . J ) ) 

XMATIkk+1 . J ) a «9* I YMUE I X ) • I PHI I.K+I )+PH(I »K| I+YMUEI J)*(PH( |.K)+ 

1  PH(I.Ji)) 

XMATlKK+l . J— 1 ) a XMAT IKK. J) 

XMATIKK. JJ+1 )a-*9«YMUE( J)*(PH( I .X )+P( I *K)-P| I  1 1 J) ) 

XMATIKK+1 (JJ+1 laCPOT/TI 1 «K ) -.9* I YMUE I X ) • I P I  I  1 .K+l )-PI | ,K)-PH( 1 . 
I  K) | — YMUE! J) *(PH( I «K ) +PI 1 »X )— PI  I  I « J) ) ) 

XMAT(XK+l.JJia.D*VMUE( J)*(P|( I .X )-#( I . J) -PHI  I .J) ) 

303  CONTINUE 
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XKAT  INj0NeS.tt«e<»0T/T(l*2>  +  *S»  2,<MP<,'2,+P,'t'*,“PH,l'31M 

YMUE< 1 J*PH(  I il)*IPII »2)“PI I • *>  * 

DO  304  J* 1 .NHZP1 

00  304  K-l.NMZMl 

XMATA ( J.K ) "XMAT <  J»K  J 

CALL  MATINVIXMATA.NK2P1.8X) 

CONE*0. 

DO  307  J*l.NWl 

JJ« J4N2M1 

Pi(i.j+n*Pi(iiJ+ii  tbx  * J  * 

PH (  I . J4 1  )«PH(  I  « J4l >*BXt JJl 
KK«JJ4NZH1 

8X(KK)/PS1(  »*•»♦»>> 

).lt.0.i  Pi«i.J*n-PMi.i> 

P  (  1  .  J-f  1  )  *PH  C  1  .  J4  1  )  4P  1  <  1  •  1  - 

CONTINUE 

PHI I .NZT1 S»PHl 1 .N20NES) 

PI ( 1 .NZT 1 )“Pl < l.NZONES) 

PIl.NZTD*  P(l.NZONES) 
nCUT»nCUT4J 
IF<NCuT-lO 1300.300.62 
IF  (  COnE-1  .E-5)6i:  .62.60 
CONTINUE 
PSIII.l)«0« 

OO  65  P»2*NZ0N£S 

OCB AB  m » r !  cl "  It  P 1 1c  A  Y*  <  P 1  «»$  »<  *  •  "  »  > /OCCON 
ABSOROa  ABSOBO>OCB AB ( P I 
CONTINUE 
VA0S*O» 

DO  39  J* l.NZONES 

VIJ).-  IP«I.JTI»  *PM.J>l»VIJ»*DO.X 
V ABS> APAX 1 C VABS*V( J) > 

CONTINUE 
VIN-Vt I ) 

LlTEB.L>TCn*l  _  . 

FLlNapM(l»l)«»V(  l  )/IBOPH*Tf  !>>)) 

ABSOeO-ABSORO^ASCON 

CNVBaO. 

I  Sun -I 

IF(INPT)  125.1*5. 120 
CONTINUE 

•BITE  (6.123) 

. . 

FOB** AT  (  1H  3C15.7) 

CONTINUE 

DO  10  paB.NZOACS 

C(l,P).C( |-l.P)^OCBAB(  P)40ELTAT 
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DHN«HC0N*(-3l  *2982+C( I .M ) # (23.5742-C ( I <M) *6.028) ) 

CAUL  PP0P(T( I  «  M )  « PH (I.M)«P|(|«M) .XMUM. XKBEO. I SUR . !ViS> 
IF1XK0ED.LE.0. >  GO  TO  22 
IF(M.EO.NZONES)  GO  TO  15 

tl  XLN- < XKL* ( T < 1 .M)-TAMB)-XKBEO*AREY*(T( I.M*1)*T( l.M-l )-2.*T< 1 «M) )) 

XLN*XLN/AREX 
GO  TO  13 

15  MlaM-l 

16  XLN» ( XXL* I T I  I .M)-TAMQ)-XKBEO*AR£Y*(T( I .Ml )-T( l «M) ) l/AREX 

13  DTHtTAa-OHN*OCBAR(  M)*CLASS*(ESSN(M)-XLN)*CLESS 

V(M>  aTI 1-1 .M)*OElTAT*DTHETA 
XSV«AB*  (l.-ABS  (T(l.M)  /V  (  Mil) 

IF( INPT.EO.O)  00  TO  IB 
WRITCl6«13>>  DHN.XuN.DTHETA. XKBEO 
131  FORMAT (5M  OHNaE 15.7* 2X  *4HXLNa(  I5.7.2X. 7H0THETA.E 13. 7 « 

1  2X • 6HXKBED* E 1 5 • 7 } 

18  CNVR«XSV*CNVR 

10  CONTINUE 

IF(INPt)  120.120.121 

121  WRITE  (6.1 22 ) FL I N .  ABSORO • L 1  TER • NCUT . CONE  »  CNVR 

1 22  FORMAT ( 6HORHOV»E 1 5. 7. 2X • 9MOCBARSUM-E 15. 7.2X . 6HL 1  TER • 1 2  *  2X . 

1  SHNCUTa 12. 2X.5HC0N£aE 15.7. 2X«5HCNVRaE 15.7) 

12C  CONTINUE 

00  2  Ma2.NZ0NES 
2  T(l.N).VIM) 

T ( I  .  1 ) »T (1.2) 

C( I  •  1 )aC( 1 .2) 

IF C LITER-301  33.33.19 
3?  IF(CNVR-l.E-S)  19.19.20 

19  CONTINUE 

03  55  M*2.NZ0NE$ 

55  OAFT  (  InOW.M  )  a  (PS  I  (  1  .Ml— PS  I  (  I  *»l  »M  )  )/{ XLAM (  I  NOW  »M  )  —XLAM!  JCAP  «M  )  ) 

CST ARaCSTAR. ABSORO *OELT AT 
WRITE  (6.1 17) TIMC.CSTAR 

1 17  FORMAT  C6M0T lMEaF7.3.2X.6HCSTARaE 15.7) 

WRITE  (6.118) 

1 ia  FORMAT ( 1M07X.2HPH. 14X.2HPI • 15X. I HP. 1 2X . 5MTME  T A . 12X.4HCBAR) 

WRITE  (6.1 19) (PHI I.J)«PI(I*J).P(l*J).T(|«U).C(|tJ)t 

1  J-1.NZ0NES) 

119  FORMAT! 1H  5E16.B) 

IF ( ( VABS.EO .0. ) .OR. ( EPOX.EQ.O • ) 1  GO  TO  3) 

DINT  aEPOX/VASS 

34  IF(OINT.GT.OCLTAT)  GO  TO  32 

IF  ((.S40ELTATl.LT.OMIN)  GO  TO  3. 

0ELTATb.540ELTAT 
GO  TO  34 

32  IFIDCLTAT.GE.OMAX)  GO  TO  31 

IFI2.40ELTAT.GT.OMAXI  GO  TO  31 
|F(2.*oELTAT.GE.0INT)  GO  TO  31 
OELTATatc*OELTAT 
31  CONTINUE 

TIME*T|MC*OELTAT 
| NOW a | NOW* 1 

C  SET  EVERYTHING  BACK 

001305  Mai.NZONES 
Cl l.M)aC(2.M) 

T( 1 .M).T(2.M) 

P( 1 «M)aP(2«M) 

PH(  |  .N ) aPHI  2  *M ) 

PI ( 1 »M) aP| ( 2.M) 

PS I ( I .MiaPSI (2.M) 

CONTINUE 
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f*  <"»  non  non  non  non  non 


|F<  INOw.LE.  ton  GO  TO  25 
I  NOW *52 
00  27  j-1.51 
JJ« J*50 

00  27  k-I.NZONES 
XL AM { J(K) bXLAMI J  J  <K  ) 

27  OAFT (  j, K  1 «0AFT ( JJ«K ) 

29  !F{ T IM£-TMAX)  23.23.22 

100  FORMAT <661 2 .5) 

101  FORMAT ( 1419) 

ENO 

*  FORTRAN  sTAB.DECK.LSTIN 

•  INCODE  jsmf 

CCOMAT 

MATRIX  INVERSION  WITH  ACCOMPANYING  SOLUTION  OF  LINEAR  EQUATIONS 
SUBROUTINE  MAT1NVI A.N.B) 

DIMENSION  IP  I VOT ( 30 ) • A  <  30 • 30 ) • B( 30  I 
EQUIVALENCE  (AMAx .PIVOT) 

INITIALIZATION 

19  00  20  J-l.N 

20  IPIvOT (  J)  aO 

30  00  990  I  *  1  .N 

SEARCH  FOP  PIVOT  ELEMENT 

*0  AMAX«0.0 

*9  00  109  J.i.N 

90  IF  (IPlVOT(J)-l)  60.  109.  60 

60  00  IOC  K«1.N 

70  IF  « IPIVOTCKl-l )  80.  100.  7*0 

80  IF  (ABS  IAMAXI-A8S  (AIJ.K)))  89.  100<  100 

89  IROW* J 

90  l COLUM *K 

95  AMAX»A<  J.K) 

100  CONTINUE 
105  CONTINUE 

IP 1 VOT ( I COLUM I  *  1 

INTERCHANGE  ROWS  TO  PUT  PIVOT  ELEMENT  ON  DIAGONAL 

130  IF  I IROw-ICOLUM)  1*0.  310.  1*0 
toO  00  200  L* I «N 
160  SwAP*A< IROW.L) 

170  A| IROW.L )«AI 1 COLUM. L) 

200  A(  I  COLUM.  t„)>  SWAP 
SWAP-BI IROW) 

81 IROW) *8C  I COLUM  > 

B( l COLUM  )*SWAP 

DIVIDE  PIVOT  ROW  BY  PIVOT  ELEMENT 
110  OO  390  L«t.N 

150  A(  ICOLUM.L)*A(  I  COUP* »L  1 /P I  VOT 

B  (  I  COUP*  )*Bl  1  COLUM  > /PI VOT 

REDUCE  NON-PIVOT  ROWS 
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380  00  500  L 1 .  1  «N 

39C  IF(Ll-ICOLUM)  400*  550.  400 

400  T.A(L1 . ICOCUM) 

IF(T,EQ.o.>  GO  TO  550 
430  00  430  L-l.N 

490  A1lI«L)*A(L1.L1-A<  ICOluM.LMT 
9<Li ) *8(L 1 1 —B  ( ICOLUM  >*T 
S50  CONTINUE 
740  RETURN 
ENO 

*  FORTRAN  STA8.0ECK.LSTIN 

•  INCOOE  JBMF 
CPROP 

SUBROUTINE  PROPITEMP.PPH.PPI.XMUM.XKBED. ISW.IVIS) 

C  SUBROUTINE  PROPCA  CALCULATES  VISCOSITY  AnO  THERMAL  CONDUCTIVITY 
C  THE  INITIALIZATION  values  NEEDED  are  (I)  MOLECULAR  WEIGHTS 
C  (2)  COLLISION  DIAMETERS  ANO  <3>  K/EPSILON  TO  CALCULATE  OMEGA 
COMMON  SIGtEPSK.XM.DENI ,D£N2. C2 « CON. SRAT  .EPS 

DIMENSION  SIG(2>.EPSK(2).XM(2).PHI(4).C0N(2).ZZU6>.  VAR( 

DIMENSION  0MEGC2) .XMU(2) .XXI 2 I «XK(2) 

DATA  EXPN.C1 .PHI ( I ).PHI (4) . (ZZI I  I. 1*1 « 16 >/• 33333333. 2.4837b. 

1  I • . 1 • .473. .2.36E6.973..2.3E6.673..2. 1SE6.773. . 2.01E6.873. • 

2  1 .BE 6 <973. • t .S1E6. 1073. .1.1 3E6* 1 173. «6t 39E5/ 

3  P«PPH*PP I 

SOT "SORT  I  TEMP) 

DO  4  1.1.2 

4  vAR< I j.TEMP/EPSKt 1 ) 

OMEGC  1  )>3.019/(ALOG(  VAR(  1  )  )0.772> 

IF(VAR(2)-3.4)S«9.6 

9  0MEG(2)>3.36/(AL0G(VAR<2) >+2.125) 

GO  TO  T 

6  OMEGl 2  )*9.0  13/  ( ALOGI  VAR(  2 1  I -*3*  772 1 

7  CO  8  1.1.2 

8  XMi( I i.CONI I ?*SOT/OMEG( I  I 
XMXM.XMUf  t  S/XMUI2) 

SOT. SORT  (XMXM4SRATI 
PHI <2|«! 1 «+SOT)4*2/OEN2 
PH 1 ( 3 ) . ( l .♦ 1 ./SOT ) **2/OEN 1 

XK(  1  >■(  (6.624-. 00081  *TEMP+Cl  >+XMU<  1  » /XK {  l  )  144.184E7 
XK(2).  C2*XMU(2) 

14  CONTINUE 

xkm.o. 

XMUM«0. 

XX ( I ) .PPH/P 
XX(2).pPI/P 
DO  9  1.1.2 

YN*XX(  t  jy'IXXl  1  )#PH  1 12*1-1  ).XX(  2  1  *PM!  12*1  )  ) 

XMUM.XMUM-.YN4XMU1  I  1 

9  XKM«XKW+YN*X1CC  I  ) 

IFMSw.EO.O)  RETURN 

K» 1 TeMP-473. )«l.E-2+3. 

IFCK.GE.8)  K.7 
IFIK.LT.11  K. 1 
C  BETWEEN  KK  ANO  KK+2 
KK»2*K-1 

XKSOL*ZZ(KK«l  >+<T£NP>22l<X)  1*1  ZZtKK»3)<~Z2fKK.t  »  »•!  .E-2 
RT»CPS*« XKS0L-X*M|/HKM*4€XMN 
XKBCO.EPS* *  WCM-XKSOL I ♦XXSOL 

10  XM3>XK8C0«*eXPN 

F  •  1 X*S0C-XK8€0 I  ✓XMS-RT 
FP. - <XKS0L*£ • 4XKBE0 » ✓ 1 3c •a*3«XKe€0 I 
0€L“-F/FP 
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XKBED"XKBED*DEL 
IF(XKBEO) t 1*11*12 

11  WRITE  (6.100) 

100  FORMAT ( 50H  Btu  CONDUCTIVITY  CALCULATED  ZERO  OR  NEGATIVE  > 

WRITE  (6*101 1XKBED«RT ,XKSOL*XKM,F«FP«OEL 

101  FORMAT ( 1P7E 12*4 ) 

RETURN 

12  IF ( ABS  ( DEL/XKBEC 1-1 »E-5 113*13*10 

13  RETURN 
END 

l  FORTRAN  stab*deck*lstin 

*  INC ODE  J6MF 
CSERv 

SUBROUTINE  SERVAL  (ALPS* ANSI 

C  THIS  SUBROUTINE  EVALUATES  The  SUM  FROM  l  TO  INFINITY  OF  THE  SERIES 
C  EXP(-ALPS*N#*2 )/N*«2 
DIMENSION  XJ(«1 
OATA  PfE/3* 1*15927/ 

DATA  XlS«XJ/25.«1..4..9.*16./ 

IF ( ALPS.GT.O. )  GO  TO  1 
ANS«PIE**2/6. 

return 

1  ALLS»ALPS*XLS 

EXS>.5»EXP( -ALLS ) 

ANS«  EXS/XLS 
DO  3  J«1.4 

3  mNS>AnS4ExP  (-ALPS  *XJ(JI)/XJIJI 

C  GAUSS  CALCULATES  THE  INTEGRAL  PART  OF  THE  S’JM  APROXIMATION 

4  CALL  GAUSS (PSUM.ALuS) 

ANS*ANS*  sort  ( ALPS)*EXS*PSUM 
RETURN 

ENO 

*  FORTRAN  ST AB. DECK , LST I N 

%  tNCOOE  IBMF 

CGAUSS 

SUBROUTINE  GAUSS (PSuM, ALLS) 

DIMENSION  A(7)«X(7) 

DATA  A  *X/. 4 093 1893*  *42183128* • 1 4712639* • 206333 14E—1 . 

1  . I0740101E-2. . 1S869464E-4* .31 703155C-7*. 19304368* I .0266649* 

2  2.5678767,4.9003331 .8.1821336. 12.73418* 19.393728/ 

P5UM-0. 

CO  1  J» 1*7 
V»X( J)+ALLS 
F«V*SORT  ( Y | 

1  PSUM.pjuM+A  (  Jl/F 

RETURN 
END 

*  FORTRAN  STAB.OECK.LSTIN 

S  INCOOE  IBMF 

CFAZCDA 

SUBROUTINE  FAZEDA(OY1«S.THETA.C) 

DIMENSION  A(3) *B(3) *0(4) 

DATA  A,B*0( 4) « THC3/-3.84 15 *38. 6433* -34.2639*-3 1298.2* 23374. 1 • 
1  -6028., 9.2821 *27.8463/ 

PRCSHaOvNES 
C  » 1 »  73 

XLEFT  •  ALOG(PRCSH)  / 2. 302383 
DO  1  Ja  1  •  3 

1  OUi*A(Jf*«(J)/TKTA 

TC2-2.*D(3l 

3  F "XLEFT- ( D( I |4C4lD(2)4C»(0(3)4C»O(4) ) » I 

FP"-(0(2)4C4(TC2*C4TnC3» ) 
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OELC>—r/FP 

C>C*OELC 

,F‘*8S  (OELC/O-l  .£-6 1998.998. 3 
EJOlL  >9*9656  -  11494.3  /  THETA 

IE  (XLEPT  -  cOUlL)  997.997.996 
*n,TE  <6. 1011 IPRESH. THETA 

PORN AT ( 29*0  OUTSIDE  DELTA  REGION.  P 

OPF 10.3.  an  KELVIN  ) 

RETURN 
Cl 
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